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ABSTRACT 

The use of computer technology to control large pointing 
systems can significantly improve performance and reduce human work 
load. The goal of this thesis was to design software for an 
inexpensive, yet accurate and efficient control system for the 36- 
inch reflecting telescope owned and operated by the Monterey 
Institute for Research in Astronomy. Within this thesis, a 
computer program is developed to automatically move the telescope 
to a set of celestial coordinates and track with an accuracy of 
one-tenth of an arc second for five minutes within 75° of the 
zenith. Set times are anticipated to be between four and thirty 
seconds. Corrections are made to celestial coordinates to account 
for precession, nutation, aberration and atmospheric refraction 
effects. The user is provided an interface to the computer-based 
system that allows storage and editing of 100 star positions, 
editing of the system parameters and display of the telescope's 
status. Manual control of the telescope is also permitted at any 
time. Safety of the telescope structure is the primary concern of 


system software. 
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I. INTRODUCTION 

Digital computing equipment, in the form of personal 
Sempucers, Can provide a powerful tool to control large 
physical systems. This thesis explored the feasibility of 
Mpgrading the control system of one such systen, an 
astronomical telescope owned by the Monterey Institute for 
Research in Astronomy (MIRA). A computer program was written 
to control the telescope's movements for star tracking to 
several arc seconds degree of accuracy. A tracking rate error 
of less than a tenth of an arc second degree of accuracy for 
five minutes within 75° of the zenith was’ required. 
Additionally, requirements to interface a personal computer to 
the existing hardware were determined. Priority was given 


throughout the design process to the safety of the telescope 


Steemucture. 


A. STATEMENT OF PROBLEM 
The MIRA telescope was to be given the capability to meet 


the following design criteria: 


* Given a known star position, slew on command at the 
maximum acceptable slew rate to within one degree of the 
Star's calculated position. 


- Once within one degree of the calculated position, set via 
the right ascension and declination stepping motors to 
within one step of the calculated position. 


- The star's calculated position is required to be within 
ten arc seconds of the star's actual apparent position. 


- Manual adjustments to eliminate error between actual 
apparent and calculated positions are to be made via 
existing hand paddles or the computer keyboard. 


- Calculated corrections to a star's tabulated position are 
to include precession, nutation, aberration and refraction 
of light through the Earth's atmosphere. 


Once error has been sufficiently reduced, track the star's 
apparent position by compensating for rotation of the 
Earth and estimated refraction of the atmosphere. 

- A user interface is required, using standard astronomical 
vocabulary and symbols. This interface shall include the 
capability to read, write and edit a listing of one 
hundred star positions from a formatted data file. 

- The telescope's status is to be displayed at all times 
Guring telescope movement. The status display is to 
provide the following information: The telescope's 
position in celestial and terrestrial coordinates, local, 
Sidereal and universal times and the state of each motor. 

- The user interface shall also display adequate error and 
warning messages to inform the user of special 
circumstances. 

- Allow guiding of the telescope manually by hand paddles or 
keyboard at any time to eliminate or reduce subsequent 
errors. 

Use of a digital computer was chosen as the simplest 
design solution. Motor controller hardware, telescope 
position data and user choices were to be interfaced to this 


computer via the computer program developed by this thesis. 


B. THESIS SCOPE 
Requirements to upgrade existing MIRA hardware to 
encompass an automatic control system for the telescope fell 


into two categories; hardware and software. This thesis 


determined the hardware required to interface the existing 
motor sets and position determination devices with a personal 
computer and developed the software which that computer would 
use to control motor movement. hie “ehacking« accuracy 
requirement of one-tenth of an arc second required that the 
relation of the Earth's movement to any target star be 
precisely determined and accounted for. Refraction of light 
through a variable atmosphere and aberration effects due to 
the Earth's motion were also considered. 

The concept of a virtual and a physical telescope is used 
throughout this thesis. The virtual telescope is an 
abstraction in software upon which operations are performed 
without effecting the movement of physical objects. The 
virtual telescope represents ideal viewing conditions, a 
stable Earth with no atmosphere and no Sun to perturb the path 
of incoming light. Coordinates for a star, corrected only for 
precession, define the desired position of the virtual 
telescope for a given time. This position is easily 
calculated by an astronomer and is used to interface with the 
computer system. In reality, however, many sources of error 
other than precession must be considered, such as nutation, 
aberration and refraction of light through the Earth's 
atmosphere. The correction of a star's coordinates for all of 
these errors defines the required position of the physical 
telescope. An astronomer, in effect, positions the virtual 


telescope based on incomplete, yet readily obtainable data. 


The control system, guided by the computer, corrects the 
pointing instructions and positions the physical telescope. 
In effect, the virtual telescope is defined as pointing to a 
star's correct geometrical coordinates. The physical 
telescope must point elsewhere in order to see an object at 
those coordinates. 

Corrections to the virtual telescope include precession, 
nutation, aberration, and atmospheric refraction effects. 
Structural considerations such as telescope flexure, fork 
twist, droop, and encoder eccentricity also cause errors which 
were not considered within this thesis. These corrections 
require careful measurements of the physical telescope and 
must be made prior to system completion. 

To prevent distortion, the telescope's main mirror is only 
rear and side defined; it is not front constrained in its 
housing. Were the telescope to become inverted, this mirror 
would lose alignment and even possibly be at risk. Safety of 


this mirror was the overriding design concern. 


Ce BACKGROUND 
ts THE MONTEREY INSTITUTE FOR RESEARCH IN ASTRONOMY 
The Monterey Institute for Research in Astronomy is 
the only privately-owned professional astronomical observatory 
in the United States to be founded in this century. MIRA is 
a non-profit corporation dedicated to research and education 


in astronomy. It maintains a close association with two 


institutions of higher learning in Monterey, the Monterey 
Peninsula College and the U.S. Naval Postgraduate School 
(NPS). A memorandum of understanding exists between MIRA and 
NPS allowing the use of government-owned equipment on MIRA 
property and the teaching of NPS students by MIRA astronomers. 

MIRA owns and operates a 36-inch reflecting telescope 
located on Chews Ridge in the Los Padres National Forest. The 
Oliver Observing Station is located at approximately 5000 feet 
altitude above sea level on land owned by the U.S. Forest 
Service. Additionally, facilities in Monterey include 
administrative offices and an engineering and receiving 


building with computer, electronics and machining facilities. 


2 CURRENT HARDWARE CONFIGURATION 

The telescope at Chews Ridge is a 36-inch reflecting 
telescope. A 36-inch f/4 paraboloid mirror reflects incident 
light to a second, smaller cassegrain hyperboloid mirror 
mounted axially in front of it (see Figure 1). This second 
mirror directs light to the Guide and Acquisition Package 
(GAP). The telescope has a back focus of fourteen inches. 
The GAP is mounted axially, with a zero to six-inch back focus 
and may contain many instruments, including a spectrometer, a 
CCD camera, photographic equipment, and an optical eyepiece. 
The entire telescope is mounted eguatorially, which means that 


it rotates about an axis parallel to the Earth's axis. 


Approximate Celestial North 
(Polar Axis) 


: 


Secondary Mirror Primary Mirror 


ff 


_ Zenith 
(Vertical from Earth's Surface) 


L = Latitude 





Figure 1 - The MIRA 36" reflecting telescope with equatorial 
mount 


Rotation about the polar axis is then equivalent to changing 
the right ascension of the telescope (see Chapter II). 
Rotation about an axis orthogonal to the polar axis is also 
possible, which changes the telescope's declination, or 
celestial latitude. Two sets of motors are used to control 
telescope movement. Two one-quarter horse power slew motors, 
one each for right ascension and declination, are used for 
large telescope movements. For finer movements, such as 


tracking a star, two additional stepping motors are used. 


An optical encoder on each of the two axes detect 
movements of the telescope. This information, once referenced 
to a known location, is decoded by existing electronics and 
sent via an RS-232 link to provide the telescope's current 
(uncorrected) celestial coordinates in the form of declination 
and local hour angle. Right ascension must be calculated from 
the local hour angle by knowing the current sidereal time. 

3. CURRENT PROCEDURES 

Prior to this project, the four motors were controlled 

electro-mechanically by means of switches mounted on hand 


paddles. mreacking the Star, that 15, adjusting the raght 
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Figure 2 - Existing hardware prior to this thesis 
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ascension of the telescope at a rate approximately equal to 
the rotation rate of the Earth, 1s automatic. Guiding the 


telescope entails making adjustments for effects not 


associated solely with the Earth's rotation, such as 
atmospheric refraction. In this case, the virtual ana 
physical telescopes are identical. An astronomer wishing to 


observe a star would keep the telescope properly pointed by 
continual minute adjustments to the stepping motors. This 
process was long, tedious and tiring which lead directly to 
the need for an automatic system. Upon completion of the 
required hardware, astronomers of the future shall be able to 
sit in a warm-room environment while concentrating on an 


experiment's data instead of telescope operation. 


II. COORDINATE SYSTEMS 


A. CELESTIAL COORDINATES 

The celestial coordinate system used iS centered on the 
Earth and extends outward to the celestial sphere. The 
celestial sphere is a hypothetical spherical surface upon 
which all stars are said to reside. Celestial North is 
aligned with the Earth's North Pole. Declination corresponds 
to latitude on the celestial sphere and right ascension is 
analogous to terrestrial longitude. Where terrestrial 
longitude is referenced to the prime meridian at Greenwich, 
England, right ascension is referenced from the prime 
celestial meridian. The prime celestial meridian passes 
through the vernal equinox, sometimes called the First Point 
of Aries. Figures 3, 4 and 5 show these relations. 

Right ascension is defined as the angle, measured eastward 
along the celestial equator, between the vernal equinox and 
the great circle passing through the poles and the body. 
Declination is the angle, measured north- or southward along 
the meridian, from the celestial equator to the body. North 
declination is taken as positive, south as negative. 

An hour circle is a great circle passing through the 
celestial poles. An hour circle passing through an observer's 


zenith is also called a celestial meridian. Hour angles are 


Celestial 
Pole. 





Figure 3 - The celestial sphere 


angles meaSured east- or westward from a celestial meridian to 
a body. Local hour angles are referenced from an observer's 
meridian to a body, with westward angles being positive and 
eastward angles being negative. 

From any point on Earth, the point directly overhead is 
known as the zenith. Its opposite point, directly through the 
Earth is the nadir. An angle from the horizon to a body is 
known as the elevation. The zenith angle of a body is simply 
the angle from the zenith to the position of the body on the 
celestial sphere. It may be recognized as the co-elevation, 
or 90 degrees minus the elevation angle. The azimuth of a 


body is the angle from local North, measured eaStward in a 


10 


Celestial 
Pole 


_ 
Vernal Equinox 





Figure 4 - Right Ascension and Declination 


circle, to the body. A body lying due west of an observer 


would therefore have an azimuth angle of 270 degrees. 


B. CONVERSION OF CELESTIAL TO TERRESTRIAL COORDINATES 
Safety of the telescope was the primary design criterion. 
For alignment and safety reasons the telescope could not be 
pointed below an elevation of 15, even during testing of the 
new system. Therefore, care is taken to predict the position 
in which the telescope would be during a movement. To 
accomplish this, the celestial coordinates of the active star 
are converted to their respective terrestrial coordinates, in 


the form of azimuth and zenith angles. The two coordinate 


a 


Celestial 
Pole 





Figure 5S - Hour Angle and Declination 


systems are related through a direction cosine approach as 
follows: 

A rectangular coordinate system is determined on the 
Earth's surface, at the latitude (L) and longitude (y) of the 
observing site. This coordinate system is designated by the 
letter A. Local north is aligned with the number one axis, 
the number three axis points to the zenith and the number two 
axis completes a right-handed orthogonal set. Figure 6 shows 
this arrangement. 

The A-frame is related to an intermediate frame, the B- 


frame, as shown in Figure 6. The number three axis of the B- 


i2 


frame points to celestial North. This relationship is shown 
Pmairection cosine notation in Equation 2. 

The B-frame 1S related to another intermediate frame, 
denoted by the letter C, by Equation 3. The angle between the 
number one axes of frames B and C is equal to the longitude of 
the observing site, measured westerly from Greenwhich. 

The final, celestial, frame of reference is designated the 
D-frame. It is related to the C-frame by the Greenwhich Hour 
Angle of Aries, which is the angle between number one axes. 
The Greenwhich Hour Angle of Aries is denoted as GHAy. This 


rotation is shown in direction cosine notation in Equation 4. 


f 
me 


® A, . B. out of page 
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Figure 6 - Definition of the A and B Frames 


IES 


Sin(e0 —cos( 
DZGe— O a 0 =e] (2) 
Gos(lL) 0 SinGe 


(3) 
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Figure 7 - Azimuth (¢), Elevation (0) and Zenith Angle (z 
= 90° - @) related to rectangular coordinates 


The rectangular coordinate vectors must be represented in 


some more useful form. The unit terrestrial coordinate vector 
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D3 = Celestial North 


‘WN 


| 





Figure 8 - Right Ascension (a) and Declination (6) related 
to celestial rectangular coordinates 


that points to a star, a, may be related to azimuth and 
elevation (or zenith angle) as shown in Figure 7. Note that 
azimuth is measured clockwise from the A, axis and the zenith 


angle is equal to the co-elevation. Specifically, 


: (6) 
cos (8) cos (6) 
=|-cos (0) sin (o) 
sin (6) 


LES, 


Similarly, the celestial coordinate vector may be broken 
into right ascension and declination as shown in Figure 8 and 


below. 


(7) 


Solving for the terrestrial coordinates, assuming that a 
and 6 are known, and noting that °c* = bree) for an orthogonal 


set, 


a= "Ga (8) 


When solving for azimuth, quadrant dependence of the 
answer must be taken into account. This solution assumes that 
the distance between the center of the Earth and the observing 
Site is negligable when compared to the distance from the 
observing site to an observed body. This is true for all 
observations of stars, but would not be acceptable for viewing 


closer objects, such as Earth-orbiting satellites. 
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III. TIME 


A. UNIVERSAL COORDINATED TIME 

Universal coordinated time (UTC) is the time standard for 
the National Bureau of Standards and is based on the solar 
Gay . This standard of time is based on the atomic second, 
which is defined as "the duration of 9,192,631,770 periods of 
the radiation corresponding to the transition between two 
hyperfine levels of the ground state of the cesium atom 133" 
(Reference 2, Volume 1) , and iS measured by a cesium beam 
atomic clock. In the United States, radio stations WWV (Fort 
Collins, CO) and WWVH (Kauai, HI) broadcast time ticks based 
on UTC at one-second intervals on a continuous basis. The 
fine timing requirements inherent in astronomical tracking 
make it desirable to know UTC directly. This was accomplished 
by reception of the WWV radio signal and decoding into time 


information. 


B. TIME SYSTEM CONVERSION 

Universal coordinated time (UTC) may be easily obtained by 
iano, Signal. Earth's rotation makes any local coordinate 
system a moving frame of reference. Since Earth makes one 
revolution in one sidereal day, this leads naturally to the 
use of sidereal time for astronomical calculations. 


Conversion between the two systems is therefore necessary. 


aby 


The Greenwich mean sidereal time (GMST) is the sidereal 
time at Greenwich, England as referenced from the mean vernal 
equinox. Therefore, GMST is also the Greenwich hour angle of 
the mean equinox. At midnight (0 hour) on any given day GMST 
may be determined by the following relation: 


GMST at 0 hour=24110.548415+8640184 .812866° 


PAG) ES MOAR, DT 0 Ihe & 


(Julian Date-2451545.0) 


jes 
4 36525 
Di is defined as the interval of time Since noon Univemeew 
Time on 1 January, 2000, expressed in Julian centuries. See 


Reference 1, page Bé6. 

Greenwich apparent sidereal time (GAST) is found by 
applying a moderate correction to GMST to account forme 
difference between the true and mean equinoxes. 

In the everyday world, we deal with local time as a 
function of the time zone in which we are. The time at our 
le6Gation Vis, therefore, the same aS our neighbors. 
Astronomers and others concerned with precise time 
measurements, however, view local time much differently. Each 
longitude on Earth has a different hour angle from Greenwich. 
Local time is defined as Greenwich time plus east longitude. 
Therefore, you and your neighbor will only share the same 
local time if you are directly north or south of one another. 


This’ 1S known as focal mean sceolcam ee ine 
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Local mean sidereal time is then Greenwich mean sidereal 
time plus east longitude and local apparent sidereal time is 
Greenwich apparent sidereal time plus east longitude. 

From the above it can be easily seen that the local hour 
angle, important to many calculations regarding stellar 
positions, may be determined by subtracting the apparent right 


ascension from the local apparent sidereal time. 


C. TIME SYNCHRONIZATION 

The Earth revolves on its axis 360 degrees per sidereal 
day. At this rotation rate, one tenth of a arc-second would 
only take six-thousandths of a second. The accuracy required 
for the control system therefore mandates an extremely 
accurate method for determining and maintaining time 
Micormation. 

In the United States, two radio stations transmit time 
standard information on a 10 MHz carrier wave. The stations 
are operated and maintained by the National Institute of 
Standards and Technology, which was previously known as the 
National Bureau of Standards. Radio stations WWV, in Fort 
Collins, Colorado, and WWVH, on Kauai, Hawaii, broadcast these 
Signals which are synchronized to Universal Coordinated Time 
by way of the U.S. cesium atomic clock. 

The method chosen to allow the computer system time clock 
to be synchronized to this time standard was an electronic 


radio receiver/decoder made to fit an IBM-compatible personal 


le 


computer. This board is manufactured by Odetics Corporation's 
Precision Time Division, formerly known as Coordinated Time 
Link,. Cf wAnahein.. Calsict onoar It integrates an AM radio 
receiver with digital signal processing electronics which 
decode the WWV/WWVH signal into time and date information. 
The board is serviced by an external dipole antenna. An 
onboard battery and clock increment the computer's system 
clock between synchronizations. 

It should be noted that failure of the CTS-10 “onbe cmd 
battery will cause the computer's system clock to be reset to 
the CTS-10 factory default (afternoon of 19 December, 1991) 


until the system has re-synchronized. 


D. CORRECTING COORDINATES FOR TIME 

Many stars are listed in catalogs which give their 
celestial coordinates for a certain date and time. These 
reference dates and times are known as epochs. Epoch is Latin 
for 'time'. Epochs differ primarily due to precession and 
nutation in the interval separating them. The plane of the 
Earth's equator, and hence the celestial equator, varies with 
precession and nutation. Referencing stars' positions to a 
standard epoch in effect references them to an inertial frame 
Of sre terrence. The two most common epochs for celestial 
reference are B1950.0 and J2000.0. As one may expect, they 
refer to the years 1950 and 2000 A.D. The time of the 


coordinates is taken as midnight on O January of that year. 
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Failure to account for differing epochs can lead to 
substantial errors in the celestial coordinates of a body. 
For example, precession alone will cause an error of 
approximately seven minutes of time in right ascension and 
five minutes of are in declination between 1992 and 2000. 

The Astronomical Almanac (Reference 1), updated every 
year, also lists bright stars, with the coordinates given at 
one-half of the way through the current year. A bright star 
listed in the 1992 Astronomical Almanac would then be 
referenced to epoch 1992.5. 

Of the corrections made to correct a set of coordinates 
for time variance, precession is the one commonly referenced 
to an epoch. The others, nutation and aberration, may be 
determined separately, through formulae which are known for a 
given interval of years. Refraction, obviously, does not 
change with the Earth's position, but is a function of the 
atmosphere and local stellar coordinates. 

Precession corrections in this thesis to the current 
epoch, i.e. the current date and time, are made by referencing 
a star to J2000.0, then back to the current epoch. This 


assures a commonality of treatment for all stars entered. 


PAM 








IV. CORRECTIONS TO CELESTIAL COORDINATES 

The Earth is not a perfect sphere of uniform density (i.e. 
a centrobaric body). This causes close celestial bodies to 
exert disturbing gravitational moments on it. Disturbing 
moments with long periods of oscillation are called 
precession, those with short periods are called nutation. 
These moments cause the Earth's pole to "wobble" in relation 
to the plane of the ecliptic. This motion must be taken into 
account for precise angular determination of stars. 

Movement of the Earth about the Sun causes apparent 
displacement of a celestial body called annual aberration. At 
the degree of accuracy required, aberration effects, too, must 
be used to correct a star's apparent position. 

The Earth's atmosphere changes in pressure and temperature 
regularly, causing changes in the refraction of light incident 
upon it. Refraction must be modeled and its effects taken 
into account. For these and other reasons, a star's tabulated 
position must be modified in order to determine where it will 
actually appear in the night sky. The sight must be corrected 
for the phenomena of refraction, precession, aberration and 


micat ion. 


ZS 


Table i 5) 


| 
Correction Typical Size || Cause Residual 


Error after 


correct ilom 












Precession: Gravitational 
- Tunisolar |f < 50 arc"/yr'lattraction of \<203) arc™7aae 
= planetary <.09 are" /yralicelestial | N/A 
bodies 
Nutation Lunar <0 alae 
Gravitatienal 
| ater ae On 
VWAberration: 
mM = anima. 20.5 are’) hmotion of li< 02) saacen 
- diurnal OW oe renarel” Al Meteo N/A 
Refraction < e201 are™ Earths <T! are! 


atmosphere 


Proper Motion of 


Met Lome L stars 


stars 





Each of these phenomena will be described in detail and 
the necessary corrections explained in this chapter. Of these 


corrections, only precession 1S applied to the virtual 
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telescope. This is because precession is readily calculable 
by an astronomer and 1S commonly calculated for each star 
position taken from a tabulated reference. Nutation, 
aberration and refraction are calculation-intensive and are 
solely the responsibility of the computer. These corrections 
are applied to the physical telescope, but are not shown to 
ine -aStronomer on the status screen. The user is thereby 
insulated from the minutia of the control system and can 
concentrate on an intuitive understanding of the telescope's 
Mesition. 


Table (1) summarizes corrections to celestial coordinates. 


A. PRECESSION 
The Astronomical Almanac, Reference 1, defines precession 


as: 


"the uniformly progressing motion of the pole of a freely 
rotating body undergoing torque from external 
gravitational forces. In the case of the Earth, the 
component of precession caused by the Sun and Moon acting 
on the Earth's equatorial bulge is called lunisolar 
precession; the component caused by the action of the 
planets is called planetary precession. The sum of 
lunisolar and planetary precession 1s called general 
precession." 


Precession of the Earth's axis 1s caused primarily by the 
effects of lunisolar moments due to the much larger 
gravitational attractions of those bodies. One precessional 
period is approximately 26,000 years and accounts for the 


apparent changing of the "North" star in different centuries. 
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Gravitational attraction would be symmetric abotieima 
perfect sphere of uniform mass distribution and no precession 
would therefore occur. Since the Earth 1S not a wWnialtenn 
sphere, attraction to distant bodies varies with density and 
distribution of mass. On the first order, Earth maybe 
pictured as having a belt, or bulge, about the equator. 

Figure 9 shows the 
"misalignment" of the 
Earth's North pole with 
the pole Or the 
EGim perc. The angle e€ 
is called the obliquity 
of the ecliptic and is 
the angle between the 
Barthus equatorial 
bulge and the plane of 
the ecliptic. Figure 


10 shows the effects of 





such an arrangement. 


Figure 9 - Effects of nutation on 
coordinate systems 


a the figure, the 
lower Weigt Or the 
ellipsoid would be more stronger attracted to the point O than 
the upper part due to a smaller distance OP. This causes a 
moment counterclockwise acting on the body G. This moment 


would cause a precession of the pole on body G. 
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Precessional effects, although of very long period, are 
noticeable in astronomical terms. Errors on the order of arc 
minutes may be observed for every few years of difference from 
the referenced epoch if precessional effects are not taken 
into account. Since the MIRA requirement of a few arc seconds 
accuracy is far below that, precession must be taken into 


AGceount. 





Figure 10 - Causes of nutation 


Star data in the Astronomical Almanac isS_ generally 
referenced to an epoch corresponding to the middle of the 
current year. Other star catalogs typically represent star 
coordinates from the standard epochs of B1950.0 or J2000.0 
(See Chapter III for an explanation of epochs). Coordinates 
need, therefore, to be referenced to a Known standard, then to 
the current date for uSe. The approach used referenced all 


Stars first to the J2000.0 epoch and then to the current date. 


Zs 


Two precessional corrections were then made for each set of 


star coordinates. 


From the Astronomical Almanac, page B19: 


FOr, Geaducel1on tov acogo 0: 


C=O MoNCI (aya TCO, 


§ =5-Ncos(a,) Ce 
a =a-— (M+Nsin(a) tan(6)) 
7 (11) 
§ ,=6- = (NCOS (Gy) ) 
For reduerlen fromuZ 00070: 
® .=a+M+Nsin(a,)tan(6,) A 
§ ,=5+Ncos («,) a 
a =a++ (M+Nsin(a) tan(8)) 
“ (13) 


8,=8+ = (Ncos (a,,) ) 


where the subscript zero refers to epoch J2000.0 and the 
subscript m refers to the mean epoch. 


The precessional constants M and N are given by: 


M=1.,,.2812323°7+ 0), 0002.8 494 = Ob Clee Omega 7c (14) 
N=0.5567530°T-0.00011785°T--C4G000116°T° 


where T = (Jullan Date - 245 1545.0)/36 525 
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The precessional coefficients of M and N change less than 
0.00000005 degrees per year around the J2000.0 epoch. Since 
it would take approximately 50 years to accumulate errors on 
the one-tenth arc second order of magnitude, no attempt was 


made to change the precession coefficients annually. 


B. NUTATION 

Nutation, or "nodding", is another name for short-period 
oscillations in the precession moments. For) thee Barth, 
approximately 18.5 years describe one nutation period. 
When compared to the Moon, the Sun has almost no short-period 
effects and they are normally neglected. The effect of the 
Moon's gravitational attraction to the Earth's ellipsoidal 
shape is therefore responsible for the "nodding" of the 


Peer’ S pole. 


A¥Y=-0.0048°sin(318.5°-0.053a) 
-0.0004°sin(198.8°+1.971a) (15) 

NE=40 0026 “Cos (3 13. 5°=0 . 053.4) 
#O.0002"9COS (198 -6°r1..97 1c) 


Julian Date - 244 7891.5 


where ad 
O 


€ = 23.44 = Obliquity of the ecliptic. 


The corresponding rotation matrix 1S 


PS 


fl -A¥cos(e) -A¥Psin(e) 
N=|AY¥cos (e) i ~Aeé (16) 
AY¥sin(e) Aeé a 


The celestial coordinates may be corrected for nutation by 
multiplying the rotation matrix and the coordinate vectomund 


rectangular form: 


S,=NSs, (17) 


The new coordinates must then be decoded from rectangular form 


to right ascension and declination. 


C. ABERRATION 

Stellar objects may be apparently displaced from their 
true, geometric positions by the motions of the observer and 
the object. The finite velocity of light may also cause the 
apparent location of a body to differ from its aemaem 
pes ht ion. The combination of these effects is called 
aberration. Of the motions involved, the only one which is 
Significant to the accuracy required iS annual aberration 
Annual aberration is that component of aberration caused by 
the motion of the Earth vaboue enesoun 

First order annual aberration corrections to right 


ascension and declination are given by: 
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Aa=_KSintA) cos (€) cos (a) 


Goo lo) (19) 
NOe-KSta A) COS (@) Sind) 
Pecos.) |COS (6) Sinta) simi o))—sin(c) cos (6) } 


where A = true longitude of the Sun, the angle from the 
First Point of Aries to the Sun, measured eastward in the 
plane of the ecliptic. 

€ = obliquity of the ecliptic 
and the constant of aberration is given by: 


i 
PC 





(20) 
=20.496 arc seconds 
where a = 1.49598 x 10” cm = 1 AU 
e = 0.01675 = mean eccentricity of the Earth's orbit 
P = 3.1558 x 10’ sec = length of sidereal year 
c = 2.997925 x 10” cm/sec = velocity of light 
The true longitude of the Sun is given by: 
A=278.926°+0.98567d (21) 


where ad = day of year + fraction of day since o” UT 

The true longitude of the Sun is an approximation good to 
about one arc-second. However, since it affects aberration 
only through a product of trigonometric expressions, the above 
approximation of aberration is adequate to the one-tenth arc- 


second level. 


eek 


D. REFRACTION 

When a photon transits from one medium to another in its 
travels, it -réfracts. In other words, the direction Of mime 
travel is alteréd. For the case of a photon traveling throug 
Space from a distant star and impacting the Earth's atmosphere 
(a denser medium), the photon is refracted toward a vector 
normal to the atmosphere. This phenomenon is shown in Figure 


11 below and described by Equation 21. 


Incoming Ray 


Refracted Ray 





Figure 11 -" Retraction ceirough a medium 
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Hp, SIN(®) =p,SIN(O) (22) 


mere fl = The index of refraction of space (= 1), 
u2 = The index of refraction of the atmosphere, 
@ = the angle of incidence between the incident ray and 


the atmospheric normal, 
© = the angle of refraction between the refracted ray 


and the atmospheric normal. 


Since the Earth's atmosphere iS constantly changing, 
refraction is dependent upon many factors and a single index 
of refraction is not possible to determine. Of these factors, 
the most prominent are atmospheric pressure and temperature 
and the wavelength of the observed light. 

The MIRA telescope operates with zenith angles of less 
than 75 degrees for safety of the primary mirror. The 
following formulae allow a reasonably close approximation of 
atmospheric refraction for zenith angles less than 75° (more 


than 15 off the horizon) (See Reference 5). 





Gan (2) 
1 0057 |, 
Zs {2+0. E france.) (24) 
R= 
27 Ste 


eS 


Ksec* (6) Sin(@#) 


«Ean (6) tan(®) +cos (H) ps 


where a = star's uncorrected right ascension 
@, = observer's right ascension 
= GHAy - observer's longitude 
P = atmospheric pressure (mmHg) 
t = atmospheric temperature (°C) 
A = center wavelength of star's light (pz) 
6 = star's uncorrected declination 
He t= estar’s qaeum ganiepiew = +a ee 
@ = observer's latitude 


. K(tan(6) cos (#) ~tan(®)) 
de tan (6) tan(@®) +cos (H) a 


This treatment of refraction can be extremely sensitive to 
changes in the atmospheric parameters. The total refraction 
correction can change on the order of one-tenth of an arc 
second in less than 1.7 seconds under nominal atmospheric 
temperature C10 1c) and pressure (633 mmHg) at a nearly zero 
zenith angle. Refraction was therefore updated every one 
second, with the clock displays. Temperature changes of a 
Single degree Celsius or pressure changes of less than a 
Single millimeter of mercury can alter the refraction 
corrections by approximately one-tenth of an arc second at 
nearly nominal atmospheric conditions at a nearly zero zenith 


angle. While manual entry of these parameters is acceptable 
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for short period observations, eventual upgrade to automatic 


temperature and pressure sensing is recommended. 


5D 








Vv. HARDWARE 


A. SYSREM OVERVIEW 

The hardware to be added to the telescope system is shown 
in Figure 12 and may be compared to Figure 1, the current 
hardware configuration. 

An IBM-compatible personal computer, based on an Intel 


80486 microprocessor, 1S to be used as the controlling 





a i | Encoders 
Se oa —— 
Telescope + “=<. == 
ee Limit . 
Spcciei: 
2 Stepping | 2 slew , Slew 
Motors | ___ Motors Paddle #1 
Motor Motor - 
Drivers Controller | 
, . | ~Set/Guide 
Paddle #1 
sn 
. | 
_ Computer . 
Keyboard a Lvicalion: 
WWV 
Synch 
Figure 12 - Proposed hardware configuration 


authority. System clock time is to be updated automatically 


via the WWV time synchronization board discussed in Chapter 
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jEEr Slew and set/guide paddles are to be interfaced to the 
computer via a TTL interface board. Motor drivers are either 
in an "on" or "off" state. The state is to be controlled by 
relays operated through an additional TTL interface board. 
Physical limit switches are installed to prevent the 
telescope from being moved beyond its limits in the case of 
software failure. Two sets of physical limit switches are 
installed. The first set of switches cuts power to the slew 
motor controller in the direction required when the 
telescope's zenith angle reaches ow Operation in the 
opposite direction is allowed to facilitate recovery. This 
set is engaged first and will stop telescope slew. The second 
set of limit switches cuts power to the stepping motor 
drivers, thus preventing the stepping motors from driving the 
telescope to an unsafe position. In this case, also, only 


movement in the affected direction is stopped. 


B. COMPUTER 

An IBM-compatible personal computer was chosen for the 
system due to availability, cost, speed, and operating 
considerations. Of these, familiarity of the MIRA staff with 
the IBM-compatible operating system and the Microsoft 
QuickBASIC programming language were the key factors. A 
machine based on an Intel 80486 microprocessor was chosen to 


allow the system to be operated in a multi-tasking 
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environment. This eliminated the need for a dedicated 
computer, cost again being a major consideration. 

Ses allow multi-tasking, the POlLows) ne hardware 
requirements are recommended as a minimum: 

- 80486 microprocessor 

* 4 MBytes RAM 

mel hard disk drive 

- 1 floppy disk drive 

wweGck monitor 

* a mathematics co-processor if a lower CPU is used 

If multitasking is not desired, the computer requirements 
may be downgraded to 1 MByte RAM and an =<Intel 80286 
microprocessor with an 80287 mathematics co-processor. 

This type of computer comes standard with eight bus 
interface slots. Of these, at least five would be needed. 
One slot would be required for each of the TTL boards for 
interface to the two control paddles and the two types of 
motor controllers. One additional slot would be required for 
the WWV board. Additionally, one RS-232 standard serial 
interface port would be required for the encoder inputs. This 
type of port also comes standard with almost any brand of IBM- 


compatible computer. 
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VI. PROGRAM DEVELOPMENT 

A computer program was developed in the Microsoft 
QuickBasic Language, version 4.00, to point the telescope to 
predetermined celestial coordinates and correct for any errors 
larger than one-tenth of an arc-second. The program uses 
rotation matrices, developed in Chapter II, to translate 
between celestial and terrestrial coordinate systems. 

The program structure 1S modular. The main program 
controls 24 subroutines which deal with data display, data 
entry, data correction, time correction, key and video 
handling, and telescope movement. Each subroutine is 
described below, within its appropriate category. 

The concept of a virtual telescope was fundamental to the 
software development. A '‘'virtual' model of the systen, 
encompassed solely in software, was developed and the effects 
of any operation on this model were evaluated prior to any 
action being taken. This 'virtual' model is seen by the user 
and represents the target star's position corrected solely for 
precession. Corrections to raw star position data for 
precession, Nutcatson, aberration, and refraction were 
completed and new motor states determined. These corrections 
and motor states determine the status of the physical 
telescope. The new position of the physical telescope was 


determined in both terrestrial and celestial coordinates. If 


41 


no errors were encountered, such as hard- or software limits 
to movement, the telescope was then issued motor commands. 
Fig 13 shows the highest-level modular flow of the 


program. 


A. MAIN PROGRAM DESCRIPTION 

The main program, beginning on page 1 of Appendix A, 
controls the flow of telescope operation and includes the main 
operating loop. A simplified flow chart of the main program 
is shown in Figure 14. 

The program begins by initializing system hardware. It 
then proceeds to allow user entry of daily atmospheric and 
observation parameters, if required. The telescope is aligned 
by pointing to a known bright star. The star's celestial 
coordinates nee entered, thereby giving the system a reference 
point. Data entry of desired stars for view is accomplished 
either by direct entry or by reading a previously written and 
formatted file. Editing of a stored file is allowed. A star 
is then selected for view. Once a star is selected, the 
program ensures that the star is visible above the current 
local horizon. If it is, the main loop is entered, allowing 
the system to find and track the star. If it is not, another 
star must be selected. 

Since safety of the telescope structure was the primary 
design criterion, several layers of protection were utilized 


in the hardware and software. At any time that a zenith angle 
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Figure 13 - Program modular concept 
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Figure 14 - Simplified flow chart of the main program 
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of 65 degrees or greater was calculated, a warning message 
would appear to the operator on the telescope status display 
screen. When a 75 degree zenith angle was reached, the 
software would automatically drop track on the target star and 
ask the operator to choose a new star for view. In the event 
of catastrophic software errors that allowed the telescope to 
move beyond a 75 degree zenith angle, physical hardware limit 
Switches were installed at a 76 degree zenith angle. Two sets 
of limit switches exist on the telescope. The first set cuts 
power to the slew motor controller. This set of switches must 
be engaged first due to the speed involved in slewing and the 
eae that the system requires to decelerate. The second set 
of switches cuts power to the stepping motor drivers. 

A simplified flow chart of the main program loop is 
provided in Figure 15. In this loop the main tracking actions 
take place. Every second, time and position information is 
updated to the telescope status screen. Pie wearget estar s 
right ascension, declination, hour angle, zenith angle and 
azimuth are displayed. Movement commands are issued to keep 
the differences between the stars' coordinates and the 
telescope's within acceptable limits. The loop also checks 
Once per loop for user input. An ESC key hit stops telescope 
movement immediately and awaits user confirmation to quit the 
program. If confirmation is received, the program is exited 
via a procedure to safely shut down all hardware. If the 


ENTER key is hit, the change-telescope-status menu appears, 


45 





_ Start D 
“——1---—_Every Second Update Time 
Displays 
Display RA, DEC, HA 
a | ames Oitnalie oes 
| | Update Refraction 
| Calculate and Display ! Corrections | 
Zenith and Azimuth = oa 
ee RTN 
“Error? ~~ Yes_| _ Error — 
= ' Window 
es nl 
No On Ca: 
re epropmate 
Limit Switch Checks | Aca 
ad ee | 


} 





awa ———____—— : 


| Slew Motor Commands. 
—_— ——————— 





—_—— - 


| Step Motor Commands | 
—— —{—— 


< a sMonitogng > 





Figure 15 - Simplified flow chart of the main loop 
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allowing the user several options including the ability to 
proceed to another star. 

Once the virtual telescope's right ascension = and 
declination, the actual telescope's azimuth and zenith angles, 
the star's hour angle and the times are ready for display in 
the main loop, the TeleStatusDisplay subroutine is called to 
display them. It updates the Pacific Standard Time, sidereal 
time, and universal coordinated time clocks once per second, 


along with the updated refraction correction. 


B. SUBROUTINE DESCRIPTIONS 
1. Display Subroutines 
Display subroutines are responsible for the display of 
data and graphics to the computer screen. 
a. InitialDisplay 
This subroutine creates the display screen graphics 
which will be utilized for the remainder of the program run. 
A user box is displayed at the bottom of the screen, in which 
all user interaction and display occurs. The user box, as 
shown in Figure 16, is simply a screen area for data exchange. 
Additionally, menus and error messages occur in windows which 
occupy other video screens. 
This subroutine also colors the background and identifies the 


program and its intended location to the user. 
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Figure 16 - The initial display screen with user box 
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b. StarDisplay 
This subroutine modifies the initial display screen 
for star data entry. It draws the individual fields in which 
the parser (see subroutine StarDataEntry) will operate and 
also prints instructions for leaving the data entry routine. 
c. Windows 
All error messages and most menus are displayed on 
video pages which are separate from the primary display page. 
These messages and menus are pre-defined data structures and 
are displayed by this subroutine. A rectangular window is 
placed at a position determined in the data structure. The 
size of this window is dependent on the number of lines of 
text to be displayed and the length of the longest line. 
To determine the video parameters and if the next 
video page is available, the subroutine VideoState is called. 
A window 1S removed by calling the WindowsPop 
subroutine. 
d. WindowsPop 
This subroutine removes the last window placed by 
the subroutine Windows from the active video page and makes 
the next lowest video page active. 
e. TeleStatusDisplay 
This subroutine modifies the user box for the 


display of telescope and star positions and time during the 


49 


main loop. The name of the active star is placed in the 
appropriate box. 
f. ShowTime 

There are three different times displayed on the 
status screen during telescope movement, universal coordinated 
time (UTC), mean sidereal time and Pacific Standard time 
(PST). These three times are updated and displayed through 
this subroutine. 

g- TeleAngles 

Using the celestial coordinates to which the 
telescope is pointed, this subroutine determines the position 
of the telescope in terrestrial terms. The terrestrial 
rectangular coordinates are decoded into azimuth and zenith 
angles. Azimuth is defined as the number of degrees measured 
eastward in a circle from North. North is therefore 0 
degrees, East is 90 degrees, etc. Zenith angle is defined as 
the number of degrees from the point directly overhead the 
telescope to the telescope pointing vector. Zenith may also 
be thought of as the co-elevation, where zenith equals 90 
degrees minus the elevation. Elevation is the number of 
degrees from the local horizon to the telescope pointing 
vector. 

The above angles are calculated by means of 
direction cosine matrices. These matrices are described in 


detail in Chapter II. 
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2. Data Entry Subroutines 
Data entry subroutines all accept user input prior to 
proceeding. User input is necessary to align the telescope, 
to provide the program with atmospheric parameters anda list 
of stars on which the telescope shall sight. 
a. TeleAlign 

The telescope must have an alignment point to which 
all future movements will be referenced. This subroutine 
accepts user input to determine the right ascension and 
declination of the telescope pointing vector at a given time. 
Coordinates from encoder output Signals are then treated as a 
difference from this baseline. 

Typically, the right ascension and declination of 
the present pointing vector will be entered from the bright 
star list of the Astronomical Almanac. Therefore, the epoch 
of the coordinates must be given and accounted for. This 
Subroutine precesses, nutates, aberrates and refracts the 
given coordinates to account for the difference between 
tabulated and apparent positions of the bright star chosen for 
alignment. 

b. Atmospheric 

Refraction through the atmosphere depends primarily 
on four variables: the zenith angle, the atmospheric 
temperature, the atmospheric pressure and the center 


wavelength of light from the star being observed. The latter 
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three must be provided to the program in order to be accounted 
for. This subroutine allows user entry of these variables. 
c. StarDataEntry 

This subroutine is a parser for accepting data for 
up to one hundred stars. Four fields are required for each 
listing: star name, right ascension, declination and the 
epoch to which the right ascension and declination refers. 
The subroutine displays five lines of the listing at one time 
and allows scrolling, page jumps and home or end-of-file 
jumps. Full word-processing ability has been included for 
appropriate keys. Insert, delete and shift keys are active, 
with their usual functions. 

No attempt was made in this procedure to ensure the 
correctness of the formatted data. If the data formats are 
not correct, the program will later display an error message 
and allow the listing to be edited. 

As the listing is updated, the matrix variables 
Star$(), RAS$(), DECS() and EPOCHS() are updated. These four 
variables are each one hundred elements in length and hold the 
star names, right ascensions, declinations and epochs, 
respectively, in string form. 

Figure 17 shows the StarDataEntry screen. 
d. ChooseStar 
Once a listing of stars has been read into the 


matrix variables Star$(), RAS$(), DECS$() and EPOCHS, one star 
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must be chosen for view. This is referred to as the "active 
star". An active star will therefore have associated with it 
a name, a right ascension, a declination and an epoch. 

This subroutine modifies the user box to show a 
listing of the one hundred star entries in the current file. 
This file will be on either a floppy or hard disk and may have 
any file name. The default name is "STARFILE.DAT". The 
subroutine allows the user to scroll, page up or down, jump 
home or end-of-file. Once the selected star is highlighted, 
it is chosen for view by pressing the F10 key. The active 
star's four elements are then passed to the main loop. 

3. Data Correction Subroutines 

As described in Chapter IV, the movement of Earth 
around the Sun, the wobble of Earth's axis, the rotation of 
the Earth, the relativistic effects of the Sun's presence and 
the refraction of light through the atmosphere must be 
accounted for. The data correction subroutines make 
adjustments to the right ascension and declination of the 
active star to account for these effects. 

a. Precession 

The effects of precession on astronomical sightings 
were discussed in Chapter IV. The Precession subroutine uses 
the approximate formulae for precession found in the 
Astronomical Almanac, page B19. The epoch of the active star 


is first precessed to the standard epoch of J2000.0, then to 
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the current epoch. The current epoch is the calculated to the 
current day. Right ascension and declination for the current 
star are then modified to account for the difference between 
the epoch entered by the user and the current’ epoch. 
Precession is accounted for at the beginning of each star 
Saiciite. 
b. Nutation 

Nutation is described in Chapter IV. The Nutation 
subroutine uses the approximate formulae for nutation found in 
the Astronomical Almanac, page B20. A rotation matrix is used 
to relate the unmodified and nutated celestial rectangular 
coordinates. The rectangular coordinates are then decoded 
into right ascension and declination. Nutation is accounted 
for at the beginning of each star sight. 

c. Aberration 

The causes of aberration are briefly discussed in 
Chapter IV. This subroutine accounts for aberration by 
determining the relation between the apparent positions of the 
Sun and the active star and modifying the right ascension and 
declination accordingly. The method used was garnered from 
Reference 11, pages 499-502. The mean longitude of the Sun 
was calculated from formulae provided in the Astronomical 


Almanacypager eZ . 


one 


d. Refraction 

The effects of refraction through the Earth's 
atmosphere were considered in Chapter IV. Refraction is 
affected by four variables: zenith angle, atmospheric 
temperature, atmospheric pressure and the center wavelength of 
light for the observation. The latter three were entered by 
the user in the Atmospheric subroutine. The first was 
calculated in the TeleAngles subroutine. 

Refraction was calculated by the method employed by 
Reference 5, pages 86-87. Refraction corrections to right 
ascension and declination were calculated in the subroutine 
and applied in the main loop. This is necessary because 
refraction 1s updated once per second and old corrections must 
be removed prior to new ones being applied. 

4. Time Correction Subroutine 
a. SiderealTimeCalc 

Sidereal time is a meaSure of time in relation to 
the stars (specifically to the vernal equinox), instead of in 
relation to the Sun. Local mean sidereal time is the actual 
number of hours, minutes and seconds from the alignment of the 
local meridian to the mean equinox. This subroutine 
calculates the local mean sidereal time and adds a minor 
correction to account for the difference between the true and 
mean equinoxes. This yields the local apparent sidereal time. 


This time is then displayed from the main loop. 
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5. Key and Video Handling Subroutines 
a. KeyCode 
For any key pressed, KeyCode returns a unique 
integer value. This value may then be used to determine the 
key or key combination that has been pressed. Once the 
subroutine is entered, the program will wait until the key 
buffer contains a character before returning control to the 
calling routine. This is especially useful in any application 
requiring the decoding of user response. The following call 
KeyCode or KeyCodeNowWait: The main program, ChooseStar, 
StarDataEntry and Windows. 
b. KeyCodeNowait 
KeyCodeNoWait is identical to the  KeyCode 
subroutine above except that the key buffer is checked and 
control is immediately returned to the calling routine instead 
of waiting for key input. The main loop uses KeyCodeNoWait to 
monitor the keyboard while the telescope is tracking a star. 
c. VideoState 
The video mode parameters of the computer on which 
the program is being run must be determined in order to place 
the error message windows and menus. VideoState determines 
the current video mode, the number of text columns and the 
number of the active display page. The Windows subroutine 
uses this data to determine if there is room to place a 


message window prior to placement. 
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6. Telescope Movement 

Four subroutines control motor commands to the 
telescope. Each subroutine issues motor commands for a 
certain type of movement. Movement is categorized into slew, 
set, track and guide. Slew uses the one-quarter horsepower 
slew motors to transit the telescope through large angles. 
This places the telescope in the neighborhood of a target 
star's position. Set, track and guide all use the smaller | 
stepping motors. Set is used to adjust the telescope's 
position after slewing and align the telescope's pointing 
vector to the target star. Once slew and set have occurred, 
track changes the right ascension of the telescope at a 
Semstant rate to account for the Earth's rotation. As 
refraction, movement and other errors cause misalignments, 
guide adjusts the telescope's position. 

To ensure structure safety, the position of the 
gravity limit switch is checked after every step of the 
stepping motors and periodically during use of the slew 
motors. 

Arrow keys are monitored from the main program before 
each track update to allow manual adjustment of the 
telescope's position. 

a. SlewCommands 

The SlewCommands subroutine, immediately after a 


target star is chosen, slews the telescope to within one 


So 


degree of the calculated position. On successive loops, 
movement is left to the stepping motors in the set (fast) or 
guide (slow) modes. 

The maximum slew motor rate is approximately 12 
degrees per second. To ensure the safety of the telescope, 
the gravity limit switch position is checked every degree of 
movement or every eight-hundredths of a second during slewing. 

b. SetCommands 

This subroutine commands the stepping motors 
following a slew to a new target star. Once the SlewCommands 
routine has positioned the telescope to within one degree of 
a calculated position, this subroutine uses the stepping 
motors, in the fast mode, to bring the position to within one 
step of the star's calculated position. The astronomer 
operating the system will then guide the telescope to the 
star's actual apparent position by operating the stepping 
motors manually, either by the hand paddles or the arrow keys. 

c. TrackCommands 

Once a star has been found, the TrackCommands 
subroutine will update the right ascension motor stepping rate 
to adjust for changes in the atmospheric refraction. The 
telescope will, at this time, be slowly sweeping in right 


ascension to compensate for the rotation of the Earth. 


60 


C. ERROR AND WARNING MESSAGES 

Error handling iS an important aspect of any computer 
program. To prevent a program from halting catastrophically 
upon encountering an error, errors must be anticipated and 
dealt with appropriately. Additionally, some situations may 
be anticipated that require that the user be warned of 
impending problems. The main program contains ten error and 
warning windows to advise the user of anticipated situations. 

1. The Quit Confirmation Window 

Following standard practice, the ESC key is used to 

abort a telescope tracking run or to exit the program 
entirely. Upon striking the ESC key telescope movement is 
halted, if the telescope is moving. The quit-confirmation 
window then appears. If the user types either an upper or 
lower case "Y", the program is terminated. Otherwise, control 
returns to that segment of the program previously running. 

2. The File Save Error Window 

When a data file is to be saved, many errors may 

occur. The appropriate disk drive may not be functioning, the 
drive may not contain a disk, the disk may be write-protected, 
the disk may be full, etc. In any of these cases, the file- 
Save-error window appears to inform the user that the data 
file could not be saved. Upon the press of any key, control 
returns to the file editor. From this point the user may 


correct the error or quit the program, as appropriate. 
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3. The File-Not-Found Error Window 
When the user requests that an old data file be used, 
the program requests the path and file name top retrieve. If 
the path does not exist or if the file name cannot be located, 
the file-not-found-error window appears. This window informs 
the user of the error and waits for any key to be pressed. 
Control is then returned to the data entry screen menu, 
allowing the user to choose an old file, re-enter the file 
name and path, or begin a new file. 
4. The Equinox Correction Error Window 
Corrections to the Greenwich mean sidereal time are in 
the form of a number, in seconds. The default values provided 
are zero seconds for both the current and next days. Should 
the user fail to except the default values and enter a 
correction that could not be interpreted aS a number in 
seconds, the equinox-correction-error window appears. once a 
key 1S pressed, the user is allowed to except the defaults or 
re-enter new numbers. 
5. The Slew Confirmation Window 
Before the telescope is slewed to a new star, the new 
target star's coordinates are calculated. The zenith angle 
anticipated is displayed, along with a warning that the Fl key 
will cause the telescope to slew. Although not displayed on 


the screen, the F10 key will cause the program to run in 
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SIMULATION mode. Any other key will return the user to the 
ChooseStar subroutine. 
6. The Change Telescope Status Window 
The change-telescope-status window is a user menu of 

seven possible actions. Tt appears when the ENTER key is 
pressed during a telescope tracking run and allows the user to 
take one of the following actions: 

- Edit the current list 

- Choose a star from the current list 

* Change to a new data file 

‘« Change atmospheric parameters 

* Park the telescope 


- Exit the menu and return to whatever tracking was 
previously taking place. 


7. The RA/DEC Format Error Window 
Right ascension and declination of a target star are 
entered from a formatted data file. If the program attempts 
to decode these coordinates and fails, the RA/DEC-format-error 
window appears. The user is then given the opportunity to 
edit the data file or choose another star. 
8. The Slew Error Window 
If the user ignores the slew-confirmation window with 
a zenith angle over 75 degrees and attempts to slew the 
telescope to an unsafe position the slew-error window appears. 


No telescope movement is allowed. Upon the press of any key 
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control is returned to the change-telescope-status menu so 
that the user may choose another star or take other 
appropriate action. 
9. The Zenith Angle Warning Window 
The maximum zenith angle allowed is 74 degrees. When 
the zenith angle reaches 75 degrees telescope movement is 
halted. To warn the user of impending cessation of movement, 
the zenith-angle-warning window appears any time the zenith 
angle is 65 degrees or over. It displays the simple warning 
"The zenith angle is X degrees". 
10. The Zenith Angle Halt Movement Window 
Should the user ignore the zenith-angle-warning window 


and allow the telescope to track to its limits, the zenith- 


angle-halt-movement window appears at 75 degrees. All 
telescope movement is halted. The change-telescope-status 
menu appears upon the press of any key. The user may then 


take action as appropriate such as choosing another star. No 
action that the user takes can cause the telescope to 
automatically track to a zenith angle greater than 74 degrees. 
11. The Telescope Status Help Window 

When the telescope status screen is shown, help for 
the allowed key strokes is available through this subroutine. 
The Fl key will display a window listing the keys which 
control manual slew or guide of the telescope. The arrow keys 


and numeric keypad is used for this purpose. 
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12. The CTS-10 Failure Window 

diate is possible that the CTS-10 WWV time 
synchronization board will fail and reset. This usually 
indicates an onboard battery failure or a poor physical 
connection to the data bus of the computer. In this event, 
the CTS-10 will reset the computer's system clock to the 
factory default setting. This condition is tested in the main 
program upon initialization. If an error is detected, an 
error window will appear, advising the user to either repair 
the CTS-10 problem or manually update the system clock to 
correspond to UTC. After any key is pressed, the program will 
halt execution and must be run again after the error is 


mMeetitLled. 
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Vil. PROCEDURES 

Once final hard- and software modifications have been 
completed, an astronomer shall have automatic control of 
telescope movements to a known set of celestial coordinates. 
He or she will remove the telescope from its stowed 
configuration by removing obstructing objects, sliding the 
enclosing roof away from the telescope and opening the doors 
covering the primary mirror. After this point, procedures 
differ drastically from the current procedures. 

The astronomer will provide power to motor controllers and 
the computer. Following systems checks, the MIRACTRL program 
will be run. The user will provide the celestial coordinates 
to which the telescope is aligned by looking at a known bright 
star. The astronomer will enter atmospheric data for the 
evening and provide the computer with a list of star positions 
either by a data file or by manual entry. He or she will 
choose a star to view and the program will check whether it is 
in view. If so, a Single key stroke will slew and set the 
telescope to its desired position. The astronomer will guide 
the telescope's position onto the star and allow the computer 
to begin tracking. The computer will continue to track until 
the astronomer guides the telescope manually, the star is out 
of view or the gravity limit switch is engaged. The 


astronomer may choose to halt the telescope, exit the program, 


C7 


thereby halting the telescope, or change the target star at 
any time. 


The telescope status display screen is shown in Figure 19. 
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VIII. PROGRAM VALIDATION 
To the date of this writing, the computer program written 
had not been tested on the telescope due to lack of funds and 
the time reguired to connect the associated hardware. 
Extensive testing of each subroutine that affected telescope 


movement was accomplished in simulation. 


A. TRACKING VALIDATION 

To ensure accuracy of the program's tracking commands and 
related coordinate systems, several simplified simulations 
were performed. Initially, the latitude of the observing site 
was set to zero degrees and a star entered with a zero 
declination. Corrections for precession, nutation, refraction 
and aberration were not included. This ensured that, with 
properly defined coordinate systems, the star would have a 
zenith angle pointing directly East or West at any time. The 
Simulation was run and meridian passage of a simulated star 
was observed. It was verified that at meridian passage, the 
local hour angle became zero and switched from a negative to 
a positive sign. It was also observed that the zenith angle 
shifted from East to West. 

In the next validation simulation, the observing site's 
latitude was changed to the correct latitude at the Oliver 


Observing Station. The above simulation was repeated and 
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meridian passage of the simulated body observed. A target 
star from the Bright Star listing of the Astronomical Almanac, 
Reference 1, was then entered and its calculated azimuth and 
elevation verified by use of a U.S. Navy Rude Starfinder. 
This process was repeated in each quadrant of right ascension 
and for both positive and negative declination while varying 
the system clock to allow for each target star to be in view. 

When the system clock was restored to its correct time and 
the target star was changed to different quadrants, the 
appropriate error messages were observed when the star was not 


in ev lew. 


B. CELESTIAL COORDINATE CORRECTION VALIDATIONS 

The correction of celestial coordinates was verified by 
allowing only a single correction to be applied at a time. 
The correction under scrutiny was observed to be within the 
expected limits set in Table 1. Precession, nutation, 
aberration and refraction corrections for a given target star 
were calculated by hand and the results compared to the 
computer solution. The direction of each correction was 
verified at concurrently. For example, refraction of light by 
the Earth's atmosphere always causes the apparent position of 
a star to be shifted toward the zenith. This was verified in 
all simulations run. 

Corrections to celestial coordinates must take place in 


the proper order, i.e. precession, nutation, aberration and 


v2 


Bettaction. Simulations were also run to verify that the 
corrections were taking place in the proper sequence and that 
the only changes to the celestial coordinates came from those 


sources. 


C. MOTOR COMMAND VALIDATION 

Motor commands could not be fully validated since the 
computer waS not connected to the motor drivers or 
controllers. Simulations to verify motor command logic were 
completed. The telescope status display screen also displays 
the motor states for all four motors. These displays were 
used to determine when motors were turned on or off during 
representative actions. Slew, guide and track motor commands 
were tested individually and together to ensure that the 
sequence of commanded motor actions was correct. Following 
hardware installation, a check should be made of signal output 
to each motor controller or driver prior to allowing telescope 
movement. Limit switch actions were also validated by 


Simibation. 
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IX. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

It has been shown that an automatic telescope control 
system is both practical and relatively inexpensive. 
Additional hardware requirements for the system discussed 
would cost approximately eight thousand dollars as of the time 
of this writing. While a personal computer based on an Intel 
80486 microprocessor was used for this design, it was 
determined experimentally that the program will run in real 
time on an Intel 80286-based system, as long as an Intel 80287 
mathematics co-processor is installed. This variation would 
require a dedicated computer, because multitasking on an 
80286-based computer is not practical with available operating 


systems or data bus speeds. 


B. RECOMMENDATIONS FOR ADDITIONAL WORK 

There remains additional work on the MIRA telescope system 
before fully automatic operation is possible. Hardware 
updates yet required include modification of installed optical 
encoders, replacement of the slew motor drivers and computer 
control of the focusing step motor. 

Completion of computer program subroutines to reflect the 
types of motor drivers eventually purchased will be necessary. 


Additionally, a warm-room computer monitor and keyboard as 
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well as slew and guide/track paddles will need to be 


installed. This work was not completed because of lack of 
funds. Figure 20 shows the desired hardware configuration 
upon completion. The pointing model, which represents the 


changing telescope structure in software, must also be 
included. Testing of the completed system will require 
extensive simulation in both hard- and software prior to 
actual telescope movement. Before the telescope is first 
moved by computer control, a rigorous testing plan must be 
developed. 

Automatic environment sensing of temperature and pressure 
is required for long duration tracking of a single body 
because of refractive effects. Less than one-tenth of an arc 
second resolution may be obtained under the current model due 
to manual input of time-varying parameters. 

Optimization of slewing movements may be accomplished by 
a "traveling salesman" approach to minimize the required slew 


times or accelerations. 


var 





APPENDIX A MIRACTRL. EXE 


A. MIRACTRL.BAS LISTING 


PKKKKKKEKKKEKKKEKEKEKEKKKEKKKEKEKKKKKKKKKKKKKKKEER 


tk Name: MIRACTRL a 
"kX Type: Program xx 
tke Module: MIRACTRL. BAS wx 
Vk Language: Quickbasic 4.50 ** 
'kx* Pie nhor : David P. Wood & 


VK KKKKKKKEKEKEKKEKKEKEKKKKKKKEKKKKEKKKKKKKKKKKEKSE 
‘ 


‘Control Core for MIRA Telescope. 


'Usage: No Command Line Parameters 

',MAK File: (none) 

‘Parameters: (none) 

'Variables: w*Texts() Text for the error windows 

: Stars () A 100-element matrix of desired 
stars for 

viewing 

RAS () A 100-element matrix of right 
ascensions 

for the stars listed in Star$() 
DECS() A 100-element matrix of 
declinations 

for the stars listed in Star$() 
Epochs() A 100-element matrix of epochs for 
the 

RA's and DEC's of the stars listed 
ipotar > ( ) 

1 


EQofEQ! () The Equation of Equinox 
Corrections to 
' sidereal time for successive dates 


(due 

: to Gorrections for Ehe Earth's 
nutation 

: EQofEQUIN! The Equation of Equinox correction 
for the 

’ current time 

w.* Windows Type variables, defines 
window 

appearance 

: ming The current minute of time 
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hours 


The current hour of time 
The current day 

The current month 

The current year 


LocalTime$ A string representing the local 


SiderealTime$ A string representing the local 


| 

day% 

' month% 

years 

q 
time 

i] 

q 

: TIMES 
to correspond 

| 

(URE) 

: Kees 
on the keyboard 

t 
KeyCodeNoWait% 

q 

' Characters 
returned 

q 

! Abort. 
keyboard buffer 

q 
user wishes to 

q 

' Ounites 


keyboard buffer 
q 


Abort% is TRUE 


'Define constants 

CONST ENTER = 13 

CONST TABHIT = 9 

CONST BACKSPACE = 8 
CONST ESCAPE = 27 

CONST SHIFTTAB = 3840 
CONST RIGHTARROW = 19712 
CONST LEFTARROW = 19200 
CONST UPARROW = 18432 
CONST DOWNARROW = 20480 
CONST NUMRIGHTARROW = 54 
CONST NUMLEFTARROW = 52 
CONST NUMUPARROW = 56 
CONST NUMDOWNARROW = 50 
CONST one = 49 

CONST three = 51 

CONST seven = 55 

CONST ehnine = 57 

CONST F1 = 15104 

CONST F2 = 152860 


apparent sidereal time 
The system clock time string, set 


to Universal Coordinated time 
A unique code for any key struck 
returned by either KeyCode% or 


Functions 
A string representing a character 


from the keyboard 
A variable assigned to a current 


input used to determine if the 


quit the program 
A variable assigned to a current 


input used to determine if the 
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CONST F3 = 15616 

CONST F4 = 15872 

CONST F9 = 17152 

CONST F10O = 17408 

GONSYT Fill = -31488 

Senst Fi2 = -31232 

CONST DELETE = 21248 

CONST FALSE = 0O 

CONST TRUE = NOT FALSE 

CONST PGDOWN = 20736 

CONST PGUP = 18688 

CONST INSERT = 20992 

CONST HOME = 18176 

CONST ENDFILE = 20224 

CONST PLUS = 1 

CONST MINUS = -1 

CONST ‘ul = 3.14159265359# 

Define color constants 

CONST BLACK = O 

Seis! BLUE = 1 

CONST GREEN = 2 

CONST CYAN = 3 

CONST RED = 4 

CONST MAGENTA = 5 

CONST BROWN = 6 

CONST WHITE = 7 

CONST BRIGHT = 8 

CONST BLINK = 16 

CONST YELLOW = BROWN + BRIGHT 

GyYEeE Regrype 
ax AS INTEGER 
bx AS INTEGER 
Cx AS INTEGER 
ax AS INTEGER 
Bp AS INTEGER 
si AS INTEGER 
di AS INTEGER 
flags AS INTEGER 

END TYPE 

TYPE RegTypex 
ax AS INTEGER 
bx AS INTEGER 
ex AS INTEGER 
ax AS INTEGER 
Bp AS INTEGER 
si AS INTEGER 
di AS INTEGER 
flags AS INTEGER 
ds AS INTEGER 
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es AS INTEGER 


END) LY PE 


TYPE WindowsType 


action AS INTEGER 
edgeline AS INTEGER 
row AS INTEGER 
col AS INTEGER 
fgdEdge AS INTEGER 
bgdEdge AS INTEGER 
fgdBody AS INTEGER 
bgdBody AS INTEGER 
fgdHighlight AS INTEGER 
bgdHighlight AS INTEGER 
fgdTitle AS INTEGER 
bgdTitle AS INTEGER 
fgdPrompt AS INTEGER 
bgdPrompt AS INTEGER 
returnCode AS INTEGER 
END’ Ty PE 
t 
'Functions 
DECLARE FUNCTION KeyCodeNoWait% (Character$) 
DECLARE FUNCTION KeyCode% (Character$) 
t 
'Subprograms 
DECLARE SUB windows (w AS WindowsType, wText$(), wTitles, 
wPrompts) 
DECLARE SUB WindowsPop () 
DECLARE SUB VideoState (mode%, columns%, page%) 
DECLARE SUB InitialDisplay () 
DECLARE SUB StarDisplay () 


DECLARE SUB StarDataEntry (Star$(), RAS(), DECS(), EPOG a 


NewFile%, wl AS WindowsType, wlText$(), wlTitle$, wlPrompt$) 

DECLARE SUB ChooseStar (Star$(), RAS({), DECS(), EPOGHS Ge 
NewFile%, wl AS WindowsType, wlText$(), wlTitle$, wliPrompts, 
currentStar?) 

DECLARE SUB TeleStatusDisplay (ActiveStar$, Simulation’) 

DECLARE SUB DayLightSaving (DayLightSavings%, hour%, day%, 
month%, year?) 

DECLARE SUB SiderealTimeCalc (SiderealTime$, min%, hour%, 
day%, month%, year%, EQofEQUIN!, GMSTOhrSEC#, JulianDate#, 
dayFraction#, DaysPassed%) 

DECLARE SUB EquationEquinoxes (EQofEQ(), EQofEQUIN! ) 

DECLARE SUB Interrupt (intnum%, inreg AS RegType, outreg AS 
RegType) 

DECLARE SUB ShowTime (min%, hour%, day%, month%, year%, 
DayLightSavings%, SiderealTime$) 

DECLARE SUB StarData (StarS() , @RAS@), BPECS(), SSEoeo ee 
EditDecision%, NewFile%) 
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DECLARE SUB Refraction (Pressure!, Temperature!, Wavelength!, 


GMSTOhrSEC#, HourAngle#, GHAaries#, RAofStar#, DECofStar#, 
RARefractCor#, DecRefractCor#, RAofMIRA#) 

DECLARE SUB Atmospheric (Pressure!, Temperature!, 
Wavelength!) 

DECLARE SUB TeleAngles (Zenith%, Azimuth%, GHAaries#, 
RAofStar#, DECofStar#, CelCoord#(), DCM#(), TerraRec#() ) 

DECLARE SUB Precession (RAS(), DECS i EPOCHS(), 
currentStar%, JulianDate#, EpochError%, RAofStar#, DECofStar#, 
RADecError?) 

DECLARE SUB Aberration (RAofStar#, DECOPStare, 
RAAberrationcCor#, DecAberrationcCor#, day baaec ont 
DaysPassed?) 


DECLARE SUB SlewCommands 


RAofStar#, 
DECLARE SUB TrackCommands () 


DECLARE SUB SetCommands 


RAofStar#, 


DECLARE SUB TeleAlign 


TeleDEC# ) 
DECLARE SUB Nutation 
Moe ol#(), S2#()) 


DE COnSieana) 


DECOESEar?# ) 


‘Data Structures 


DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 


‘Arrays 


DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 


w9 


WindowsType 
WindowsType 
WindowsType 
WindowsType 
WindowsType 
WindowsType 
WindowsType 
WindowsType 
WindowsType 


(RAofStar#, 


wl0O AS WindowsType 
wll AS WindowsType 
wl2 AS WindowsType 


wabextS(1 TO 1) 
w2Text$(1 TO 1) 
w3Text$(1 TO 1) 
marexts (1 TO 1) 
mwolexts(] TO 7) 
w6Text$(1 TO 9) 
warexts(1 TO 1) 
w8Text$(1 TO 2) 
w9Text$(1 TO 3) 
wl0Text$(1 TO 2) 
wiibextS(i TO 13) 
w1l2Text$(1 TO 5) 
Star$(1 TO 100) 
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(Simulation, 


(Simulation?, 


(GMSTOhrSEC, 


DECofStar#, 


TeleRA#, TeleDEC#, 


TeleRA#, TeleDEC#, 


GHAaries#, TeleRA#, 


JulianDate#, 


DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 


RAS(1 TO 


100) 


DECS$(1 TO 100) 


EPOCHS (1 
EQOEE@.( 1 
CelCoord# 
DCM#(1 TO 
TerraCoor 
TerraRec# 
S1#(1 TO 
S2#(1 TO 
N+ ClO” 3 


TO 100) 

TO 2) 

(1 TO 3) 

4 TO 3) 
da#(1 TO 3) 
(1 To: 3) 
>) 

3) 
, 1 TO 3) 


Define error windows 


The quit program confirmation window 


wl. 
Wee 
wl. 
wl 
wl. 
Wale 
Won 
wl 


w1T 
wl1lT 


. bg aBody 
PEGGH IT ChiraG 
.-bgdHighlig 
.fgdaTitle = 
.bgdaTitle = 
,tGdPronpt 
Og al romp t 


action = 1 
edgeline = 
row = 5 


-col = 5 


fgdEdge = 
bgdEdge 
roadBoday 


it hkescu— a) 
ext$(1) = 


if 


CYAN + BRIGHT 
RED 
BRIGHT + WHITE 
RED 
ht = 
ht = 
CYAN 
RED 
YELLOW 
= RED 


O 
0 


WE Oune se roguan: —(Y Otel) 


WEPronpes =)" 


The file save error window 


w2 
w2 


w2 
w2T 
w2T 


fact 1on «= 
.edgeline = 
-row = 5 
neo = 35 
. £fGdEdge 
.bgdEdge 
. fgdBody 
. bgdBody 
Sige fels(veloul ie: 
.bgdHighlig 
.fgdTitle = 
.bgdTitle = 
Jar romp ec 
~Sgal emia 


ities. —— | 
ext$(1) = 


1 


CYAN + BRIGHT 
RED 

BRIGHT + WHITE 
RED 
ht = 0 
ht = 0 

CYAN 
RED 

YELLOW 

RED 


mt Save error 


w2Prompt$ = "Press any key to continue" 


'The file-not-found error window 


ws.action = 1 

w3.edgeline = 1 

w3.row = 5 

wo.col = 5 

w3.fgdEdge = CYAN + BRIGHT 
w3.bgdEdge = RED 
w3.fgdBody = BRIGHT + WHITE 
w3.bgqdBody = RED 


w3.fgdHighlight = 
w3.bgdHighlight = 
w3.fgdTitle = CYAN 
w3.bgdTitle = RED 
w3.fgdPrompt = YELLOW 


O 
0 


w3.bgdPrompt RED 

metaitlies = 

Temexeo( i) = " Error —- File not found or path incorrect " 
w3Prompt$ = "Press any key to continue" 


'The Equinox correction error window 


w4.action = 1 

w4.edgeline = 1 

w4.row = 5 

w4.col = 5 

w4.fgdEdge = CYAN + BRIGHT 
w4.bgqdEdge = RED 
w4.fgdBody = BRIGHT + WHITE 
w4.bgqdBody = RED 


w4.fgdHighlight = 0 
w4.bgdHighlight = 0 
w4.fgdTitle = CYAN 
w4.bgdTitle = RED 
w4.fgdPrompt = YELLOW 
w4.bgdPrompt = RED 


Ware les = 

w4Text$(1) = " Error - Equinox corrections must be a NUMBER 
in seconds " 

w4Prompt$S = "Press any key to continue" 

| 


'The slew confirmation window 
w5.action = 1 

w5.edgeline = 1 

w5.row = 5 

w5.col = 5 

w5.fgdEdge = CYAN + BRIGHT 
w5.bgdEdge = RED 

w5.fgdBody = BRIGHT + WHITE 
w5.bgdBody = RED 
w5.fgdHighlight = 0 
w5.bgdHighlight = 0 
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Woe pet eke 
w5.bgdTitle 
w5.fgdPrompt 
w5.bgdPrompt 


w5Titles = 
w5Texts (1) 
w5Texts (2) 
w5Texts (3) 
+ vf deg. 0 

w5Text$ (4) 
w5Texts (5) 
w5Text$(6) 


TT 


WHITE + BLINK + BRIGHT 
RED 

YELLOW 
= RED 
WARNING 


This action will slew the telescope. 
" The initial zenith angle will be " + Zeniths 
" Any key other than Fl exits this procedure. " 
" The ESC key will halt telescope movement at " 
" any time." 


w5Texts$ (7) es 

w5Prompt$ = "Press the Fl key to slew" 
The Change telescope status window 
w6.action = 1 

w6.edgeline = 1 

w6.row = 5 

w6.cole— 3 

w6.fgdEdge = CYAN + BRIGHT 

w6.bgdEdge = BLUE 

w6.fgdBody = BRIGHT + WHITE 
w6.bgqgdBody = BLUE 

w6.fgdHighlight = 0 

w6.bgdHighlight = 0 

w6.fgdTitle = WHITE + BRIGHT 
w6.bgdTitle = BLUE 

w6.fgdPrompt = YELLOW 

w6.bgdPrompt = BLUE 

w6TitleS = " Change of Telescope Status " 


w6Texts (1) 
w6Texts (2) 
w6Texts (3) 
w6Texts (4) 
w6Texts (5) 
w6Texts (6) 
w6Texts (7) 


w6Texts (8) 
w6Texts (9) 


il 


" Fl: Edit the current jis 

" F2: Choose a star from the current list " 

" F3: Change to a new data file " 

" F4: Change atmospheric parameters " 

"RS: Park Ghe telescopes. 

" ESC: Exit the program and halt the telescope 


Any other key aborts the change " 


w6Prompt$ = "Choose an option" 
f 


'The RA/DEC format error window 


W7 


— 
— 


- row S) 
Pu et @ al 5 
. fGdEdge 
. bgdEdge 
. £fgdBody 


<AcEIoOn = 
.edgeline 


it 


ul 


CYAN, +] BRIGHT 
RED 
BRIGHT + WHITE 


86 


w7.bgdBody = 
w7.fgdHighli 
w7.bgdHighli 
w7.fgdTitle 

w7.bgdTitle 

wee ogarrompt 
w7.bgdPrompt 
w7TitlesS = " 
w7Text$(1) = 
eoes Star ™ 


RED 
ght = 0 
ght = 0 
= CYAN 


ll il 
a 
es 
Oo 


" Error - RA or DEC is formatted incorrectly for 


w7Prompt$ = "Press any key to continue" 
1 

'The slew error window 

w8.action = 1 

w8.edgeline = 1 

w8.row = 5 

w8.col = 10 

w8.fgdEdge = CYAN + BRIGHT 

w8.bgdEdge = RED 

w8.fgdBody = BRIGHT + WHITE 

w8.bgdBody = RED 

w8.fgdHighlight = 0 

w8.bgdHighlight = 0 

w8.fgdTitle = CYAN + BLINK 

w8.bgdTitle = RED 

w8.fgdPrompt = YELLOW 

w8.bgdPrompt = RED 

w8TitleS = " NOTICE " 

w8Text$(1) = " This star is not in view or has a zenith 
angle " 


w8Text$(2) = " 
another star." 
w8PromptS = "Press any key to continue" 
1 


greater than 75 degrees. 


'The zenith angle warning window 


Please choose 


w9.action = 0 

w9.edgeline = 1 

w9.row = 7 

w9.col = 57 

w9.fgdEdge = CYAN + BRIGHT 
w9.bgdEdge = YELLOW 
w9.fgdBody = BRIGHT + WHITE 
w9.bgdBody = YELLOW 
w9.fgdHighlight = 0 
w9.bgdHighlight = 0 
w9.fgdTitle = CYAN + BLINK 
w9.bgdTitle = YELLOW 
w9.fgdPrompt = YELLOW 
w9.bgdPrompt = YELLOW 
w9Title$S = " NOTICE " 
w9Text$(1) = " The zenith " 


87 


" angle is" 


) = STRS (Zenith%) + degrees 1" 


w9Texts (3 

w9Prompts = "" 

{ 

'The zenith angle halt movement window 
wl0O.action = 1 

wl10.edgeline = 1 

w1l0.row = 5 


w10.col = 5 


w10.fgdEdge CYAN + BRIGHT 


w10.bgdEdge = RED 
w10.fgdBody = BRIGHT + WHITE 
w10.bgdBody = RED 


w10.fgdHighlight = 0 
wl10.bgdHighlight = 0 

wl0.fgdTitle = CYAN + BLINK 

wl0.bgdTitle = RED 

w10.fgdPrompt = YELLOW 

w10.bgdPrompt = RED 

wl0Titles = " NOTICE " 

wl0Texts (1) " Telescope movement has been halted. " 
wlOTexts (2) " The zenith angle has exceeded 75 degrees. " 
wlOPrompt$ = "Press any key to continue" 

t 

'The Status help window 

Wii ci@ieeb liga 

wll.edgeline = 1 

Wll.row = 5 
wll1l.col = 3 
wl1.fgdEdge 


CYAN + BRIGHT 


w1l1l.bgdEdge = BLUE 
w11.fgdBody = BRIGHT + WHITE 
wl1l.bgdBody = BLUE 


w1l1.fgdHighlight = 0 

wll bodHighlight — 0 
wll.fgdTitle = WHITE + BRIGHT 
wll.bgdTitle = BLUE 
w11.fgdPrompt = YELLOW 
wl1l.bgdPrompt = BLUE 


wl1Title$S = " Telescope Status Help" 

wllTextS(1) = ™ Slow step: Fast Steger 
f 

w11Text$(2) = " aS 
w11Text$(3) = " up arrow = Increase DEC 8 = Increase 
DEG” 

w11Text$(4) = " down arrow = Decrease DEC 2 = Decrease 
DEC 

w11Text$(5) = " left arrow = Increase RA 4 = Increase 
RA wt 

wllText$(6) =" right arrow = Decrease RA 6 = Decrease 
RA Ww 
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w11Texts$ (7) eee "1 


wlliTexts(8s8) =" Slew: " 

wl1Text$(9) =" " 
wliText$(10) = " 1 = Increase DEC " 
w11Text$(11) = " 3 = Decrease DEC " 
w11Text$(12) = " 7 = Increase RA " 
wl1lText$(13) = " 9 = Decrease RA tL 
wllPrompt$ = " Press any key to exit help " 

1 

'The CTS-10 error window 

wl2.action = 1 

w1l2.edgeline = 1 

wl1l2.row = 5 


tae. cOl = 5 

w12.fgqdEdge CYAN + BRIGHT 
wl2.bgdEdge = RED 
wl2.fgdBody = BRIGHT + WHITE 
w1l2.bgdadBody = RED 
weezer Et GadHighlight = O 
wl2.bgdHighlight = 0 
wl2.fgdTitle = CYAN + BLINK 
w12.bgaTitle = RED 
w1l2.fgdPrompt = YELLOW 
w12.bgdPrompt = RED 


wl2Title$ = "Error" 

wl2Text$(1) = " The CTS-10 WWV Time Syncronization board is 
mew 

w1l2Text$(2) = " responding. Take one of the following 
acierons: " 


Wee hexts(3) =" " 

w12Text$ (4) u a. Check the CTS-10 onboard battery " 
wl2Text$(5) = " b. Manually update the system time to 
UTE." 

w1l2Prompt$ = "Press any key to exit" 

1 


'Initialize the display 
SCREEN 0, , 0, O 

WIDTH 80 
InitialDisplay 

i] 


VKKKKKEKKKEKKEKKEKEKKEKKKEKKEKEKEKEKKKEKEKEEKKEEKRE 


'kx Temporary Debugging Statement ** 
DO 
LOOP UNTIL INKEY$ <> "" 
V KKK KK KKK KK KK KK KKK KK KEK KKKKKKKEKKKKKEKE 
t 
‘Initialize the rest of the system (describe) 
VK KH KEKKKKKKKKEKKKEKKKEKEKEKEEKEEKKEEKKKEKEKKKEKEEKKEKKEKKEKEKKKKEKEKEEE 
KREEKKKEKKKKKEKK 
Vax 
*x* 
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Ik x More Work Needed 
** 


'"POKE (address of external boards), 0 
kk 


"POKE (address of Key paddle buffer), 0 
kk 
PH HKKHKKKKKKKKEKKEKEKKEKKKEKKKEKKKEKKEKKEKKEKKKKKKKKKKKK KKK KKK KKK KKKKKK 


Kak Kea KRKRKRKRKRKEKESK 
’ 


'Check the CTS-10 WWV Time Synch board for response 

year% = VAL(RIGHTS(DATES, 4)) 

IF year%’ < 1992 THEN windows wl2, wi2Text$(), wl2Titles, 
w1l2Prompt$ 

SY ove 

t 

'Perform a one-time calculation of the Equation of Equinox 
correction to 

'sidereal time 

ON ERROR GOTO EQofEQError: 

EquationEquinoxes EQofEQ(), EQofEQUIN! 

ON ERROR GOTO 0O 

t 

'Determine current time and date information for use in 
DayLightSaving, 

'SiderealTimeCalc SUBRoutines and Status DO LOOP, below 

min%’ = VAL(MIDS(TIMES, 4, 2)) 

hour%S = VAL(LEFTS$(TIMES, 2)) 

day% = VAL(MIDS(DATES, 4, 2)) 

monthS = VAL(LEFTS(DATES, 2)) 

year% = VAL(RIGHTS(DATES, 4)) 
t 


'Check whether Daylight Savings time applies to local time 
DayLightSaving DayLightSavings%, hour%, day%, month%, year% 
t 


ChooseDecision% = TRUE 
AtmosChange% = FALSE 
DataEntry% = TRUE 
RADecError% = FALSE 
firstloop% = TRUE 
WindowOn% = FALSE 
Simulation% = FALSE 


DatafileOld$ = "STARFILE. DAT" 

Pressure! = 633! ‘mmHg @ 5000 ft altitude default 
Temperature! = 10! 'c @ 5000 ft altitude default 
Wavelength! = .5 ‘microns @ 5000 ft altitude default 
1 


'Allow the user to change the default atmospheric conditions 
used to calculate refraction 


Atmospheric Pressure!, Temperature!, Wavelength! 
t 


'Perform sidereal time calculations 
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SiderealTimeCalc SiderealTime$, min%z, hour%, day%, month?2, 
years, EQofEQUIN!, GMSTOhrSEC#, JulianDate#, dayFraction#, 
DaysPassed% 

t 


'Align telescope in software 
TeleAlign GMSTOhrSEC#, GHAaries#, TeleRA#, TeleDEC# 
1 


DO 
'Perform sidereal time calculations 
SiderealTimeCalc SiderealTime$, min%, hour%, day%, 
month?, years, BOOEEQUIN(, GMSTOhrSEC#, JulianDate#, 
dayFraction#, DaysPassed% 
t 


te Databmery.)— TRUE THEN 
DataEntry% = FALSE 
'Accept star data entry from disk or by hand 


at 
COECGK WHITE. BLUE, - BLUE 
FOR r% = 14 TO 21 
bL@GALE re 6,0 
PRINT SPACES (68) 
NE KR rs 
EditDecision% = FALSE 
NewFile%t = TRUE 
DO 
COCK WHEE UBLUE. -BLUE 
LOCATE. 17), 110), 0 
PRINT SPACES (45) 
LOCATE: 17, 10,0 
PRINT "Enter star data by: (F1) Disk" 
POCKET 1S, 31.10 
PRINT "(F2) Keyboard" 
kee% = KeyCode%(Character$) 
SELECT CASE kee% 
CASE 15104 'F] 
BadPathType% = FALSE 
NewFile% = FALSE 
FOR 13° =\‘1 70-2 
mOCA UE. 2G. i 
eS, 0 
PRINE 
SPACES (34) 
NEXT 1% 
MOG ny eA 7 ee al 
PRINT DatafileOld$ 
LOCA TR 7 Or al 
INPUT "Please enter 
the path and file name: ", DatafilesS 
IF LEN(Datafile$) = 0 
aN 


ei 


Datafileolds 


= Datafiles 


THEN 
BadPath: 


Datafile$ FOR INPUT AS #1 


O 
BadPathType%s = FALSE THEN 
1 


WHILE NOT EOF(1) 


INPUT #1, Star$(1%), RAS(13), 
13 = 1% + 1 


<j 00 “THEN 


FOR h% = i% TO 100 
Stars(hs) =o" 
RAS(h%) = "" 
DEGS (hs) “=sie 
EPOCH$(h%) = "" 
NEXT h% 
i, ii °C 


SPACES (63) 
ge Oe, © 
"Do you wish to edit the file" 


10, O 


OZ 


DatafileS = 


ELSE 
DatafileOlds 


END ALE 
IF LEN(Datafiles) >0 


ON ERROR GOTO 
OPEN 
ON ERROR GOTO 


LE 


DEC$S(i%), EPOCHS(i%) 


LOOP 
IF1i% 


END Se 
RESEG 
LOCATE 
PRINT 
LOCATE 
PRINT 


LOCATE 


PRINT 
Perior to proceeding? (Y or N) " 


COLOR 
WHITE, BLUE, BLUE 
LOCATE 
ieee 40, 0 
PRINT 
LB Na 
LOCATE 
18, 40, 1 
kee% 
= KeyCode% (Characters) 
1 
UCASES(Character$) = "Y" THEN EditDecision% = TRUE 
Stops% 
= TRUE 
ELSE 
‘no 
action, return to beginning of loop 
END IF 
Boe 
no. aceon: 
return to beginning of loop 
END LF 


CASE oeiGO me eer 2 
EditDecision% = TRUE 
NewFile%t = TRUE 
Stops% = TRUE 
CASE ELSE 
Stops% = FALSE 
END SELECT 
LOOP UNTIL Stops% = TRUE 
BND IF 
i] 
IF EditDecision% = TRUE THEN 
EditDecision% = FALSE 
'Change the screen to display readiness for 
Stam Gata entry 


StarDisplay 

‘Accept star data entry 

StarDataEntry Star$(), RAS ()y DECS () x 
EPOCHS(), NewFile%, wl, wlText$(), wilTitile$, wilPrompt$ 


'Save data to file 
IF NewFile% = TRUE THEN 
NewFile% = FALSE 
COLOR WHITE, BLUE, BLUE 
ROR-r> = 14 TO 21 
LOCATE r%, 6, O 
PRINT SPACES (68) 
NEXT rs 


2 


DataFile2s$ 


= Date abies > 


#1 


DEH@GS(12) ,»@EPOCHS (14) 


ELSE 


DataFile2s$ 


= DataFile2s$ 


+] 


COLOR BLUE EEA Ghe wr ieuiE 
LOCATE 23 2 aa 

PRINT STRINGS(30, 177) 

Datafiles = "StarFilewpar 

COLOR WHITE, BLUE, BLUE 

LOCATE 15, 37, 0 

PRINT Datafiles 

LOCATE ally aeGy an 

INPUT "Please enter the filename: ", 


LOCATE 16, 10, 0 
COLOR CYAN + BLINK, BLUE SBEue 

PRINT "SAVING" 

IF LEN(DataFile2$) > 0 THEN Datafiles 


ON ERROR GOTO NewFileError: 
OPEN Datafile$S FOR OUTPUT AS 


ON ERROR GOTO O 
FOR 1% = 1 TO UBOUND(Star$) 
WRITE #1, Star$(i%), RAS(aiee 


NEXT? 1% 
RESET 


'lf an old file, store torwdicr 
COLOR WHITE, BLUE, BLUE 
FOR r% = 14 TO 21 
LOCATE r%, 6, 0 
PRINT SPACES (68) 
NEXT r% 
COLOR BLUE, BLACK, BLUE 
LOCARE W237 52 eo 
PRINT SERINGS (e057 7) 
COLOR WHITE, BLUE, BLUE 
LOCATE 15°" 37570 
PRINT Datafiles 
LOCATE oes Oe 
INPUT "PleaSe enter the filenames 


LOCATE 16, 10, O 
COLOR CYAN + BLINK, BLUE, BLUE 

PRINT "SAVING" 

IF LEN(DataFile2$) > 0 THEN Datafiles 


ON ERROR GOTO OldFileError: 
OPEN Datafile$ FOR OUTPUT AS 


ON ERROR GOTO O 


IF Filekrrors = FALSE THEN 
FORwuis = I TO 100 
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WRikkE Si Stars (is) / 
RA$(i%), DECS$(i%), EPOCH$(i$) 
NEXT 1% 
BN Date 
RESET 
FOR r% = 14 TO 21 
HOGATE rs, 6, 0 
PRINT SPACES (68) 
NEXT r3% 
END ft 
END IF 


'Allow the user to update the atmospheric parameters 
Siewhich refraction 1s calculated 
IF AtmosChange% = TRUE THEN 
AtmosChange% = FALSE 
Atmospheric Pressure!, Temperature!, 
Wavelength! 
pul Oe se) 
t 


‘Allow the user to choose the active star on which the 


telescope will track 
IF ChooseDecision% = TRUE THEN 


ChooseDecision% = FALSE 

Gneesesctas Sears kas ()e) DECGS() 4 "EPOCHS.( ); 
NewFile%, wl, wlText$(), wlTitleS$, wlPrompt$, currentStar% 

ActiveStar$ = Star$(currentStar%) 


'Calculate Precession and Nutation effects for 
the day 

Precession RAS(), DECS() | EPOCHS(), 
currentStar%, JulianDate#, EpochError%, RAofStar#, DECofStar#, 


RADecError%s 
| 


'RA and DEC a@isplayed are corrected for 
precession only 

RAofStarDisplay# = RAofStar# 

DECofStarDisplay# = DECofStar# 

t 


Nutation RAofStar#, DECofStar#, JulianDate?;, 
Nemo ol#(), S2#() 
i] 


Aberration RAofStar#, DECofStar#, 
RAAberrationCor#, DecAberrationCor#, dayFraction#, DaysPassed% 
RAofStar# = RAofStar# + RAAberrationcCor# 
DECofStar# = DECofStar# + DecAberrationcor# 
t 
Refraction ERmeccune |) Temperature!, 
Wavelength!, GMSTOhrSEC#, HourAngle#, GHAaries#, RAofStar#, 
DECofStar#, RARefractCor#, DecRefractCor#, RAOfMIRA# 
RAofStar# = RAofStar# + RARefractCor# 
DEGotEStar# = DECoOfStar#?#? + DecRefractCcor#? 
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RARefractCorOld# = RARefractCor# 
DecRefractCorOld# = DecRefractCor# 
DO WHILE RAofStar# >= (2# * pi) 
RAofStar# = RAofStar# - (2# * pi) 
LOOP 
DO WHILE RAofStar# < 0 
RAofStar# = RAofStar# + (2# * pi) 
LOOP 
DO WHILE DECofStar# > (pi / 2#) 
DECofStar# = pi - DECofStar# 


LOOP 

DO WHILE DECofStar# < -(pi / 2#) 
DECOLStard — =(pi + BDEGOESUai) 

LOOP 


TeleAngles Zenith%, Azimuth%, GHAaries#, 
RAofStar#, DECofStar#, CelCoord#(), DCM#(), TerraRec#() 
Zenith$ = STRS$(Zenith3) 
w5Text$(3) = " The initial zenith angle will 
be! + -Zeniehs. +" dequass 
Windows w5, w5Text$(), w5Title$, w5Prompt$ 
ESE 
w5.reCuwmncode.,— 504) ri 
END IF 
i] 
‘If the user enters F10 vice Fl, simulation of the 
control process occurs 
IF w5.returnCode = 17408 THEN Filo 
w5.returnCode = 15104' Fil 
Simulation% = TRUE 
ENDSLE 
SELECT CASE w5.returnCode 
CASE 1530425 ra 
WindowsPop 
TimeSec% = VAL(RIGHTS(TIMES, 2)) 
q 
"Display eeleseope Status Secrecy 
TeleStatusDisplay ActiveStar$, 
Simulation’ 
q 
Quits = FALSE 
FirstPass% = TRUE 


DO 
‘Calculate and display local 
sidereal and UTC times 
‘and calculate Refraction 
COrrections 
er TimeSec% = 59 AND 
VAL(RIGHTS (TIMES, 2)) = 0 THEN 


ShowTime minz%, hour%, 
day%, month%, year%, DayLightSavings%, SiderealTime$ 
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TimeSec% = 
VAL (RIGHTS (TIMES, 2) ) 
RefractionPressure!, 
Temperature!, Wavelength! , GMSTOhrSEC#, HourAngle#, GHAaries#, 
RAofStar#, DEG®LSOtar?, RARGCErFactCor#:, Decketractcor?, 
RAofMIRA# 
RAGEStar# = RAOLStar- 
= RARefractCcorOld# 
DECOLrS ear? = 
BeeerStar# - DecRefractCoro0ld# 
RAofStar# = RAofStar# 
+ RARefractCor# 


DEGCOGEStar#F = 
DECofStar# + DecRefractCor# 
RARefractcorold# = 
RARefractCor# 
DecRetractCorold? = 
DecRefractCor# 
Eber ae TimeSec% < 


VAL(RIGHTS(TIMES$, 2)) THEN 
ShowTime min%, hour?, 
day%, month%, year%, DayLightSavings%, SiderealTime$ 
TimeSec% = 
VAL(RIGHTS(TIMES, 2)) 
RefractionPressure!, 
Temperature! , Wavelength! , GMSTOhrSEC#, HourAngle#, GHAaries#, 
RAofStar#, DECofStar#, RARefractCor#, DecRefractCor#, 
RAofMIRA# 
RAofStar# = RAofStar? 
ere tr ractCorOld# 
DECOtTStar? = 
DECoOfStar# - DecRefractCoroOld# 
RAofStar# = RAofStar# 
+ RARefractCor# 


DECoTStar? = 
DECofStar# + DecRefractCor# 
RARefractCorold# = 
RARefractCor# 
DecRefractCorOld# = 
DecRefractCor# 
END et 
DO WHILE RAofStar# >= (27 * 
pi) 
RAofStar# = RAofStar# 
oamee: * pi) 
LOOP 
DO WHILE RAofStar# < O 
RAofStar# = RAofStar# 
feet |6* 6 Pi) 


LOOP 


a | 


DO WHILE DECofStar# > (pi / 


27) 
DECofStarg = pila 
DECGES Eas. 
EOOr 
DO WHILE DECofStar# < “[(jgam 
Za) 
DECofStar# = —=(¢enee 
DEGofStar,) 


EOOr 
t 
'Display current RA, DEC and 
Hour Angle of the star 
COLOR WHITE.” ELUE, VBL 
RAGEStarInHovies? 


RAofStarDisplay# * 24# / (2# * pl) 
RAHOurs#% = 

INT (RAofStarInHours#) 

RAMinutes% 
INT ( (RAofStarInHours# - CDBL(RAHours%)) * 60?) 

RASeconds! = 
( (RAofStarInHours# - CDBL(RAHours%)) * 60# - CDBL(RAMinutes?%) ) 
* 60% 


LOCATE “15> )t0Fae 
PR ie 

RTR IMS Gi Rents, ST R SagR teh ones 2mm ee wos ts 
RTRIMS(LTRIMS(STRS(RAMinutes’))); mop: 
RTRIMS (LTRIMS (STRS(INT(RASeconds! * 100) / 100))); SPACES$(5) 

DECinDeg# = DECofStarDisplay# 
* 3604 7c «Foe, 

DECinDeg! = CSNG(DECinDeg#) 

DECDegrees% = INT(DECinDeg!) 

DECMinutes% = INT( (DECinDeg# 
- CDBL(DECDegrees%)) * 60#) 


DECSecondss® 
CSNG((( (DECinDeg# - CDBL(DECDegrees%) ) x 60#) = 
CDBL(DECMinutes%)) * 60#) 
LOCATE S > 7.0; 0 
P Ree 
RTRIMS (LYRIMS(STRS (DEGCDeGg reese CHRS ( 2 43a 
RTRIMS (LTRIMS ( SMR eb ee iM tmustecis see Be 
RTRIMS (LTRIMS(STRS(INT(DECSeconds! * 100) / 100))); “"t'"; 
SPACES (5) 
I 
'Update the star's hour angle 
and display 
HourAngle# = RAofMIRA# - 
RAofStar# 'in radians 
DO WHILE HourAngle# >= (2# * 


pi) 
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HourAngle# = 
HourAngle# - (2# * pi) 


LOOP 
DO WHILE HomrAngle?: < 0 
HourAngle# = 
hemnrAngle# + (24 * pi) 
LOOP 


HourAngleInDeg# = HourAngle# 
peer / (2% * pi) 
IF HourAngleInDeg# > 180 THEN 
HourAngleInDeg# = 
HourAngleInDeg# - 360# 
END 
HourAngleInHr# = 
HourAngleInDeg# / 15#'Conversion of degrees to hours of time 
HourAngleHours$ = 
STRS$ (FIX (HourAngleInHr#) ) 
IF VAL(HourAngleHourssS) = 0O 
AND HourAngleInDeg# < O THEN HourAngleHours$ = "-0oO" 
HourAngleMin# 
ABS ((HourAngleInHr# - CDBL(FIX(HourAngleInHr#)))) * 60# 
HourAngleMinutes$ 


STRS$ (FIX (HourAngleMin#) ) 
IF LEN(HourAngleMinutes$s) = 1 
THEN HourAngleMinutes$ = "0" + HourAngleMinutes$ 
HourAngleSec# = (HourAngleMin# 
- CDBL(FIX(HourAngleMin#))) * 60# 
HourAngleSecMod% = 
FIX (HourAngleSec# * 100#) 


HourAngleSeconds$ 
STRS (ABS (CSNG (HourAngleSecMod%) / 100!)) 
ROGAIRE, (20 a0 pea0 


[eer S(RTRIMS (HourAngleHourss$) ) ; ta 
Mee LMS(RTRIMS (HourAngleMinutesS$) ) 
LTRIMS (HourAngleSeconds$) ; SPACE$(10) 

t 


> ° 
f e s 


'Calculate and display Zenith 
Angle and Azimuth Angle on status screen 

TeleAngles Zenith%, Azimuth%, 
GHAaries#, RAofStar#, DECofStar#, CelGoorda? (), DCM#(), 
TerraRec# () 

LOCATE 159.37, 0 

PRINT Zenith; 

LOCATE 15, POS(0) - 1, 0 

PRINT CHRS(248); SPACES(5) 

MOCATE .18.,. 36,7 0 

PRINT Azimuth%; 

LOCATE 18, POS(0) - 1, O 

PRINT CHRS$(248); SPACES(5) 

IF RADecError%é = TRUE THEN 


go 


RADecErrors = FABoE 
keet = ENTER 
windows w7, w7TextS(), 
w7Title$, w7Prompts$ 
WindowsPop 
ELSEIF Zenith >= 75 AND 
firstloop% = TRUE THEN 
firstloop% = FALSE 
windows w8, w8Texts$(), 
w8Title$S, w8Prompts 
WindowsPop 
Qulit% = TRUE 
ELSEIF Zenith? >= 75 AND 
firstloop% = FALSE THEN | 
windows wl1O, 
wlOText$(), wlOTitle$, wlOPrompts 
WindowsPop 
Quit% = TRUE 
Pia & 
IF Zenith? > 65 AND 
WindowOn% = FALSE THEN 
WindowOn% = 
TRUE 
ZenithoOld% = 
Zenith% 
w9Text$(3) = 
STRS(Zenith.) +" degrees “ 
windows’ w9, 
w9Text$(), w9Title$S, w9Prompt$ 
END ok 
IF Zenith > 653 
WindowOn% = TRUE THEN 
IF Zeni €ieeee: 
ZenithOld% THEN 
ZenithOld% = Zenith% 
WindowsPop 
w9Text$(3) = STR$(Zenith?) + " degrees " 
windows w9, w9Text$(), w9Title$, w9Prompt$ 
END ik 
jen DE La 
1a Zenith% <> 
Zenithold% AND WindowOn% = TRUE THEN 
Zenithold%s = 
Zenith% 
WindowsPop 
WindowOn% = 
FALSE 
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display boxes 


BLEUE 

oc, 0 
SPACES (1) 
Pes) 6+ 650, O 
SPACES (1) 


position data from encoders and operate 


TeleDEC# in RADIANS 


VK KKKKKKKKKKEKKKKKKKKKKKKKKKKKKK 


JES IDR be os 


“Clear Chemotor state 


COLOR WHITE, BLUE, 


FOR i = 1 TO 2 
LOCATE Ss a. 10 
PRINT 


NEXT if 
DORMS —— she 22 
BOGATE 27)-(10 
PRINT 


NEXT a3 
1 

‘Read 
on them 
'todefine TeleRA# and 


in Telescope 


tkk Read 1n encoder 


data xk 
Ve keKekeke aK KK KK KKK KKK KKK KKKKKKKKKK 
RAGiff# = TeleRA# - 
RAofStar# 
DECdiff# = TeleRA# - 
PEe@orStar# 
TF Fircerasce.-— TRUE 
Ta EN 
She 
ABS (RAGiff#) > .0175 OR ABS(DECdiff#) > .0175 THEN 
SlewCommands Simulation%, TeleRA#, TeleDEC#, RAofStar#, 
DECofStar# 
END Le 
RAdGif f # = 
TeleRA# - RAofStar# 
DECd1f fF = 
TeleRA# - DECofStar# 
are 
ABS (RAdGiff#) < .0175 OR ABS(DECdiff#) < .0175 THEN 
SetCommands Simulation?, TeleRA#, TeleDEC#, RAofStar#, 
BEecorStar# 
ENE sat 


Od. 


FirstPass% 
FALSE 
END IF 
TrackCommands 
kee% 


KeyCodeNowait% (Characters) 
END IF 
SELECT CASE kee% 
CASE. ENTER 
Quits = TRUE 
CAck. ESCAFE 
windows wl, 
w1Text$(), wlTitle$, wlPrompts 
WindowsPop 
i 
wl.returnCode = 121 OR wl.returnCode = 89 THEN Abort% = TRUE 
IF Aborraw. 


TRUE THEN 
PKK KKEKKKKKEKKEKEKEKEKEKEKKEKEKEKKEKEKKKKEKKKEKKEEK 
tak* 
System Shutdown Procedures ** 
tae 
go here!! * 
PKKKEKKKKEKEKEKKEKEKEKKKEKEKKKEKEKKKKKKKKEKKKEEK 
CES 
SYST a 


BND iF 
GASE 1 
windows wll, 
w11Text$(), wl1l1lTitle$, wi11Prompts 
WindowsPop 
CASE LEFTARROW 
‘Change RA 
step multiplier to SLOW speed, INCREASE right ascension 
1 


'POKE (address 


OF RASStep Controller). x 

'Updatemotor 
display box 

COLOR RED? 


BLUE, BLUE 
LOCATE 3 yegn 


0 
PRI 


CHRS (222) 
CASE RIGHTARROW 
' Change RA 
step multiplier to SLOW speed, DECREASE right ascension 
1 
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'POKE (address 
et RA step controller), xX 
'Updatemotor 
display box 
COLOR RED, 


BLUE, BLUE 
LOCATE 3, 20, 


O 
Pea 
CHRS (222) 
CASE NUMLEFTARROW 
'Change RA 
step multiplier to FAST speed, INCREASE right ascension 
i] 


'POKE (address 
Gee Step controller), X 
'Updatemotor 
display box 
COLOR Ree. 
BEUE, BLUE 
EOCAIES 2 0) 
O 
PRINT 
Shino (222) 
CASE NUMRIGHTARROW 
‘Change RA 
step multiplier to FAST speed, DECREASE right ascension 
i] 


'POKE (address 
of RA step controller), X 
'Updatemotor 
display box 
COLOR Risd; 
beor BLUE 
LOCATE 35.20; 
O 
PRINT 
@imen( 2 2 2 ) 
CASE UPARROW 
‘Change DEC 
step multiplier to SLOW speed, INCREASE declination 


'POKE (address 


meer c Step controller), X 

'Updatemotor 
display box 

COLOR RED? 


BLUE, BLUE 
LOCATE 3, 70, 


OQ 


Les 


PRIN? 
Crips 2) 
CASE DOWNARROW 
‘Change DEC 
step multiplier to SLOW speed, DECREASE declination 
{ 
'POKE (address 
Oe DEC Step, controller). xX 
'Updatemotor 
display box 
COLOR RED, 
BLUE, BLUE 
LOCATE 37703 
O 
PRIN 
CHRS (222) 
CASE NUMUPARROW 
‘Change DEC 
step multiplier to FAST speed, INCREASE declination 
t 
'POKE (address 
Of "DEG Seep ceontrolver), xX 
'Updatemotor 
display box 
COLOR RED; 
BLUE, BLUE 
LOCATE 377.3% 
0 
PRIS 
CHRS(222) 
CASE NUMDOWNARROW 
‘Change Bee 
step multiplier to FAST speed, DECREASE declination 
t 
'POKE (address 
Of, DEC "Step controlterm),, = 
‘'Updatemoues 
adieplay box 
COLOR RED; 
BLUE, BEGE 
LOCATE 33a0 
O 
PRiWs 
CHRS (222) 
CASE seven 
‘Manually slew 
to INCREASE right ascension 


‘'Updatemotor 
display box 
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BLUE, BLUE 
0 
CHRS (222) 


of RA slew controller), 


to DECREASE right ascension 


display box 
peer, BLUE 
O 


Gio (222) 


Of RA slew controller), 


eer Ne REASE 


display box 


BLUE, BLUE 
O 
CHRS (222) 


of DEC slew 


memOECREASE 


display box 


BLUE, BLUE 
O 
CHRS (222) 


of DEC slew 


declination 


CONnErOl ler). 


deciination 


eomero ler). 


x 


X 


COLOR RED, 
LOGATE 3105 


PRN 


"POKE (address 


CASE nine 
"Manually slew 


‘'Updatemotor 
COLOR RED, 
LOCATE 3.107 
PRINT 

'POKE (address 


CASE one 
‘Manually slew 


'Updatemotor 
COLOR’ RED, 
LOCALE. 3: 66, 
PRINT 

'POKE (address 


CASE three 
‘Manually slew 


‘'Updatemotor 
COLOR RED 
LOCALE 34760; 

PREN. 


'POKE (address 


IOS 


GASE Elst 
AR lic slew 


motors OFF once no key is pressed 
J 


"POKE (address 
of RA slew controller), X 
'POKE (address 
of DEC slew controller), X 
END SELECG@E 
LOOP UNTIL Quit% = TRUE 
IF WindowOn% = TRUE THEN 
WindowOn% = FALSE 


Windows Pop 
EN De eer 
Quits = FALSE 
Windows’7 w6, w6Texts$(), w6Titles, 
w6Prompt$ 
WindowsPop 


SELECT CASE wé6é.returnCode% 
CASE 125104 “srt 
EditDecision%’ = TRUE 
ChooseDecision% = TRUE 
firstloop% = TRUE 
CASE 5c 60 SEZ 
ChooseDecision% = TRUE 
firstloop% = TRUE 
CASE 156067" F3 
DataEntry% = TRUE 
ChooseDecision%’ = TRUE 
firstloop% = TRUE 
CASE 15872 ' F4 
AtmosChange% = TRUE 
ChooseDecision% = TRUE 
firstloop% = TRUE 
CASE 16123 = S55 


i, 2 a> at db db a> a> at a> a> a> db a> a> a> dt at db dy dy dt a db a> a> a> ab 
Peeks Park the 


telescope xx 


DKK KKK KKEHKKEKKKKKKKKKKKKKKKKKSE 
CASE ESCAPE 
windows wl, wlTextS$(), 

wlTitle$, wlPrompt$ 

Windows Pop 

IFwl.returnCode=121 
OR wl.returnCode = 89 THEN Abort% = TRUE 

IF Abort% = TRUE THEN 


PK KKKKKKEKKEKEKKKKKEKEKKEKEEKEKEKEKKEKEKEKEKEEE 


ROG 


kK System 
Shutdown Procedures ** 


kk 
go here!! * 
PKK KKK KKKKKKKKKKKEKKEKKKEKKKKKKKKKKKKEK 
crs 
SvoluEM 
ele JDe Ane 
CASH OBESE 


END Spike CYT 
t 
CASE Elio Er: 
WindowsPop 
windows w6, w6Text$(), w6Titles, 
w6Prompts 
WindowsPop 
SELECT CASE w6é.returnCode% 
CASE 15104 Ei 


Edatbeeheion. = “PRUE 
ChooseDecision% = TRUE 

CASH 15260 " .F2 
ChooseDecision% = TRUE 


CASE 156216 "° a3 
Databnitie.e— TRUE 
ChooseDecision% = TRUE 
CASESU58 72.) 2na 
AtmosChange% = TRUE 
ChooseDecision% = TRUE 
CASE 1626" "E5 


VK KKK KKK KKK KKK KK KKK KK KKKKKKE 
kx Park the 
telescope x 


VKKKKEKKKEKKEKKEKKKKKKEKKKKKKKEKEEEEK 
CASE 2S GAPE 
windows wl, wlTexts$(), 


wlTitle$, wlPrompts$ 
WindowsPop 
IFwl.returnCode=121 


Okwwl-returnCode = 89 THEN Abort% = TRUE 
IF AbDorts = TRUE THEN 


PKK KKKKKEKEKKKKKKKKKEKKKKKKKKKKKKKKKKEK 
1k System 


Shutdown Procedures ** 
kek 


go here!! * 


kK KK KK KKK KK KKK KKK KKK KKKKKKKKKKKKEKSEK 


G7 


CLS 
SYolEM 
PND aa 
CASE ELSE 
ENDOSELECY 
END SELECT 


LOOP 

4 

'Error Handling and time updates 

GOTO EndLabel: 

{ 

NewFileError: 
RESET 
BEEP 
windows w2, w2Text$(), w2Title$, w2Prompts$ 
WindowsPop 
FileError% = TRUE 
RESUME NEXT 

{ 

OldFileError: 
REoEYT 
BEEP 
windows w2, w2TextS$(), w2Title$, w2Prompts 
WindowsPop 
FileError% = TRUE 
RESUME ShiE cr 

{ 

BadPath: 
RESE® 
BadPathType% = TRUE 
BEEP 
Windows w3, w3Text$(), w3Title$, w3Prompts 
WindowsPop 
RESUME NEAT 

{ 

EROGCrFECE LEC ER: 
BE te 
windows w4, w4Text$(), w4Title$, w4Prompts 
WindowsPop 
REOOME 

{ 

EndLabel: 

END 


PKKKKKKKKEKKKKEK KEK KKK KKK KKK KKK KKK KKKKKKEKKE 


rls Name: Aberration xx 
'kx Type: Subroutine x x 
'k* Module: MIRACTRL. BAS xx 
'k* Language: OQureckbasice 4.50 * 


PKKEKKKKKKKKKKEKKEKKEKKEKKEKKKKKEKEKKKKKKEKKKKKER 
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| 
fealeculates the aberration correction fo) celestial 


coordinates. 
1 


'Example one use: Aberration RA@fStar ys -, DEGOLS tans, 
RAAberrationCor#, DecAberrationCor#, dayFraction#, DaysPassed% 
1 

SUB Aberration (RAofStar#, DECofStar#, RAAberrationcCor#, 
DecAberrationCor#, dayFraction#, DaysPassed%) 

1 

'Calculate the mean longitude of the Sun, L# 

d# = dayFraction# + CDBL(DaysPassed@%) 


L# = 298.926 + .985647 * d# ' in degrees 

ie iy * (2% * pi) / 3607 ' Convert to radians 

t 

'Define the constant of aberration, k# 

k# = 20.496 ' in arc seconds 

k# = k# * (2# * pi) / (360# * 3600#)' convert to radians 
1 

awetrime the obliquity of the ecliptic 

epsilon# = 23.440332# - .00000036# * d# 'degrees 
epSilon# = epsilon# * (2# * pi) / 360# 'radians 


| 
Teme wlate the corrections to RA and Declination due to 
aberration 


RAAberrationCor# = (-k# * SIN(L#) * SIN(RAofStar#) - k# x 
COS(L#) * COS(epsilon#) * COS(RAofStar#)) / COS(DECofStar#) 
DecAberrationCor# = -k# * SIN(L#) * COS(RAofStar#) * 
SIN (DECofStar#) + k# * COS) * (COS (epsilon#) * 
SIN(RAofStar#) * SIN (DECofStar#) - SIN (epsilon#) x 
SeamphPcofStar#) ) 

i] 
END SUB 

PKKKKKKKKKKKEKKEKKKKKKEKKEKKKEKKKEKKKKKKKKKKEKK 

Vk x Name: Atmospheric * 

tx Type: Subroutine x x 

lk x Module: MIRACTRL. BAS kx 

1k Language: Quickbasic 4.50 xx 


PK KKKKEKKKKKKKKKKEKKEKKEKEKKKEKKKKKKKKKEKEKKKEEKSE 
1 
'Allows user entry of atmospheric temperature, pressure and 


the center 
'wavelength of light for an observation. These data are used 


oeeetraction 


"camculations. 

1 

‘Example oof § use: Atmospheric Pressure!, Temperature!, 
Wavelength! 

1 

SUB Atmospheric (Pressure!, Temperature!, Wavelength! ) 


NO Dg, 


‘Clear the body of the status box 

COLOR WHE SoU Sebi 

FOR (bs) =—14 TOR21 

E@CATLEMG 4... Grmae 
PRINT SPACES(68) ; 

NEXT r% 

t 

LOCATE ela > 5) 0 

PRINT Temperature! 

LOCATE io, © 

INPUT "Please enter the atmospheric temperature (C): ", 
NewTemp$ 

IF LEN(NewTempS) > 0 THEN Temperature! = VAL(NewTemp$S) 
LOCATE Dire 55., 0 

PRINT Pressure! 

LOCATEMIY . eherane 

INPUT "Please enter the atmospheric pressure (mmHg): ", 
NewPres$ 

IF LEN(NewPres$) > 0O THEN Pressure! = VAL(NewPres$) 
LOCATE 20 -4igr.0 

PRINT Wavelength! 

LOCATE 1275 10, 0 

PRINT "Please enter the center wavelength of light" 
LOCATE 20, 10, O 


INPUT "for the chosen star (microns): ", NewWaves$ 

IF LEN(NewWaves) > O THEN Wavelength! = VAL(NewWave$S) 
t 
ENDVsSvs 


PAKKKKKKEKKKKKKKKKKEKKKKKKEKKKEKKKEKKKEKKKKKKEKSE 


Va Name: ChooseStar sts 
ex AEN) Sy=)5 Subroutine kx 
Vk Module: MIRACTRL. BAS aK 
tax Language: Qulekbasic 450 ws 


Ve KH KKK KKK KKK KKK KKK KEK KKK KKK KKK 

| 

'Allows choice of stars from list. Sets active star focus 
and "track 


‘of telescope. 
’ 


‘Example of use: StarDataEntry ~Star5(), RAS(), DEGlae 
Epochs(), NewFile%, wl, wlText$(), wlTitleS, wlPremere 
CurrentStar% 

‘Parameters: (none) 

'Variables: Staged) A 5-element matrix of desired 
Stars tei 

viewing 

RAS () A 5-element matrix of right 


ascensions 
for the stars listed in Star$() 


1BEo 


‘ DECS () A 5-element matrix of declinations 
EOonetme, stars 1isceed in Stars () 
Epoch$() A 5-element matrix of epochs for 
the 
RA's and DEC's of the stars listed 
ime stars () 
g%, h%, 1% Looping indices 
' A flag activated when the 
Backspace key is hit 

: 1% A flag activated when the Shift 


LJ 
oe 


hase Up or Down 
Arrows are hit 
: row%, col%é The current row and column markers 


fom the 

: cursor 

StOpS A string representing whether or 
not the 

character entry loop is to be 
exited 


' kee% A unique code returned by the 
KeyCode% or 

KeyCodeNoWait% Functions for each 
key hit 

: Character$ A string read from the keyboard 
buffer 

Currents A string representing the entry of 
iEnvemcurrent 

field. This variable will become 
one of the 

: Eel lewsndg: Star$S(), RAS(), 
DECS(), EPOCHS() 

t 


'Module Level 

BEGEARATIONS: DECLARE SUB StarDataEntry (StarS(), RAS(), 
DEG () , 
Epochs(), NewFile%, wl, 
wlText$(), wlTitle$, wlPrompt$, CurrentStar%) 

t 


SUB ChooseStar (Star$(), RA$(), DECS(), EPOCHS(), NewFile%, wl 
AS WindowsType, wlText$(), wlTitle$, wlPrompt$, currentStar%) 
i] 

‘Initialize variables 

current$ = Star$(1) 


pews = 17 
eos = 9 


PprestRows = 17 
bastRows = 21 
LastCol% = 62 


ple Lines = 1 
TopFileLineShown% = 1 
currentStar% = 1 


Va 


ChangeFileLine% = FALSE 
HomeOrEnd% = FALSE 
Stops% = FALSE 
QuitEntry% = FALSE 

1 


‘Erase Box 
COLOR WHITE, BEUE] SBLUE 
FOR rs = 14a8© 21 
LOCiIRnmr = Grane 
PRINT SPACES(68); 
NEXT rz 
1 
'Draw bottom edge of the status box 
LOGATRER22 eS. 0 
PRINT CHRS$(200); STRINGS$(68, 205); CHR$(188); 
t 
'Place lettering 
COLOR WHITE, BLUES SeEvVE 
LOGATEe4 227. 0 
PRINT "Choose a star to view and press F10." 
LOCATE 15, 33, 0 


PRINT "RA DECI 
LOCATE IG, 110. 0 
PRINT "Star Name hie Mon Sec 


RGCAT EG NlaS nO 
PRINT CHRS(248 ) acto a: 

LOCATE 16,63, G 

PRE Ee Oe. 7, 

t 

'Complete data fields by reading in one key at a time, ESC 
exits 

DO 


| 
rte 


IF TopFileLineShown% < 1 THEN TopFileLineShown% 
IF TopFileLineShown% > 96 THEN TopFileLineShown% = 
IF FileLine% < 1 THEN FileLine% = 1 

IF FileLine% > 100 THEN FileLine% = 100 


| 
\O 
rep) 


IF row% < FirstRow% THEN row% = FirstRow% 
IF row% > LastRow% THEN row% = LastRow% 


'Update the row numbers and fields 
FoR ag: — lero 5 
COLOR WHITE) SLUE  SEUE 
LOCATE 16 + g%, 6, O 
PRINT TopFileLineShown% + g% - 1 
COLOR WHITE.) SBEACK,) SEvUE 
EOGATE 6c: son mG 
PRINT Star$(TopFileLineShown + g% - Wij 
SPACES(19 — LEN(Stars(toprilebineshown. 150 45- on 
LOCATE 16 + g%, 33, 0 
PRINT RAS(TopFileLineShown%’ + g% - 1) + 
SPACE$(10 - LEN(RAS(TopFileLineShown% + g% - 1))) 
LOCATE, 167+ 935 46). 00 


eZ 


Pia wi DECS(TODribebeneshowns + Gs = 1) + 
SPACE$(10 - LEN(DECS$(TopFileLineShown% + g% - 1))) 
LOCATE 16 + g%, 63, 0 


PRINT EPOCHS(TopFileLineShown%’ + g% - 1) + 
SPACE$S(10 - LEN(EPOCHS(TopFileLineShown% + g% - 1))) 


NEXT g% 
COLOR BEAGCK, WHITE, BLUE 
LOCATE row%, 10, O 
PRINT Star$(FileLine?) + SPACES(19 - 
Pe Star$S(FileLine?) ) ) 
ROGATESGOW>s, 33, 0 
PRINT RAS (FileLine$) 0 SPACE$(10 - 
LEN (RAS (FileLine$) ) ) 
LOCATE row%, 48, 0O 
PRINT DRGs (Eibeh ines) + SPACES$(10 - 
LEN (DECS$ (FileLine3) )) 
LOCATE rows, 63, O 
PRINT EPOCHS(FileLine%) + SPACES (10 - 
LEN (EPOCHS (FileLine$) ) ) 
' 


COLOR WHITE, BLACK, BLUE 
'Read the next character and operate on it 
DO 

RCGAUESGews,. COls + 1, 1 

kees = KeyCode%(Character$) 

SELECT CASE Kee% 

CGRSera ENTER 
IF row% = LastRow% AND col% 


Basecols THEN 


row% = row% + 1 
IF row% = LastRow% + 


eer eis N 


TopFileLineShown% = TopFileLineShown% + 1 
FileLines = 
FileLine% + 1 
END TE 
ELSE 
mo act lon 
END iF 
Stops% = TRUE 
CASE ESCAPE 
windows wl, w1Texts(), 
wlTitle$, wlPrompts$ 
ifs wt. returncode = .12)]. (OR 
wl.returnCode = 89 THEN Abort%é = TRUE 
IF Abort% = TRUE THEN 


DK KKK KKK KKK KKK KKK KK KKK KKKKKKKKKKKKK 
'kx System Shutdown 
Procedures ** 


das 


"kx go here!! 
xx 


UHKKKEKKKKKEKKKEKKEKKKEKKEKEKEKKEKKEKEKEKKKKEEEEKSE 


WindowsPop 

CLs 

S Youve 
ladys 

WindowsPop 
BNE Sie 


CASE UPARROW 
IF row% = FirstRow% THEN 
TopFileLineShown% = 
TopFileLineShown% - 1 
BND. IF 
row% = row% - 1 
ChangeFileLine% = MINUS 
Stops% = TRUE 
CASE DOWNARROW 
IF row% = LastRow% THEN 
TopFileLineShown% = 
TopFileLineShown% + 1 
END IF 
row% = row% + 1 
ChangeFileLine% = PLUS 
Stops% = TRUE 
CASE PGDOWN 
IF TopFileLineShown% > 95 THEN 
ChangeFileLine% = PLUS 
*x* (96 - TopFileLineShown$) 
TopFileLineShown% = 96 


Rigs 
TopFileLineShown% = 
TopFileLineShown% + 5 
ChangeFileLine% = PLUS 
<a 
BNDeoiae 
Stops% = TRUE 
CASE SPGuUFr 
IF TopFileLineShown% < 6 THEN 
ChangeFileLine% = 
MINUS * (TopFileLineShown% - 1) 
TopFileLineShown% = 1 
Pile & 
TopFileLineShown% = 
TopFileLineShown% - 5 
ChangeFileLine% = 
MINUS 4 —5 
BNDeerr 
Stops% = TRUE 
CASE HOME 


114 


TopFileLineShown% 

HomeOrEnd% = 1 

Stops... — TRUE 
CASE ENDET EE 

TopFileLineShown% 

HomeOrEnd% = 100 

Stops% = TRUE 
CASE 174065 7710 


lI 
a) 


96 


Guintemestars =| Filebines 


QuitEntry% = TRUE 
Stops% = TRUE 
CASE. ELSE 


PKK KEKKKKKEKKKKKKKEKKKKEKEKEKKRKKKEKKKKKKKKKKKKEEK 
'kx* 
*k* 


Ik xk Error Handling 


Needed kx 
tex 


PKK KEKKKKKKEKKKEKEKKKKKKEKEKKKEKKEKKEKEKEKKEKKKKKKKEKE 


i BY Ses) yl liy sh Ci & 
LOOP UNTIL Stops% = TRUE 
Stops% = FALSE 
IF FileLine% < 1 THEN FileLine% = 
IF FileLine% > 100 THEN FileLine% 
IF ChangeFileLine% <> FALSE THEN 


LOO 


FileLine%’ = FileLine% + ChangeFileLine% 


IF FileLine% < 1 THEN FileLine% = 


IF FileLine% > 100 THEN FileLine? 
ChangeFileLine% = FALSE 
currentS = Star$(FileLine%) 
END IF 
| 
IF HomeOrEnd% <> FALSE THEN 
IF HomeOrEnd% = 1 THEN 
row% = FirstRow% 
ELSEIF HomeOrEnd% = 100 THEN 
rows = LastRow% 
END IF 
FileLine%’ = HomeOrEnd% 
HomeOrEnd% = FALSE 
COLOR WirlLE ~ BEUE. .edasue 
LOCATE LastRow%, 9, 0O 
PRINT SPACES (1) 
current$ = Star$(FileLine%) 
END IF 


LOOP UNTIL QuitEntry% = TRUE 


1S 


J 
= 100 


END SUB 


PKK KKKKEKKKKE KEK KKK KKK KKKEKEKKEKKKEKKKKKKKKKKEKK 


tk * Name: DayLightSaving x x 
"kx Type: Subroutine x * 
kx Module: MIRACTRL. BAS x * 
tk x Language: Quickbasic 4.50 * x 


PKKKKKKKEKEKKKEKEKKKEKKEKEKEKEKEKEKEKKEKKKEKEKEKEKKKKEKEKK 
t 


‘Determine if Daylight Savings Time applies to local time 
{ 


'Example of use: DayLightSaving DayLightSavings%, hour%, 
day%, month%, year% 

‘Parameters: (none) 

'Variables: DayLightSavings% A varaible representing 


whether or not 

daylight savings time is 
im ef freer 

GdayOfWeek$ A string representing the day of 
the week 

SundayPassed% A varialbe representing 
whether or not 

Sunday has passed in the 


current week 


’ hewrs< The current hour of teime 

day% The current day 

: month% The current month 

years The current year 

t 

'Module Level 

t 

' DECLARATIONS: DECLARE SUB DayLightSaving 
( 

t 


DayLightSavings%, hour%, day%, month%, year%) 


SUB DayLightSaving (DayLightSavings%, hour%, day%, month%, 
years) 
PKK KKEKKKK KK KKK KK KKK KEKEKKEKKEKKKEKKEK KKK KKKEKKKKKEKKEKKEKEKHK 
'xx Determine the day of the week (=dayOfWeekS) ** 
'x* and if Sunday has passed (=SundayPassed$) x x 
PKK KKK KKK KKK KKK KKK KKK KKK KKK KKKKEKEKEKHEKEKEKEKKKK KKK 


IF month% > 10 OR month% < 4 THEN DayLightSavings% = 


FALSE 
IF month% > 4 AND month% < 10 THEN DayLightSavings% = 
ts 
LF months = 4. eh 
LE dave = 25) Set 
IF dayOfWeek$S = "SUNDAY" THEN 
IF hours > 1 THE 
DayLightSavings% = TRUE 
ELSEIF SundayPassed%* = TRUE THEN 


DayLightSavings% = TRUE 


iG 





END iF 
END IF 
END IF 
IF months = 10 THEN 
IF day% > 24 THEN 


IF dayOfWeek$S = "SUNDAY" THEN 
ie hours > 1 TREN 
DayLightSavings% = TRUE 
ELSEIF SundayPassed%® = TRUE THEN 
DayLightSavings% = TRUE 
END DE 
END@a Fr 
BND IF 
END SUB 
VKKKKKEKEKKEKEKEKEKEKKEKEKEKEKKEKKEKEKKEKEKEKEKKEKKEKKKKKKEEE 
Vk * Name: EquationEquinoxes ** 
Vk & ioe: Subroutine * 
'kx* Module: MIRACTRL. BAS ** 
Vx Language: Quickbasic 4.50 * 


PKKKEKKEKKKEKKEKEKKEKEKKEKKKEKKKKEKEKKKKKKKKKKKKEKE 

| 

'Allows input of time correction argument from the Equation of 
Equinoxes 

'(See page B6-B15, Astronomical Almanac) (Correction needed due 
permutation 

pote th ) 

i] 


'Example of use: EquationEquinoxes EQofEQ(), EQofEQUIN! 
"Parameters: (none) 

'Variables: EQofEQ() An array of the the equation of 
equinox 
results for the next two days 
EQofEQUIN! The equation of equinox 
Cemmeecti1on for 

the current time 

NextDay$ A string representing either 
the current 

Or next day 

NextDate$ A string representing either 
the current 

’ or the next total date (month, 
day, year) 

| 


sroaule Level 

DECLARATIONS: DECLARE SUB EquationEquinoxes (EQofEQ(), 
EQofEQUIN! ) 

i] 


SUB EquationEquinoxes (EQofEQ(), EQoOfEQUIN! ) 


'Input the equation of eguinox corrections for the 
Current day 


UY, 


‘and the next day 
COLOR WHITE, BLUE, BLUE 
FOR 1% = 1 TO 2 
NextDay$ = LTRIMS (STRS (VAL(MIDS (DATES, 473uz5m 


NextDate$S = LEFTS(DATES, 3) + NextDay cam 
REGHTS (DARE S oy 

LOGATE 16 + 1%;,8573—0 

PRINT EQofEQ(1%) 

LOCATE 16 + 1%, 10, O 

PRINT "Enter the Equation of Equinoxes for "; 


NexzEDatresi 6 "s 
POGARE 01707 hoo nO 
PRINT "(in seconds)" 
LOCATE SG soe 
INPUT "") EFQOfEQ(i3) 
NEXT 1% 


'Interpolate to find the equation of equinox 
COrrect 1on for the 

‘current time 

EQofEQUIN! = (TIMER * (EQofEQ(2) - EQofEQ(1))y 
86400!) + EQofEQ(1) 

i] 


'Erase the above cues 
FOR 1% = 1 TO 4 
LOGATE? 16 + 12722 eo 
PRINT SPACES (45) 
NEXT 13% 
END SUB 


PARK KKKKEKKKEKKEKKKKKKKEKKKEKKKKKKKKEKKKEKKKKEKEKSE 


tks Name: InitialDisplay * 
tke Type: Subroutine x 
kx Module: MIRACTRL. BAS ** 
lk Language: Quickbasic 4.50 * 


VKKKKKKEKKKEKKEKEKEKEKEKEKKKEKEKEKEKEKKKKKEKKKKKKKKKEKSK 
' 


'Creates the initial display screen graphics 
q 


'Example of use: InitialDisplay 
‘Parameters: (none) 
'Variables: LS jeeias Looping indices 


'Module Level 

DECLARATIONS: DECLARE SUB InitialDisplay 
1 

SUB Dnt taliban plan 

CLs 

‘Color the baeckaround 

FOR 13 = 1] T0725 


ITs 


COLOR BLUE, BLACK, BLUE 
PRINT) STRENGS (80, 177) 
NEXT 1% 
LOCATE 24, 1, O 
PRINT STRINGS(80, 177); 
MOGATE 25, 1, 0 
PRINT STRINGS$(80, 177); 
t 
Betace the lettering 
MOCATE 4, 18, O 
COLOR WHITE, BLUE, BLUE 
PRINT "Monterey Institute for Research in Astronomy" 
LOCATE 6, 28, 0 
PRINT "Oliver Observing Station" 
HOCATE 10, 30, O 
PRINT "Automated Telescope" 
MOCHTE 11, 32, 0 
PRINT "Control PacKage" 
t 
Dimaw top edge of status box 
BESATE 13, 5, O 
COLOR WHITE, BLUE, BLUE 
PRINT CHR$(201); STRINGS(68, 205); CHR$(187); 
t 


'Draw the body of the status box 
Boke rs = 14 TO 21 
BocATE rz, 5,0 
ERENT GHRS (186); 
PRINT SPACES(68) ; 
Brand "CHRS (186) ; 
NEXT r% 
1 
‘Draw bottom edge of the status box 
WeekTE 22, 5, O 
PRINT CHRS$(200); STRINGS(68, 205); CHRS$(188); 
t 
iesolay "Initializing" status in box 
SenorR WHITE + BLINK, BLUE, BLUE 
POeSTE 17, 25, O 
Pree? “Initializing" 
meenTE 17, 37, 0 
SeerokR WHITE, BLUE, BLUE 


PRINT " - Please stand by" 
END SUB 

PK KKK KEKE KK KKEKKKKEKKKEKKEKKEKKKEKKKKKKKKKKKKKEKK 
kk Name: KeyCode% * 
tk* Type: Funet ron * 
kx Module: MIRACTRL. BAS * 
"kk Language: Quickbasic 4.50 xx 


KKK KKK KK KKEKEKKKKKKKKKKKKKKKEKKKK KKK KKKKEK 


Eo 


'Returns a unique integer for any key pressed 
q 


'Example of Use: kee%t = KeyCode% Character$ 
‘Parameters: (none) 
'Variables: Characters A string representing the 


current Keyboard 
buffer value 
{ 
‘Module Level 
DECLARATIONS: DECLARE FUNCTION KeyCode% (Character$) 
' 
FUNCTION KeyCode% (Character$) STATIC 

DO 

Characters = INKEYS 
LOOP UNTIL Characters <— 
KeyCode% = CVI(Character$ + CHRS$(0) ) 

ENDREFUNCT ION 


Pk KKK KKK KE KKK KKKKKKKKKKKKKEKKKKKKKKKKKKKKE 


tk x Name: KeyCodeNoWait% * x 
= Type: FunecETen x 
tk x Module: MIRACTRL. BAS * 
kx Language: Ourckbasie. 4450 * x 


PKK KKKKKKKEK KEKE KEKKKEKKKEKKKEKEKEKEKEKEKKEKKKKEKKEEEE 
' 


'Returns a unique integer for any key pressed, but does not 
wait 

'for a key to be pressed 

q 


'Example of Use: kee% = KeyCodeNoWait% Character$ 
‘Parameters: (none) 
‘Variables: (none) 
‘Module Level 
DECLARATIONS: DECLARE FUNCTION KeyCode% (Character$) 
q 
FUNCTION KeyCodeNoWait% (Character$) STATIC 
Character$ = INKEYS 


IF Characters. <>. U"  THEN 

KeyCodeNoWait%’ = CVI(Character$ + CHR$(0) ) 
Ee oie: 

KeyCodeNowWait% = 0 
END DE 


END FUNCTION 


PKK K KKK KKK KEK KKK KKK KKK KKK KEKEKEKKKEKKKKKKEKK 


Ik x Name: Nutation * x 
1k x Type: Subroutine * 
none Module: MIRACTRL. BAS xx 
kx Language: Quickbasic 4.50 * x 


PAKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKE 


Jup2i0: 


feaenlates the nutation corrections to celestial coordinates 
Perea Given 


'day and time. 
{ 


'Example of use: Nutation RAofStar#, DECofStar#, JulianDate#, 
N#(), S1#(), S2#() 
t 


'Note: Astronomical Almanac notation is used wherever possible 
t 


SUB Nutation (RAofStar#, DECofStar#, JulianDate#, N#(), S1#(), 
S2#()) 
{ 


‘Unmodified celestial coordinates converted to rectangular 
form 


eae (1) = COS(DECofStar#) * COS(RAofStar#) 
S1#(2) = COS(DECofStar#) * SIN(RAofStar#) 
—eo EM DECOLSLaAr se ) 


Sa) 


'Nutation Parameters 
Gee— JUulianDate# - 2447871. 5% 
epsilon# = 23.440332# * (2# * pi) / 360# 


DeltaPsi# = -.0048#4 * SIN((279.9# - .053# * d#) * (2# * pi) 
/ 360#) - .0004# * SIN((197.9# + 1.9714 * a#) * (2# * pi) / 
3604) 


DeltaPsi# = DeltaPsi# * (2# * pi) / 360# 

{ 

DeltaEpsilon# = .0026# * COS((279.9#4# - .053# * d#) * (2# * 
mpm es60#) + .0002%# * CcOS((197.94 + 1.971# * d#) * (24 * pi) 
/ 360#) 

DeltaEpsilon# = DeltaEpsilon# * (2# * pi) / 3603 
’ 


'Nutation Rotation Matrix 


fee 1) = 1 
N#(1, 2) = -DeltaPsi * COS(epsilon?) 
N#(1, 3) = -DeltaPsi * SIN(epsilon#) 
N#(2, 1) = DeltaPsi * COS(epsilon#) 
N#(2, 2) = 1 
N#(2, 3) = -DeltaEpsilon# 
N#(3, 1) = DeltaPsi * SIN(epsilon#) 
N#(3, 2) = DeltaEpsilon# 

= l 


N#(3, 3) 


'Celestial Coordinates modified for nutation in rectangular 
form 

eeeael)) = N#(1, 1) * S1#(1) + N#(1, 2) * S1#(2) + N#(1, 3) * 
Siz 3) 
S24(2) 
S1# (3) 


Pee ee om jet 2 > 2) -* SlA(2) + N#(2, 3) # 


Pan 


S2#(3) = N4#(3, 1) * S1#(1) + N#(3, 2) * S1#(2) + N#(G eee 
S1#(3) 


'Decode S1l1#() elements from rectangular form to RA and DEC 
IF S2#(1) >= 0# AND S2#(2) < O# THEN 


alpha# = ATN(=S2#(1) / S2#(2)) + (3% * pi) 7 22 
ELSEIF S2#(1) < O# AND S2#(2) < O# THEN 

alpha# = ATN(S2#(2) / S2#(1)) + pi 
ELSEIF S2#(1) < O# AND S2#(2) >= O# THEN 

alpha# = ATN(-S2#(1) / S2#(2)) + pi / 2# 


ELSE alpha# = ATN(S2#(2) / S2#(1)) 
END IF 


delta# = ATN(S2#(3) * COS(alpha#) / S2#(1)) 


RAofStar# = alpha# 
DECofStar# = delta# 
1 


END’ SUE 

VK HKKKKKKKEKKKKKEKKKRKKKKKKKR KKK KKK KKK KKK 
"kx Name: Precession kx 
Vk Type: Subroutine * x 
'k* Module: MIRACTRL. BAS kk 
"kx Language: Onieckbascic 4750 *& 


PHKKKKKEKEKEKKEKEKEKKKEKEKKKKKEKKKKKKKKKKKKKKKKKEK 
' 


'Calculates the precession corrections ES celestial 
coordinates for a given 

'day and time. 

1 

'Example of § use: Precession RAS(), DEESiO), EPOCHS(), 
currentStar%, JulianDate#, EpochError%, RAofStar#, DECofStar#, 
RADecError’% 


'Note: Astronomical Almanac notation is used wherever possible 
i] 


SUB Precession (RAS$(), DECS$(), EPOCHS$S(), currentStar%, 
JulianDate#, EpochError%, RAofStar#, DECofStar#, RADecError?%) 
i] 


'Decode RAS() and DECS() into radians 


RAFirstColons = INSTR(RAS (ecUrrentStar a3.) 
RASecondColon% = INSTR ( (RAF@rSsteolon-. + ae 


RAS (currentStar%), “:") 
IF RAFirstColon% = 0 OR RASecondColons = 0 THEN 
RADecErrors = 


TRUE 
Ries E 
RAofStar# = VAL(LEFTS (RAS (currentStancoe 
(RAFirstColon%’ - 1))) * 3600 + VAL(MIDS(RAS(currentStar%), 
(RAFirstColon% + 1), 2)) * 60 + VAL(RIGHTS (RAS (currentStar%), 
(LEN (RAS (CurrentStar%)) - RASecondColon$) ) ) 


Zc 


RAofStar# = RAofStar# * 2# * pi / 86164.090537 
'Seconds of time converted to radians of sidereal movement 
DO WHILE RAofStar# >= (2# * pi) 
RAofStar# = RAofStar# - (2# * pi) 
LOOP 
DO WHILE RAofStar# < O 
RAofStar# = RAofStar# + (2# * pi) 


LOOP 
poo IF 
DECFirstColon% = INSTR(DECS(currentStar%), ":") 
DECSecondColon% = INSTR((DECFirstColon% + i 
Dees(currentStarz), ":") 


IF DECFirstColon% = 0 OR DECSecondColon% = O THEN 
RADecError% = TRUE 


Bbokr 
DE COLDS aur # = VAL UCLEEE US (DECS (CeumrencStars) ; 
CoEeCrirstColons = Hee) + VAL(MIDS$(DECS$(currentStar%) , 
M@eeecFirstcolons + 1), 2 1) i 6 0 + 
VAL(RIGHTS (DECS (currentStar%), (LEN (RAS (currentStar%) ) - 
DECFirstColon%))) / 3600 
DECofStar# = DECofStar# * (2# * pi) / 360#'Degrees 


converted to radians 
BORVHEEE DE@otstary > (pil 7 27) 
DECofStar# = pi - DECofStar# 


LOOP 

DO WHILE DECofStar# < -(pi / 2#) 
DECofStar# = -(pi + DECofStar#) 

LOOP 


END IF 

'Account for precession between the epoch chosen by the user 
and J2000.0 

! 


'Determine the Julian date at last midnight 
year! = VAL(EPOCHS$(currentStar?%) ) 
IF year! = 0 THEN 
EOCenEaror. = TRUE 
ESE 
WPOtEpecn! = (47123! + year!) #2365.25 
] 


'Using the approximate formulae for precession in the 
Astronomical Almanac: 
‘Account for precession FROM the star's Known epoch TO 


J2000.0 

T1# = (CDBL(JDofEpoch!) - 2451545#) / 365254 

M1# = 1.2812323# * T1# + .0003879# * Tl# * 2 + 
-0000101# * T1# * 3' degrees 

i= se oio7 SS OOUOC0OOus « Tig=— .0001185e# * T1# * 2 
meeeoOOOll64 * T1# * 3! degrees 

t 

alpha# = RAofStar# * 360# / (2# * pi1)' RA in degrees 


L235 


delta# = DECofStar# * 3607 7 wi2q * Ee Dec in 
degrees 

alphaInRad# 

deltaInRad# 


RAofStar# 
DECofStar# 


'M suffix refers to mean epoch 

alphaM# = alpha# - .5# * (M1# - N1# * SIN(alphalInRad#) 
* TAN (deltaInRad#) ) 

alphaMinRad# = alphaM# * (2# * pi) / 360# 

deltaM# = delta# - .5# * N1l# * COS(alphaMinRad#) 

deltaMinRad# = deltaM# * (2# * pi) / 360# 

1 

'O suffix refers to J2000.0 epoch 

alphaO# = alpha# - M1l# - N1l# * SIN(alphaMinRad#) * 
TAN (deltaMinRad#) 

alphaOInRad# = alphaO# * (2# * pi) / 360# 

delta0O# = delta# - N1l# * COS(alphaMinRad#) 

deltaOInRad# = deltaO# * (2# * pi) / 3604 


‘Account for precession FROM J2000.0 TO the current 


date 

T2# = (CDBL(JulianDate#) - 2451545#) / 365254 

M2# = 1.2812323% * T2# + .0003879% * T24 ue 
.0000101# * T2# * 3' degrees 

N2¢ = .5567530000000001# * T2# - .0001185% * T27e 
= Z000011Gs * = i2 ws a degrees 

1 


'M suffix refers to mean epoch 

alphaM# = alpha0# = oF * (M2? - N2# * 
SIN(alphaOInRad#) * TAN(delta0InRad#) ) 

alphaMinRad# = alphaM# * (2# * pi) / 360# 

deltaM# = deltaO# - .5# * N2# * COS(alphaMinRad#) 

deltaMinRad# = deltaM# * (2# * pi) / 360# 

1 


'The new alpha and delta are the CORRECTED right 
ascension and 


‘declination IN DEGREES 
alpha# = alphaO# - M2# - N2# * SIN(alphaMinRad#) * 
TAN (deltaMinRad# ) 
delta# = deltaO# - N2# * COS(alphaMinRad#) 
RAofStar# = alpha# * (2# * pi) / 360# 
DO WHILE RAofStar# >= (2# * pi) 
RAofStar# = RAofStar# - (2# * pi) 
LOOP 
DO WHILE RAofStar# < 0 
RAofStar# = RAofStar# + (2# * pi) 
LOOP 
1 
DECofStar# = delta# * (2# * pi) / 360# 
DO WHILE DECofStar# > (pi / 2#) 
DECOLStar?; = Pl — De COrorary 
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LOOP 
DO WHILE DECofStar# < -(pi / 2#) 
DECofStar# = -(pi + DECofStar#) 
LOOP 
t 
BND: IF 
BEND SUB 


PKKKKKKKEKKKEKKKEKKEKKEKEKEKEKEKKEKEKKKKKKEKKKKKKKKEK 


x Name: Refraction x 
tak Type: Subroutine kk 
“x x Module: MIRACTRL. BAS x 
tex Language: Quickbasic 4.50 xx 


VHKKKKKEKEKKKKKEKKKKKKKKKKKEKKKKKKKKKKKKKKKEKE 
1 


'Calculates the refraction corrections to Gelestial 
coordinates for a given 

'day and time. Uses data from the Atmospheric Subroutine. 

i] 

'Example of use: Refraction Pressure!, Temperature!, 
Wavelength!, GMSTOhrSEC#, HourAngle#, GHAaries#, RAofStar#, 
DECofStar#, RARefractCor#, DecRefractCor#, RAofMIRA# 

t 

SUB Refraction (Pressure!, Temperature!, Wavelength!, 
GMSTOhrSEC#, HourAngle#, GHAaries#, RAofStar#, DECofStar#, 
RARefractCor#, DecRefractCor#, RAOfMIRA#) 


hhymedalize variables 


MIRAlatitude# = .63908139# 'radians, (36 deg 37 min 
NOrech) 

MIRAlongitude# = 2.12639281# iead tans. i( le deg. SO min 
West) 

sec% = VAL(RIGHTS(TIMES, 2)) 


) 
min% = VAL(MIDS(TIMES, 4, 2)) 
hour% = VAL(LEFTS(TIMES, 2)) 
t 
'Calculate Greenwhich Mean Sidereal Time 
'Determine the day fraction that has passed since midnight 
fee oe CONDS 


Seas = secs 
hours& = hour% 
mansk = mins 


dayFraction# = CDBL((hours& * 3600) + (mins& * 60) + sec&) 
q 


'Determine the Greenwhich Mean Sidereal Time at the current 
'time in SECONDS 

GMSTinSEC# = GMSTOhrSEC# + dayFraction# 
i] 
'Calculate Greenwhich Hour Angle of Aries 

GHAaries# = GMSTinSEC# * 2# * pl / 86164.09053099999# 'Seconds 
of solar time converted to radians 


2 


DO WHILE GHAaries# > (2# * pi) 
GHAaries# = GHAaries# - (2# * pi) 
LOOP 
DO WHILE GHAaries# < 0O# 
GHAaries# = GHAaries# + (2# * pi) 
LOOP 
1 
'Calculate the current star's hour angle for this moment 
RAOfMIRA# = GHAaries# - MIRAlongitude# 
i] 
HourAngle# = RAofMIRA# - RAofStar# iy. rad hans 
DO WHILE HourAngle# > (2# * pl) 
HourAngle# = HourAngle# - (2# * pi) 
LOOP 
DO WHILE HourAngle# < 0O# 
HourAngle# = HourAngle# + (2# * pi) 
LOOP 
1 
'Calculate the refraction constant in seconds 
k# = 21.3 * Pressure! * (1 + (.0057 / (Wavelength! * 2))) / 
(273 + Temperature! ) 
'Note: Temperature in degrees C, Pressure in mmHg, Wavelength 
in microns 
K# in seconds of arc 
i] 
‘Calculate the refraction corrections in seconds 
RARefractCor# = k# * SIN(HourAngle#) / ((COS(DECofStar#) * 2) 
* (TAN(DECofStar#) * TAN(MIRAlatitude#) + COS(HourAngle#) )) 
DecRefractCor# = k# * (TAN(DECofStar#) * COS(HourAngle#) - 
TAN (MIRAlatitude#)) / (TAN(DECofStar#) * TAN(MIRAlatitude#) + 
COS (HourAng1le#) ) 
RARefractCor# = RARefractCor# * (2# * pi) / 3600#' Convert 
seconds of arc to radians 
DecRefractCor# = DecRefractCor# * (2# * pi) / 3600#' Convert 
seconds of arc to radians 


END SUE 

VK KKK KKK KKK AK KA K KKK KEK KK KEK KEK KKK AK KKK 
Bees Name: SetCommands kx 
kx Type: Subroutine * x 
'kx* Module: MIRACTRL. BAS xx 
(kx Language: Quickbasic 4.50 * 


VRKREKKEKKKKKKKKKKKKKKKKEKKKKKEKKEKKEKKKKKKKEKEK 

1 

‘Controls motor commands to set the Telescope following a slew 
to a new 


'target star. Uses the stepping motors in high speed mode. 
1 


120 


'Example of use: SetCommands Simulation%, TeleRA#, TeleDec#, 
RAofStar#, DECofStar# 


SUB SetCommands (Simulation%, TeleRA#, TeleDEC#, RAofStar#, 
DECOLS tar #F } 
i] 


'Set as appropriate in Right Ascension 
DO 
RAGiff# = TeleRA# - RAofStar# 
'Update motor display box 
GCECKUKED? ~ BLUE, BLUE 
POCATE 25. 20, 0 
PRINT CHRS$ (222) 
‘POKE (address of RA Set controller), X 
PKKKKKEKKKEKKKKKKKKKKKKEKEK 


IF Simulation% = TRUE THEN 


IF RAGiff# > O THEN TeleRA# = TeleRA# - .00001 
IF RAdiff# < 0 THEN TeleRA# = TeleRA# + .00001 
ELSE 
‘Read new TeleRA# from encoders 
195 ee © 2 


'Check limit switch 

'PEEK (address of limit switch) 

PK KKKKKKEKKKKKKKKKEKKKKKEKKKKKKKKEKK 

'k* Change RAdiff# < .0003 to ** 

'kkx acceptable limit * 
BOOP UNTIL RAdiIff# < .0003 


PHRKKKKKKKKKKKKK KK KK KKKKKKKKKKKKK 
i) 


'Set aS appropriate in Declination 
DO 
BDEGAITEFe— TeleDEC? — DECofStar# 
'Update motor display box 
COLOR RED, BLUE, BLUE 
BCCAITE. 3, 70,0 
PRINT CHRS$ (222) 
'POKE (address of DEC Set controller), X 
PKKKKKKEKKKEKKKKKKKKKKKKEKEK 
IF Simulation% = TRUE THEN 
LE wwpEeCadtiti? = 70) THEN TelebpDECr 


TeleDEC# - 
OOOO | 


TF DECdiff# < O THEN TeleDEC# 


li 


TeleDEC# + 
- 00001 


ELSE 


‘Read new TeleDEC# from encoders 
END IF 


] 
VWeheck limite Switch 


'PEEK (address of limit switch) 
VHKKKKKEKKKKKKKKKKK KK KKK KKK KKKKKKEK 


ee 


'x* Change DECdI fis =< 7000s stom. 
'x* acceptable limit mts 


LOOP UNDPTLE DEGGIE t?Y soo Ge 
ESE SC CCCCCCCCCCLCCCSCCSSCCCSO. F 


END SUB 

PK KKK KKK KK KK KK KKK KKK KKK KKK KKK KKEKKKKKKKEKK 
1k x Name: ShowTime * x 
tk Type: Subroutine ne 
kx Module: MIRACTRL. BAS ** 
kx Language: Quickbasic 4.50 x 


VK KK KKK KKK KKK KEK KKK KEKE KKK KKK KKEKKEKKKKKKKEK 

s 

'Allows entry of star position data aS an initial user input 
i] 

'Example of use: ShowTime min%, hour%, day%, month%, year%, 
DayLightSavings%, SiderealTime$ 

‘Parameters: 

'Variables: SiderealTime$ A string representing the 
local 

apparent sidereal time 
DayLightSavings% A variable showing the 
State of daylight 


Module Level 
DECLARATIONS: DECLARE SUB ShowTime (min%, hour%, day%, 
month%, year%, DayLightSavings%, SiderealTime$) 
t 
SUB ShowTime Ouse 5 lovey eee. - days, month%, years, 
DayLightSavings%, SiderealTime$) 
'Display local time in accordance with daylight 
Savings time. 
IF DayLightSavings% = FALSE THEN 
LocalTimesS = STRS(VAL(LEFTS(TIMES, 2)) - 5) + 
RIGHTS (TIMES, 6) 
je Ley 


Savings time - on or off 
' min% The current minute of time 
hour% The current hour of time 
day% THe esr ent mela, 

: month% The current month 

years The current year 

f 

f 

| 


LocalTimeS = STRS(VAL(LEFTS(TIMES, 2)) —3 
RIGHTS(TIMES, 6) 
END. TF 


{ 
IF VAL(LEFTS$(LocalTime$, 2)) < 1 THEN 
LocalTime$ = STRS(VAL(LEFTS$(LocalTimes$, 2)) + 
24) + RIGHTS( Local Times aaa 
END IF 


is 


Mipseeclonkecs — INSTR(SiderealTimesS, ":") 


LastColonLoc% = INSTR((FirstColonLoc% + i) 
SaderealTimes, ":") 

LengthSec% = 2 

BeginMin% = FirstColonLoc% + 1 

t 


SiderealSec$ = LTRIMS(STRS$(VAL(RIGHTS (SiderealTimeS, 
LengthSec%)) + 1)) 
SiderealMin$ = MID$(SiderealTime$, BeginMin?, 
LastColonLoc% - FirstColonLoc% - 1) 
SiderealHr$ = LEFTS$(SiderealTime$, (FirstColonLoc% - 
) 
t 
IF LEN(SiderealSec$) = 1 THEN SiderealSec$ = "0" 4+ 
SiderealSec$ 
IF VAL(SiderealSec$) > 59 THEN 
SiderealMin$ = STRS$(VAL(SiderealMin$) + 1) 
SiderealSec$ = "00" 
END IF 
IF VAL(SiderealMin$) > 59 THEN 
SiderealHr$ = STRS(VAL(SiderealHr$) + 1) 
SiderealMins$ = "oOo" 
PND TE 
IF VAL(SiderealHr$) > 23 THEN SiderealHr$ = "oOo" 
SiderealTime$ = LTRIMS$ (SiderealHr$) + ce + 
LTRIM$S(SiderealMins$) + ":" + LTRIMS(SiderealSec$) 
COLOR WHITE, BLUE, BLUE 
LOCATE 15, 58, 0 
PRINT LTRIMS (LocalTime$) 
BOCATE 15,556, 0 
PRINT LTRIMS(SiderealTime$) ; SPACE$(1) 
LOCATE 21, 58, 0 
PRINT TIMES 
END SUB 


Vaekk kk kkk eK KKKKKKK KK KKK KK KKK KKKEKE 


kk Name: SiderealTimeCalc ** 
kx Type: Subroutine xx 
ei Module: MIRACTRL. BAS kk 
tk x Language: Oimehe so aiomea 4.5.0 * x 


Dake kkk keke KKK KKK KKK KKK KKK KKK KKKE 
t 


‘Calculate Local Apparent Sidereal Time 
t 


'Example of use: SiderealTimeCalc SiderealTime$, min%, hour?, 
day%, month%, year%, EQofEQUIN!, GMSTOhrSEC#, JulianDate#, 
dayFraction#, DaysPassed% 


‘Parameters: (none) 
"Variables: SiderealTime$ A string representing the 
ieca ] 


corrections to 
t 

dates (due 

t 

nutation 

| 


Since a known 
t 


above in string 


the value of 
t 


not the 
t 


t 
Or NOG noon 
t 


the current 
t 


seconds 
t 


calculation of 
t 


ato N@eiar 
t 


Astronomical 
t 
t 
at oO Nour som 
t 
t 


current day that 


at the cunmrrent 


length when a 
t 


reset 


apparent sidereal time 


EQofEQ! () The Equation of Equinge: 
sSidereal time for successive 
to corrections for the Earths 
ming The current minute of time 
houses The current hour of time 
day% The current day 
month% The current month 
years iie Selieine Miemy Cae 
yearDiff! The number of years passed 


leap year (1992) divided by 4 
yearDiffStrg$ The decimal portion of the 


form. yearDiff! then becomes 


this string 


LeapYears$ A string indicating whether or 
current year 1S a leap year 
noonPassed$ A string representing whether 
has passed on the current day 
daysPassed% The number of days passed in 
year 
JulianDate# The current Julian date in 
Tu# A variable used in the 
Greenwhich Mean Sidereal Time 
following the convention in the 
Almanac 
GMSTOhrSEC! Greenwhich Mean Sidereal Time 
age Qahiciceicne tclsky 
GdayFraction! The decimal fraction of the 
has passed by the current time 
GMSTInNSsEC! Greenwhich Mean Sidereal Time 
time in seconds 
OKtoQUITS A string that becomes of some 


loop is to be exited, then is 
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'Module Level 


; DECLARATIONS: DE CEA E SUB SiderealTimeCalc 
(SiderealTime$, min%, hour%, day%, month%, year%, EQofEQUIN!, 
GMSTOhrSEC#H, GMSTOhrSEC#H, JulianDate#, dayFraction#, 


DaysPassed%) 
i] 


SUB SiderealTimeCalc (SiderealTime$, min%, hour%, day%, 
Months , years, EQofEQUIN!, GMSTOnrSEC#, JulianDate#, 
dayFraction#, DaysPassed%) 

'Determine if the current year 1S a leap year 
yearDiff! = CSNG(year% - 1992) / 4! 
yearDiffStrg$ = REGHTS(S TRS (yearDifft!:)., 

(LEN(yearDiff!) - INSTR(STRS$(yearDiff%), "."))) 

yearDiff! = VAL(yearDiffStrg$) 
IF yearDiff! > O THEN 
LeapYearS = "FALSE" 
ese 
LeapYear$S = "TRUE" 
END aE 
i] 
'Determine if noon has passed on the current day 
ite hours. = 2 TREN 
noonPassedS = "FALSE" 
BLSE 


noonPassed$ eERUEY 
jC Be 
'Determine days that have passed in the current year 


SELECT CASE month% 


CASE: <1 
DaysPassed% = day% - 1 
CASE. 2 
DaysPassed% = day% + 30 
CASE 3 
IF LeapYearS = "TRUE" THEN 
DaysPassed% = day% + 59 
ELSE 
DaysPassed% = day% + 58 
END er 
CASE 4 
IF LeapYearS = "TRUE" THEN 
DaysPassed% = day% + 90 
ELSE 
DaysPassed% = day% + 89 
EN DoE 
GASES 
IF LeapYear$ = "TRUE" THEN 
DaysPassed% = day% + 120 
jdalrey o 
DaysPassed% = day% + 119 
1d) B Us 3 


CASE 6 
IF heapYears = “TRUE 
DaysPassed% = 
1h poy 6 
DaysPassed% = 
END i: 
CASE, 7 
IF LeapYearS = "TRUE" 
DaysPassed% = 
Pio ee 
DaysPassed% = 
END IF 
CASE 8 
IF LeapYear$ = "TRUE" 
DaysPassed% = 
JES I 
DaysPassed% = 
END VEE 
CASE 9 
IF LeapYear$ = "TRUE" 
DaysPassed% = 
EGSE 
DaysPassed% = 
END 
CASE 10 
IF LeapYear$S = "TRUE" 
DaysPassed% = 
Pe E 
DaysPassed% = 
JHINID) 2018 
CASEML 
IF LeapYearS = "TRUE" 
DaysPassed% = 
ELSE 
DaysPassed% = 
ENE Sir 
CASE 2 
IF LeapYears$S = "TRUE" 
DaysPassed% = 
ELSE 
DaysPassed% = 
END IF 
CASE ELSE 


END SELECT 


2 


THe 
day% 


day% 


THEN 


° 


day% 
day% 


THEN 


° 


day% 
day% 


THEN 


° 


day% 
day% 


THEN 


° 


day% 


day% 


dels) 
day% 


days 


TaN 
day% 


° 


days 


'Determine the Julian date at last midnight 
PKK KK KKK KKK KKK KKK KKK KK KK KKK KKK KKK KEKE KKK KKK KEKE KKKEKKKKKKEKSEK 
Change to be good for all years, not just 1992 ** 


JulianDate# = CDBL(2448621.5# + DaysPassed% + 1) 
PKKKKKKKKK KKK KKK KKK KKK KKK KKK KKKKKKKE KK KKK KEKEKEKKKEKKKKKKEKSE 


aa 


Loe 


Toe 


180 


EAN 


2 dy 


243 


242 


273 


Zoe 


304 


308 


334 


333 


WeEC 


'Determine the Greenwhich Mean Sidereal Time at O hour 


rin SECONDS 
Tu# = (JulianDate# - 2451545#) / 36525! 
GMSTOnNrSEC# = 24110.54841# + (8640184.812866# * Tu#) 


meeeesslO4 * (Tu) ~*~ 2) ~- (6.2 * (10 ~*~ (-6)) * (Tu) “* 3) 


'Determine the day fraction that has passed since 


midnight 
'in SECONDS 
hounrsé = hours 
mins& = min% 


current 


SECONDS 


dayFraction# = CDBL((hours& * 3600) + (mins& * 60)) 
t 
'Determine the Greenwhich Mean Sidereal Time at the 


'Cime in SECONDS 
GMSTinSEC# = GMSTOhrSEC# + dayFraction# 


'Calculate Greenwhich Apparent Sidereal Time in 


GASTInSEC#] — GMsSTinSEC?, + CDBL(EQofEO!) 
'Add or subtract multiples of 24 hours as necessary 
DO 
Le GASPinsehe7 = 0: THEN 
GASTinSEC# = GASTinSEC# + (24! * 3600!) 
ELSEIF GASTinSEC# > 86400 THEN 
GASTinSEC# = GASTinSEC# - (24 * 3600) 
END: LF 
IF GASTinSEC# >= 0O AND GASTinSEC# <= 86400 


THEN OKtoQUIT% = TRUE 


LOOP UNTIL OKtoQUIT% = TRUE 
SkKEOOUITs = TRUE 
t 


GASThr& = INT(GASTinSEC# / 3600) 
GASTmin%’ = INT((GASTinSEC# - (GASThr& * 3600)) / 60) 


GASTsec% = INT(GASTinSEC# - (GASThr& * 3600) - 
(GASTmin% * 60)) 

SiderealTime$ = LTRIMS(STRS(GASThr&)) + ":" + 
LTRIMS(STRS(GASTmin%)) + ":" + LTRIMS(STRS(GASTsec$%) ) 


END SUB 


PHKKKKKKKKKEKKKKKEKKKEKKEKKEKKKKEKKKEKKKKKKKKKEKEK 


tax 
tex 
tk 
tk & 


Name: SlewCommands xk 
Type: Subroutine * x 
Module: MIRACTRL. BAS * x 
Language: Quickbasic 4.50 * x 


PAKKKKEKKEKKKEKKEKKKKKKKKKKKKKKKKKKKKKKKKKKKK 


‘Controls motor commands to slew the Telescope to a new 


‘target 


Star. Uses the slew motors. 
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‘Example of use: SlewCommands (Simulation%, TeleRA#, 
TeleDec#, RAGLStarF, DECOrseans 
t 


SUB SlewCommands (Simulation%, TeleRA#, TeleDEC#, RAofStar#, 
DECofStar# } 
t 


'Slew as appropriate in Right Ascension 
DO 
RAdiff# = TeleRA# - RAofStar# 
DO WHILE RAdiff# >= (2# * pi) 
RAdGiff# = RAdiff# - (2# * pi) 
LOOP 
DO WHILE RAdiff# <= -(2#4 * pi) 
RAGiff# = RAdiff# + (2# * pi) 
LOOP 
‘Update motor display box 
COLOR RED, BLUE, BLUE 
hOGA NRIs 3, ssiKOnmne 
PRINT CHR$(222) 
‘Turn RA slew motor ON 
'POKE (address of RA slew controller), X 
VaekKeKK KK KKK KKKKKKKKKKEK 
IF Simulation% = TRUE THEN 
IF RAdiff# > O THEN TeleRA# 
IF RAdGiff# < 0 THEN TeleRA# 


TeleRA# - .0001 
TeleRA# + .0001 


EBSE 
'Read new TeleRA# from encoders 
END aie 
i] 
'Check limit switch 
'PEEK (address of limit switch) 
Va KKKKKKKKKKKKR KKK KKKK KKK KKKKKKKEK 
'kx*x Change RAdiff# < .0003 to ** 
'xx acceptable limit aia 
LOOP UNTIL RAdiff# < .0003 
VHKKKKKKKKKKKKKKKK KK KK KKK KKK KKKKSK 
‘Update motor display box 
COLOR REDe BLUE, BLUE 
LOCATE 3 alee © 
PRINT SPACES (1) 
q 
'Slew as appropriate in Declination 
DO 
DECdiff# = TeleDEC# - DECofStar# 
‘Update motor display box 
COLOR RED, BLUE, EEUE 
LOCATE 3, 60, O 
PRINT CHR$(222) 
‘Turn DEC slew motor ON 
'POKE (address of DEC slew controller), X 
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VK KKK KKK KK KKK KKK KKK KKKKK 
TF Simulations = TRUE THEN 
hE DECaatf? > O FHEN TeleDEC# = TeleDEC# - 
ZOO 1 


IF DECdiff# < O THEN TeleDEC# TeleDEC# + 
SOOO] 

16 Grea 

‘Read new TeleDEC# from encoders 

END LE 

1 

Sehecko Limit Swheen 

'PEEK (address of limit switch) 

Va KkKk KKK KKK KK KKK KK KKK KKK KK KKK KKK 

'x* Change DECdiff# < .0003 to ** 

'x*x acceptable limit * x 
ifeoP UNTIIM DECdiIff£f# < .0003 

Vaekkkekke kk Keke RK KK KK KKK KK KK KKK KKK 

‘Update motor display box 

CGOECRSRED; SELUE, BEUE 

LOCA 3, 60, 0 

PRINT SPACES(1) 


END SUB 


CK KKK KK KKK KEKKKEKKKKKKKEKKKKKKKKKKKEKKKKKKK 


2a tas Name: Secansvacabncny xx 
kx Type: Subroutine x 
kx Module: MIRACTRL. BAS xx 
kx Language: Quickbasic 4.50 x 


PKK KK KEKEKKKKKKKKKEKKKKKKKEKKKEKKKKKKKKKKKKEK 
t 


‘Allows entry of star position data aS an initial user input 
1 

Peeeomole of use: StarDataEntry Stars(), RAS(), DECS(), 
Epoch$(), NewFile%, wl, wiText$(), wlTitle$S, wlPrompts 


‘Parameters: (none) 

‘Variables: Seae > () A 5-element matrix of desired 
Stars for 

viewing 

: RAS () A 5-element matrix of right 
ascensions 

: for the stars listed in StarS$() 
DECS () A 5-element matrix of declinations 
' for the stars listed in StarS$() 

; Epechs() A 5-element matrix of epochs for 
the 

! RA's and DEC's of the stars listed 
mimo tars ( ) 


g%, h%, 1% Looping indices 
2 A flag activated when the 
Backspace key is hit 


13 


1 
Tab wp cr sbDown 
{ 


ove 


A flag activated when the Shift 


Arrows are hit 
row%, cole The current row and column markers 


for the 

CULESO: 

Stops A string representing whether or 
not the 

character entry loop is to be 
exited 


kee% A unigue code returned by the 
KeyCode% or 
KeyCodeNoWait% Functions for each 
key hit 
Character$ A string read from the keyboard 
buffer 
Current$ A string representing the entry of 
the current 

field. This variable will become 
one of the 

following: StarS(), RAG@e 
DECS(), EPOCHS() 

1 


"Module Level 

DECLARATIONS: DECLARE SUB StarDataEntry (Star$(), RAS(), 
DEC Sty 
Epoch$(), NewFile%, wl, 
wlText$(), wiTitle$, wiPrompts) 

{ 


SUB StarDataEntry (Star$(), RA$(), DECS$(), EPOCHS(), NewFile3, 
wl AS WindowsType, wlText$(), wilTitle$, wilPrompt$) 
t 
'Initialize variables 
IF NewFile% = TRUE THEN 
FOR he = 1 TO 100 
Stars({(hsy = "" 


RAS (h%) — FON 
DECS$(h%) = "" 
EPOCHS(h%) = "" 
NEXT h% 

ENDS LE 

current$ = Star$(1) 

J = 0 

L% = 0 

1% = 1 

rows = 17 

col% = 9 

InsertOn% = FALSE 

FirstRows = 17 


LaStROws  =]—2. 
LastCol. = 62 
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FileLine%’ = 1 
TopFileLineShown% = 1 
ChangeCol% = FALSE 
ChangeFileLine% = FALSE 
HomeOrEnd% = FALSE 
Stops% = FALSE 
QuitEntry% = FALSE 
t 
“praw bottom edge of the status box 
MOCATE 22, 5, 0 
PRINT CHR$(200); STRINGS$(68, 205); CHR$(188); 
t 
'Complete data fields by reading in one key at a time, ESC 
exits 
DO 
IF TopFileLineShown% < 1 THEN TopFileLineShown% = 1 
IF TopFileLineShown% > 96 THEN TopFileLineShown% = 96 
IF FileLine%’ < 1 THEN FileLine% = 1 
IF FileLine% > 100 THEN FileLine% = 100 


IF row% < FirstRow% THEN row% = FirstRow% 
IF row% > LastRow% THEN row% = LastRow% 


‘Update the row numbers and fields 
FOR g% = 1 TO 5 
COLOR “WHET E BLUE, BLUE 
LOCATE 16 + g%, 6, O 
PRINT TopFileLineShown% + g% - 1 
COLOR WHITE, “BELACK, BLUE 
LOCATE 16 + g%, 10, 0 
PRINT Star$(TopFileLineShown% + g% - 1) + 
SPACES(19 - LEN(Star$(TopFileLineShown% + g% - 1))) 
LOCATE 16 + g%, 33, 0 
PRINT RAS(TopFileLineShown%’ + g% - 1) + 
SPACE$(10 - LEN(RAS$(TopFileLineShown% + g% - 1))) 
LOCATE 16 + g%, 48, 0 
PRINT DECS$(TopFileLineShown% + g% - 1) + 
SPACES(10 - LEN(DECS(TopFileLineShown% + g% - 1))) 
LOCATE 16 + g%, 63, O 
PRINT EPOCHS (TopFileLineShown% + 
SPACES(10 - LEN(EPOCHS (TopFileLineShown% + g% - 1 
NEXT g% 
COLOR WHITE, BEAGCK, “BLUE 
"Read the next character and operate on it 
DO 
ROGATE “YOWcem cols — 22, 21 
IF 1% = O THEN 


1% = 1 
J% = 0 
END IF 
IF 3% = 1 THEN 
PRInT 2." 


IF LEN(current$) > 1 THEN 


inSe7 


currentS = LEFTS (cunmerims 
(GEM CGUGrenes) = lo) 


ELSE 

currents = "™" 
END EE 
1% = 0 


END ak 
LOGATE rows, Gollcemeie, 1 
kee% = KeyCode% (Character$) 
SELECT CASE kee% 
CASE ENTER 
IF row% = LastRow% AND col% 


LastCol?.32HEN 
rows = LOW. 
IF row% = LastRow% 


+ 


1 THEN 


TopFileLineShown% = TopFileLineShown% + 1 
FileLine% = 
FileLineS + 1 
BND LE 
ELSE 
'MoOsacti on 
END TE 
IF col% = 9 THEN 


ELSE 


ND Le 
Stops% = TRUE 
CASE TABHIT 


JEQIEYS dos 
CASE ESCAPE 


windows wl, wlText$(), 
wWITitles ,wWEPronprs 


IF wl.returnCode = 121 OR 
wlsreturnCode. = 89 THEN Abert. s—— Rue 


IF Abort% = TRUE THEN 


PH KK KI HKIKKKHKKKKEKKEKKKKKHEKHaKKKEKKKEKEK 


'kk System Shutdown 
Procedures = 


ek go here!! 
x * 


Taekkkkkk kkk kkkkk kkk kk kkkkkkke kk kK KKK 


WindowsPop 
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CLS 


Se Geel ul 
Jap lis) le 
Windows Pop 
SNE DY 2) a 
CASE BACKSPACE 
eS 
i) oie 


GASE SHIFT BAB 
tieeols = 32 THEN. cols = 9 
Lheee ls 47 SHEN cols 3.2 
tieeeowic, — 62 THEN cols = 947 
Stops% = TRUE 

CASE UPARROW 
IF row% = FirstRow% THEN 

TopFileLineShown%$ = 


TopFileLineShown% - 1 
BESE 
rows = row% - 1 
END iF 
ChangeFileLine%’ = MINUS 
Stops% = TRUE 
CASE DOWNARROW 
IF row% = LastRow% THEN 
TopFileLineShown% = 
TopFileLineShown% + 1 
BESE 
rows = rows + 1 
PhD oie 
ChangeFileLine% = PLUS 
Stops% = TRUE 
CASE PGDOWN 
IF TopFileLineShown% > 95 THEN 
ChangeFileLine% = PLUS 
* (96 - TopFileLineShown$) 
TopFileLineShown% = 96 


Blok 
TopFileLineShown% = 
TopFileLineShown% + 5 
ChangeFileLine%’ = PLUS 
me 
ae EF 
Stops% = TRUE 
CASE PGUP 
IF TopFileLineShown% < 6 THEN 
ChangeFileLine% = 
MINUS * (TopFileLineShown% - 1) 
TopFileLineShown% = 1 
Papo. 


TopFileLineShown% = 
TopFileLineShown% - 5 


1339 


ChangeFileLine% 


MINUS * 5 

END EE 
Stops% = TRUE 

GASE HOME 
TopFileLineShown% = 1 
HomeOrEnd% = 1 
Stops% = TRUE 

CASE ENDFILE 
TopFileLineShown% = 96 


HomeOrEnd% = 100 
Stops% = TRUE 
CASE LEFTARRKOW 
1% = 1% - 1 
CASE RIGHTARROW 
IF LEN(currents$) > 1% THEN 
i$ = 1% + 1 
END IF 
CASE, DEE EE 
TF 1% > i AND 1s < 
LEN (CuRBem@es) + 10 tREN 
GCUrrents = 
LEFT$(current$, 1% - 1) + RIGHTS(current$, LEN(current$) - 1%) 
ELSEIF 1% = 1 AND i2e 
LEN(currents) =a] THEN 


Curren = 
RIGHTS (current$, LEN(current$) - i%) 
ELSEIF i% = 1 THEN 
Currents <1) 
Bboe 


Currents 
LEPTS (eurrents 1.5 — i) 
END IF 
LOCATE DrOw2 Co. 5 lee 
IF col% = 9 THEN 
PRINT Currents, 
SPACES(19 - LEN(currents) ) 
ELSE 
PRINT currents; 
SPACES(10 = LEN(currents) ) 
END IF 
CASE INSERT 
IF InsertOn% = TRUE THEN 
Inserton% = FALSE 


LOCATE , , , 7 
ELSE 

InsertOn% = TRUE 

LOGATE *, (3) 3) eae 
BNDe tis 


CASE 17408 'F10 
QuitEntry% = TRUE 
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Stops% = TRUE 
GASE 97 TO 122 

AcceptChar% = TRUE 
GASEV6o TO, 90 

AcceptChar% = TRUE 
CASE 48 TO 57 

AcceptChar% = TRUE 
CASE 58 TO 59 

AcceptChar% = TRUE 


CASE 34 

AcceptChar% = TRUE 
CASE 39 

AcceptChar% = TRUE 
CASE 44 

AcceptChar% = TRUE 
CASE 46 

AcceptChar% = TRUE 
CASE 32 

AcceptChar% = TRUE 
CASE ELSE 


PKHKKKKKKKKEKKKKKKKKKEEKKEKEKKKKKKKKKKKKRKKKKEKE 


tke 
zxk* 
tak Error Handling 
Needed x 
tek 
xk* 


PK KKKKKEKKKKEKEKEKEKEKKKKEKEKEKKKEKEKKKKKKKKKKKKKEEEE 
END SELECT 
IF col% = 9 AND i% = 20 THEN 
Stopes = TRUE 
ChangeCol% = TRUE 
ELSEIF col% > 9 AND 1% = 11 THEN 
Stops% = TRUE 
ChangeCol% = TRUE 
END IF 
IF AcceptChar% = TRUE THEN Stops% = TRUE 
LOOP UNTIL Stops% = TRUE 
Stops% = FALSE 
IF AcceptChar% = TRUE THEN 
IF InsertOn% = TRUE THEN 
ie ae > 1 AND@is <SBEN(currents) + 1 
THEN 
currents = LEFTS(currents, 1% 
- 1) + CharacterS + RIGHTS(currentS, LEN(current$) - 1% + 1) 
ELSEIF i% = 1 AND 1% < LEN(currents) 


tae LHEN 
currents = Characters + 
RIGHTS (current$, LEN(currents$) - 1% + 1) 
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ELSEIF 1% = 1 THEN 
currentS = Characters 
FRESE 
currents = LEFTS (currentS 73a 
- 1) + Characters 
END IF 
ELSE 
IF 1% > 1 AND i% < LEN(current$) THEN 
currentsS = LEFTS$S(currents, i% 
- 1) + CharacterS + RIGHTS(currentS, LEN(currents) —- 12) 
ELSEIF 1% = 1 AND i% < LEN(current$) 


THEN 
currents = Characters + 
RIGHTS (current$, LEN(current$) - 13%) 


ELSEIF 1% = 1 THEN 
currentS = Characters 
EMSe 
currents = LEFTS(currents, 1% 
- 1) + Characters 
END IF 
END IF 
LOCATE rowzt-eColic ati) 
IF col% = 9 THEN 


PRINT currents; SPACES (19 = 
LEN (currents$) ) 
a Sy 
PRINT Currents; SPACES (10 = 
LEN (Current$) ) 
PND Le 
1% = 1% + 1 
AcceptChar% = FALSE 
PND 
IF FileLine% < 1 THEN FileLine% = 1 
IF FileLine% > 100 THEN FileLine% = 100 
SELECT. GaSe. cols 
CASE 9 
Star$(FileLine%) = currents 
CASE [32 
RA$(FileLine%) = currents 
CASE 47 
DECS (FileLine%) = currents 
CASE 62 
EPOCH$(FileLine%’) = currents 
CASE ELSE 


PKKKEKKKKKEKEKKEKEKEKEKKKEKEKEKEKEKEKEKEKEKKEKKEKKEKKKEKE 
kx 


tx Error Hanging 
Needed x 
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t*k* 
*k* 


DKK KKK KKK KK RKKKKKEKKAKKKKKEKK EKER KKKKKKKEKKS 
END SELECT 


1 

IF ChangeFileLine%’ <> FALSE THEN 
FileLine% = FileLine% + ChangeFileLine% 
IF FileLine% < 1 THEN FileLine% = 1 
IF FileLine% > 100 THEN FileLine% = 100 
ChangeFileLine%’ = FALSE 
SELECT CASE col% 


CASE 9 

current$S = Star$(FileLine?%) 
CASE 32 

currentS = RAS$(FileLine?) 
CASE 47 

currentS = DECS$(FileLine?%) 
CASE 76:2 

currentS = EPOCHS (FileLine?) 
CASE ELSE 


PKKKKEKKKEKKEKEKKEKKEKKEKEKEKEEKEKEKEKKEKKEKKEKKEKEREKEEKER 


fax 
Kx 
tk* Error Handling 
Needed xx 
xx 
kk 


PKKEKKKKKKEKKKEKKEKKEKEKKKEKEKEKEKEKEKEKKEKEKKKEKRKEKKEKEKEER 


END SELECT 
END iF 
| 
IF HomeOrEnd% <> FALSE THEN 
IF HomeOrEnd% = 1 THEN 
row% = FirstRow% 
ELSEIF HomeOrEnd%é = 100 THEN 
rows = LastRow% 
END IF 
FileLine% = HomeOrEnd% 
HomeOrEnd%s = FALSE 
COLOR WHITE, BLUE, BLUE 
LOCATE LastRow%, 9, O 
PRINT SPACES (1) 
SELECR VeASE “cols 


CASE 9 

currentS = Star$(FileLine?%) 
CASE 32 

current$ = RAS(FileLine?%) 
CASE 47 
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GUrrenc> DECS (FileLine%) 


CASE 62 


currents 
CASE ELSE 


EPOCHS (FileLine?%) 


PK KK KKKKKKEKKKKKKKKEKKKKEKKKEKKEKKKKKKEKKKEKEE 


kx 
*kx 
tk x Error Handling 
Needed ** 
tk 
xx 


PKK KKKKKKKKKKKKKKKEKKKKKKKKEKKKKKKKEKKKKEKKEKK 


END SELEGr 
END IF 
| 
IF Changecol% = TRUE THEN 
SELECT CASE col% 


CASE 9 
col% = 32 
current$S = RAS(FileLine%) 
CASE SZ 
col% = 47 
current$ = DECS$(FileLine%) 
CASE 47 
col%é = 62 
current$ = EPOCHS(FileLine%) 
CASE OZ 
col% = 9 
rows = rows + 1 
FileLine% = FileLine% + 1 
IF FileLine%’ < 1 THEN FileLine% = 1 
IF FileLine% > 100 THEN FileLine% = 
100 
currentsS = Star$(FileLine%) 
CASE ELSE 


PKKKEKKKEKKKKEKKEKEKKKEKEKKEKEKKKEKKKEKKKKKKKKKKKKEKK 


tax 
** 
tk x Error Handling 
Needed ** 
xx 
xx 


DK KK KK KKK KK KEKE KK KEK KKK KKKKKEKEKKKKKKKKKKKEK 
END SELECT 


1% = 1 
Changecol% = FALSE 
END IF 
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LOOP UNTIL QuitEntry% = TRUE 


END SUB 

DK HKKKHKKKKKKKEKKEKKKKEKKKKKKKEEKKEKEKEKKEKEKEEE 
Vex Name: StarDisplay kx 
Vex Type: Subroutine ** 
Vex Module: MIRACTRL. BAS ** 
Vex Language: Outekbasice 4.50 ** 


PHKKKEKKEKKKEKEKEKKEKEKKEKEKKEEKKEKEKEKKEKEKEKKEKEKEKKKEEEE 
LJ 


'Creates the display screen graphics for star data entry 
t 


'Example of use: StarDisplay 


‘Parameters: (none) 
'Variables: 1% Looping index 


‘Module Level 
' DECLARATIONS: DECLARE SUB StarDisplay 


SUB StarDisplay 
§ 


'Clear field 
COLOR WHITE, BLUE, BLUE 
FOR r% = 14 TO 21 
LOGATE 36% ,.26,. 0 
PRINT SPACES (68) 
NEXT r% 
1 
'Fill in star data entry fields 
COLOR WHITE, BLUE, BLUE 
M@GATE 14, 31, 0 
PRINT "Please enter stars" 
LOCATE 15, 33, 0 


PRINT "RA DEC" ; 
mock 16, 10, O 
PRINT "Star Name heemins sec: 


LOCATE 16, 48, O 
Biren CHRS(248); ":':t'u: 
BeOGCATE 16, 63, O 
PRINT "Epoch"; 
1 
'Draw fields for star data input 
FOR 1% = 1 TO 5 
COLOR WHITE, BLACK, BLUE 
E@GATE f6.+ 122) 10,70 
PRINT SPACES (19) 
LOCATE NIGe+ 1%, 33,0 
PRINT SPACES (10) 
LOCATE 16 + 1%, 48, O 
PRINT SPACES$(10) 
LOGATE 16 + 1%, 63, 0 
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PRINT SPACES (10) 
NEXT 1% 
LOCATE 23 a5 aro 
COLOR WHITE, BLUE, BLUE 
PRINT "F10 - SAVE and exit data entry" 


t 


END SUB 

TSEC CSCC ST SCC CSCS TC CSCCSCOCS CSCIC CSCC SC. Tr 
kx Name: TeleAlign * 
kx Type: Subroutine ** 
1 kx Module: MIRACTRL. BAS kx 
kx Language: Quickbasic 4.50 * 


PKK KKKKKKKKEKKKKKKEKEKEKKKEKKKEKEKKKKKKKKKEKRKEER 

1 

'Allows user input of the telescope's Hour Angle and 
Declination to align 

'the telescopein the local coordinate system. 

1 

'Example of use: TeleAlign GMSTOhrSEC#, GHAaries#, TeleRA#, 
TeleDEC# 

SUB TeleAlign (GMSTOhrSEC#, GHAaries#, TeleRA#, TeleDEC#) 

| 


‘initialize Variables 


MITRAlatitude# = .63908139# ‘radians, (36 deg 37 min 
North) 

MIRAlongitude# = 2.12639281# 'radians, (121 deg 50 min 
West) 


sec% = VAL(RIGHTS(TIMES, 2)) 
min% = VAL(MIDS(TIMES, 4, 2)) 
hour% = VAL(LEFTS(TIMES, 2)) 
t 
‘Calculate Greenwhich Mean Sidereal Time 
'Determine the day fraction that has passed since midnight 
'in SECONDS 


sec& = sec% 
heurss& = hours 
mins& = min% 


dayFraction# = CDBL((hours& * 3600) + (mins& * 60) + sec&) 
| 


‘Determine the Greenwhich Mean Sidereal Time at the current 
'time in SECONDS 

GMSTinSEC# = GMSTOhrSEC# + dayFraction# 

| 

'Calculate Greenwhich Hour Angle of Aries 
GHAaries# = GMSTinSEC# * 2# * pi / 86164.09053099999# 'Seconds 
of solar time converted to radians 

DO WHILE GHAaries# > (2# * pi) 

GHAaries# = GHAaries# - (2# * pi) 
LOOP 
DO WHILE GHAaries# < 0O# 
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GHAaries# = GHAaries# + (2# * pi) 
LOOP 
i 
‘Calculate the telescope's hour angle for this moment 
RAofMIRA# = GHAaries# - MIRAlongitude# 


$ 
PK HKKKKKKKKKEKKKEKEKKKEKKKKEKKEKEKEKEKEKEKKEKKKKKKEKEKEKEKKKKEKE 


'Read telescope's Hour Angle from the encoders ** 
'Store to variable HourAngle# in radians * 
PKK KKH KKK KKKKKKKKKHKKKKKEKEKKKKKKEKKEKKEKKEKKEKEKEKKEKKKKEKSE 
TeleRA# = RAofMIRA# - HourAngle# 'in radians 
DO WHILE TeleRA# > (2# * pi) 

TeleRA# = TeleRA# - (2# * pi) 

LOOP 
DO WHILE TeleRA# < O# 

TeleRA# = TeleRA# + (2# * pi) 

LOOP 


§ 
PK KKKKKKKKEKKKEKKKKEKEKKKEKEKKEKEKEKKEKEKKEKEKEKEKKEKEKKEKKKKKEE 


'Read telescope's Declination from the encoders ** 


'Store to variable TeleDEC# in radians x 
eee eee ek kkk ke KKK KKK KKEKKR KKK KK KERKKKKERER KS 


END SUB 

eee c reece eee eS eee ee eS eee ee eee ee ee 
'kx* Name: TeleAngles * 
tak Type: Subroutine * 
kx Module: MIRACTRL. BAS xx 
kx Language: Quickbasic 4.50 * 


VK KK KKK KKK KKK KKK KKK KK KKK KKK KKK KKK KEK KKK 
$ 

‘Calculates the telescope's zenith and azimuth angles in the 
‘terrestrial coordinate frame. 

‘ 

'Example of use: TeleAngles Zenith%, Azimuth%, GHAaries#, 
RAofStar#, DECofStar#, CelCoord#(), DCM#(), TerraRec#() 

i] 


SUB TeleAngles (Zenith%, Azimuth%, GHAaries#, RAofStar#, 
DECofStar#, CelCoord#(), DCM#(), TerraRec#() ) 


‘'Tnitialize variables 


MIRAlatitude# = .63908139# ‘radians, (36 deg 37 min 
North) 
MIRAlongitude# = 2.12639281# 'radians, (121 deg 50 min 


West) 
§ 


‘Define a celestial coordinate vector in RECTANGULAR 
coordinates 

CelCoord# (1) 
CelCoord# (2) 


COS (RAofStar#) * COS(DECofStar#) 
SIN(RAofStar#) * COS(DECofStar#) 
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CelCoord#(3) = SIN(DECofStar#) 
1 


'Define the direction cosine matrix which relates celestial 


and terrestrial vecormdinates 


DCM#(1, 1) = (-SIN(GHAaries#) * SIN(MIRAlongitude#) 
SIN(MIRAlatitude#) - COS(GHAaries#) * COS(MIRAlongitude#) 
SIN(MIRAlatitude}#) ) 

DCM#(1, 2) = (-SIN(GHAaries#) * COS(MIRAlongitude#) 


SIN (MIRAlatitude#) + COS(GHAaries#) * SIN(MIRAlongitude?#) 
SIN (MIRAlatitude#) ) 
DCM#(1, 3) = COS(MIRAlatitude#) 


DCM#(2, 1) = (-COS(GHAaries#) * SIN(MIRAlongitude}#) 
SIN(GHAaries#) * COS(MIRAlongitude#) ) 
DCM#(2, 2) = (-SIN(GHAaries#) * SIN(MIRAlongitude#) 


COS (GHAaries#) * COS(MIRAlongitude?#) ) 

DCM#(2, 3) = 0 

DCM# (3, 1) = (COS (GHAaries#) * COS (MIRAlongitude#) 
COS (MIRAlatitude#) + SIN(GHAaries#) * SIN(MIRAlongitude}#) 
COS (MIRAlatitude#) ) 

DCM#(3, 2) = £(SIN(GHAaries#) * COS(MIRAlongitude?#) 
COS (MIRAlatitude#) - COS(GHAaries#) * SIN(MIRAlongitude#) 
COS (MIRAlatitude#) ) 

DCM#(3, 3) = SIN(MIRAlatitude}#) 

t 


‘Calculate the terrestrial rectangular coordinates 


TerraRec#(1) = DCM#(1, 1) * CelCoord#(1) + DCM#(1, 2) 
CelCoord#(2) + DCM#(1, 3) * CelCoord# (3) 
TerraRec#(2) = DCM#(2, 1) * CelCoord#(1) + DCM#(2, 2) 
CeliGoord#(2) + DCMi(2, 3) * Cellecoara 13) 
TerraRec#(3) = DCM#(3, 1) * CelCoord#(1) + DCM#(3, 2) 


CelCoord#(2) + DCM#(3, 3) * CelCoord# (3) 


'Calculate the azimuth and zenith angles 
IF TerraRec#(1) < 0O# AND TerraRec#(2) < O# THEN 
Azimuth# = ATN(TerraRec#(1) / TerraRec#(2)) + (pi 
2 #) 
ELSEIF TerraRec#(1) >= O# AND TerraRec#(2) < O# THEN 
Azimuth# = ATN(-TerraRec#(2) / TerraRec#(1) ) 
ELSEIF TerraRec#(1) < O# AND TerraRec#(2) >= O# THEN 
Azimuth# = ATN(-TerraRec#(2) / TerraRec#(1)) + pi 


+ 


% 


ELSE Azimuth# = ATN(TerraRec#(1) / TerraRec#(2)) + ((3# * pl) 


J 2#) 
END IF 
Azimuths = INT(AzZimuehn+ +6960 F 97 a 27 ne 
DO WHILE Azimuth% < 0O 
Azimuth% = Azimuth% + 360 
LOOP 
DO WHILE Azimuth% > 360 
Azimuth% = Azimuth% - 360 
LOOP 
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Elevation# = ATN(TerraRec#(3) / SQR(1 - TerraRec#(3) * 2#)) 
Zenith% = 90 - INT(Elevation# * 360# / (2# * pi)) 


END SUB 

tekekkkkekkkekekkekekkekkekeekeekekeKRKRKRKRKRKRKR KKK 
tk* Name: TeleStatusDisplay ** 
tex Type: Subroutine kx 
kx Module: MIRACTRL. BAS xx 
kx Language: Quickbasic 4.50 * 


PK KKEKEKKKKEKEKKKEKEKEKEKKKEKKKKKEKKKKK KKK 
t 


'Creates the display screen for various elements of telescope 
‘status and updates them while the telescope is tracking or 
slewing. 

] 


'Example of use: TeleStatusDisplay ActiveStar$ 

‘Parameters: (none) 

'Variables: Star$() A 5-element matrix of desired 
stars for 

: viewing 

RAS () A 5-element matrix of right 
ascensions 

: for the stars listed in Star$() 
DECS () A 5-element matrix of declinations 
: for the stars listed in Star$() 

: ActiveStar$ A string naming the star for 
which a tracking 

; solution is currently being 


calculated 

5 Pe ak Looping indices 

‘Module Level 

: DECLARATIONS: DECLARE SUB TeleStatusDisplay 


(ActiveStar$) 
| 


SUB TeleStatusDisplay (ActiveStar$, Simulations) 
Crs 
‘Color the background 
FOR 1% = 1 TO 25 
COLOR BLUE, BLACK, BLUE 
PRINT STRINGS(80, 177) 
NEXT 1% 
MOCATE 24, 1, O 
PRINT STRINGS(80, 177); 
LOCATE 25, 1, 0 
PRINT STRINGS(80, 177); 
g 


‘Draw top edge of star name box 

COLOR WHITE, BLUE, BLUE 

ROCATE 10, 29, O 

PRINT CHRS$(201); STRINGS(19, 205); CHRS$(187); 
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COLOR WHITE, BLUE, BLUE 

LOCATE 10Ope > eu 

PRINTS (Siva! 

t 

'Draw the body of the star name box and place the name 
COLOR WHITE, BLUE, BLUE 

LOCATE, Lie 2o, © 

PRINT CHRS$(186); STRINGS$(19, 0); CHRS$(186); 
LOCATE 11, INT((80 — LEN (ACElVeStars) a7 ez) no 
PRINT ActiveStar$ 

t 

'Draw bottom edge of the star name box 

COLOR WHITE, BLUE, BLUE 

LOCATE .i2Zy29, 0 

PRINT CHRS (200) ; SPTRINGS{( 19> 205) 7 -CHRo (ho cr, 
t 


'Place the lettering 

COLOR WHITE, RED, BLUE 

LOGALEM2 27,0 

PRINT " ESC key halts telescope " 
COLOR WHITE, BLUE, BLUE 

LOCATE) 6728, 6 

PRINT "Oliver Observing Station" 
LOCATE 7, 32, 0 

PRINT "Telescope Status" 


‘Draw top edge of status box 

EOGATE  1375.5,..0 

COLOR WHITE, BRUE.. BLUE 

PRINT CHRS(201)G 7StRINGS (6 Cpee2 Os Chi onens 7 a; 
t 


‘Draw the body of the status box 
FOR r% = 14 TO 21 
LOCGATETE 4705570 
PRINT CHRS(186); 
PRINT SPACES (68) ; 
PRINT CHR (lea ye; 
NEXT r% 
{ 
‘Draw bettom edge of (Ene stacus box 
ECCATEVZ2 9.57 20 
PRINT CHRS(200); STRINGS(68)> 205) >) CHR ene 250; 
LOGE 2 2 7.82 370 
PRINT " ENTER to change telescope status " 
t 


'Place the help instructions 
LOCATE 2307 .culam0) 

PRINT "Fl for help screen" 
'Place the lettering 

COLOR WHITE, BLUE, BLUE 
LOCATE atdye sO e0 


iL's).0, 


BRENT RA" 


LOCATE 
PRINT " 
LOCATE 
PRINT " 
LOCATE 
PRINT “ 
LOCATE 
Pern? ™ 
LOCATE 
PeeNT ™ 
LOCATE 
ERINT " 
LOCATE 


ile ees a ee) 
DEC" 

XG, Me, ow 
Hour Angle" 
14, Sages 
zenith Angle" 
a SB Ore U 
Azimuth" 

14, 58, 0 

PST" 

ik 7) 8 sis amee 
Sidereal Time" 
20, 58, 0 


EraNnr "UTC" 


'Show motor state boxes 
COLOR RED, BLUE, BLUE 


LOCATE 
PRINT " 


1335, 0 
Simulation" 


COLOR WHITE, BLUE, BLUE 


LOCATE 
SRiNT " 
LOCATE 


ee 0 
RA Slew" 
5 lz 0 


PRINT “RA Step" 


LOCATE 
ER@NT ™ 
LOCATE 
EraNT “ 
FOR 1% 


NEXT 1% 
FOR 1% 


NEXT 1% 
t 
END SUB 
VKKKKKKK 


Ue x 
tkx* 


Dy D7 0 

DEC Slew" 

Sno 6 70 

DEC Seep” 

= i ToO°2 

LOCATE 2 uta. sa) — 17 6 

PRENT Chbsi20l) > CHRS (205) ; CHRS(187) 
LOCATE 3, 10 * 1% - 1, O 

PRIME CHRs (186)> SPACES(1); CHRS(186) 
EOGATE 4, Os 1% = 7159 0 

PRINT CHRS (200); CHRS (205); CHRS ss) 


= 1 TO 2 

LOCATE 2, (10 * iS - 1) + 50, 0 
PRINT CHR$(201); CHRS$(205); CHR$(187) 
LOCATE 3, (10 * i - 1) + 50, 0 
PRINT CHRS(186); SPACE$(1); CHR$(186) 
LOCATE 4, (10 * i - 1) + 50, 0 
PRINT CHRS(200); CHRS$(205); CHRS$(188) 


KEKE KKK KKK KKK KEKE KHE 
Name: TrackCommands K* 
Type: Subroutine * * 


boat 


kk Module: MIRACTRL. BAS bal 


tx Language: Quickbasic 4.50 * 
PKK KKKKEKKKEKEKEKKEKEKKKKEKEKKEKKKEKKEKKKKEKKKEKKEEEE 
| 


‘Controls the right ascension stepping motor only to 
compensate 

'for local coordinate drift due to the rotation of the Earth. 
i 


'Example of use: TrackCommands 
| 


SUB TrackCommands 
| 


'Track as appropriate in Right Ascension only 
'Update motor display box 

COLOR RED, BLUE, BLUE 

LOCATE 3, 20, 0 

PRINT CHRS$ (222) 

'POKE (address of RA track controller), X 


END SUB 

PKK KKEKKKKEKKKKKKEKKEKK KKK KKK KKK KEK 
"kx Name: VideoState * 
kx Type: Subprogram * x 
kx Module: BIOSCALL. BAS x 
kk Language: Quickbasic 4.50 * 


VK KKKKKEKEKKEKEKEKEKKEKEKEKKKEKEKKKRKKKKKKKKKKKKKEEEE 
| 


‘Determines the current video mode parameters. 
i] 


'Example of use: VideoState mode%, columns%, page% 


S RegType) 
DECLARE SUB VideoState (mode%, columns%, page%) 


'Parameters: mode% Current video mode 

columns% Current number of text columns 
page% Current active display page 
‘Variables: reg Structure of type RegType 
‘MODULE LEVEL 

DECLARATIONS: TYPE RegType 

: ax AS INTEGER 

: bx AS INTEGER 

cx AS INTEGER 

: ax AS INTEGER 

' Bp AS INTEGER 

: Si AS INTEGER 

: di AS INTEGER 

' flags AS INTEGER 

END fy 

| 

: DECLARE SUB Interrupt (intnum%, inreg AS RegType, outreg 
. 

| 


SUB VideoState (mode%, columns%, page%) STATIC 
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AND &HFFOO) \ 256 
D &HFFOO) \ 256 


DIM reg AS RegType 

reg.ax = &HFOO 

Interrupt &H10, reg, reg 

mode% = reg.ax AND &HFF 

columns% = (CLNG(reg.ax) 

pages = (CLNG(reg.bx) AN 
END SUB 
PK KKKKKKKKEKKEKKEKKEKKEKEKKKEKKEKKEKKEKEKEKKEKKEKKEEEKSE 
kx Name: Windows 
'kx* Type: Subprogr 
kx Module: WINDOWS. 
Vek Language: Quickbas 


xx* 
am xx 
BAS x 
eu. 50) xx 


PKKKKEKKKKKEEKKKEKKKEKKEKKKEKKEKEKKEKKEKEKKKKEKKKEKEEKEE 


'Displays a rectangular window for information display 


'or menu s 
t 


'Example o 
t 


Creer lon. 


f use: Windows wl, 


‘Parameters: wl 

; wText$ () 
' wTitles 
: wPrompts 
| 

'Variables: mode% 

' columns% 
columns 

pages 

' cursorRow% 
' cursorcols 
: newpage% 
: lbText% 
lines 

' ubText% 
lines 

' nes 

' maxlen% 
display 

length% 
' row2% 
corner of window 

: EGo612s 
corner of 

| 

f ulZ 
character code 

: urs 
character code 

; 11% 
character code 


wText$(), wTitle$, wPrompts$ 
Structure of type WindowsType 
Array of strings to be displayed 
Pipke. string 

Prompt string 


Current video mode 
Current number of character 


Current video page 
Saved cursor row position 
Saved cursor column position 
Next video page 
Lower boundary of array of text 
Upper boundary of array of text 


Looping index 
Length of longest string to 


Length of each array string 
Row number at bottom right 


Column number at bottom right 


window 
Upper left corner border 


Upper right corner border 


Lower left corner border 


LoS 


eli 


! 
co 
t 
| 
ea 
t 
| 
fu 
'M 
t 
wT 
| 
| 


SU 
Se 
' 


Nes Lower right corner border 


aracter code 
v1s Vertical border character code 
hl% Horizontal border character 
de 
r% index to each line of text 
EGhneoices Set of unique characters for 
ch menu line 
tmp$ Work string 
kees Key code returned by InKeyCode% 
Netaon 
ODULE LEVEL 


DECLARATIONS: SUB Windows (w AS WindowsType, 
itles, 
wPrompt$) STATIC 


B windows (w AS WindowsType, wText$(), wTitles, 
ATIC 

Key code numbers 

CONST DOWNARROW = 20480 

CONST OEN@TER — 13 

CONST ESCAPE = 27 

CONST UPARROW = 18432 


Determine current video page 
VideoState mode%, columns%, page% 


Record current cursor location 
cursorRow% = CSRLIN 
cursorCol% = POS(0) 


Window will be on the next page, if available 


wText$S(), 


wPrompt$) 


newpage% = page% + 1 

IF newpage%t > 7 THEN 
SCREEN 207 0 
PRINT "Error: Windows - not enough video pages" 
SYSTEM 

END IF 


Copy current page to new page 
PCOPY page%, newpage% 


Show the current page while building window on new page 


SCREEN , , newpage%, page% 


Determine array bounds 
lbText% = LBOUND(wText$) 
ubText% = UBOUND(wText$) 


Check the text array bounds, lower always 1, upper > 0O 


IF lbText% <> 1 OR ubText% < 1 THEN 
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SCREEN , , 0, 0 
PRINT "Error: Windows - text array 


imcorrec tly” 


eyes ay 
END 1: 


Determine longest string in text array 
maxLen% = O 
FOR 1% = lbText% TO ubText% 

length% = LEN(wTexts(i%) ) 

IF length% > maxLen% THEN 

maxLen% = length% 

END IF 

NEXT 1% 


Determine the bottom right corner of window 
row2% = w.row + ubText% + 1 
Gol2% w.col + maxLen% + 3 


Gaeck that window fits on screen 


dimensioned 


me w.row < 1 OR w.col < 1 OR row23% > 25 OR col2% > columns? 


THEN 


Sere EEN... 0, 20 
PRINT "Error: Windows - part of window is off screen" 
PRINT "columns% = "; columns% 
SYSTEM 
BND IF 


Set the edge characters 
SELECT CASE w.edgeline 


CASE O 
ul% = 32 
ur% = 32 
11% = 32 
1lr% = 32 
vl% = 32 
hls = 32 

CASE 1 
ul% = 218 
ur% = 191 
11% = 192 
1r% = 217 
vl% = 179 
his = 196 

CASE 2 
ul% = 201 
ur% = 187 
11% = 200 
1r% = 188 
v1% = 186 
hl3 = 205 


Lo 


CASE ELSE 


SCREENS GaGa & 
PRINT “Error: Windows - Edge line types incorrect" 
SYSTEM 


END SELECT 
i] 
'Draw top edge of box 
LOCATE wW.TGOWw, w.col, 0 
COLOR w.fgqdEdge, w.bgqdEdge 
PRINT CHRS$(ul%); STRINGS (maxLen% + 2, h1%); CHRS$(ur%) ; 
t 


'Draw the body of the window 

FOR r% = w.row + 1 TO row2% - 1 
LOCATE r%, w.col, 0O 
COLOR w.fgdEdge, w.bgqdEdge 
PRINT CHRS(v13%); 
COLOR w.fgdBody, w.bgdBody 


tmp$ = LEFTS(wText$(r% - w.row) + SPACES(maxLen%), 
maxLen$) 
PRINT Th W ; tmps$; Th th 


COLOR w.fgdEdge, w.bgdEdge 
PRINT CHRS$(v1%); 
NEXT r% 
t 
'Draw bottom edge of the box 
LOCATE row2%, w.col, 0 
COLOR w.fgdEdge, w.bgdEdge 
PRINT CHRS$(11%); STRINGS (maxLen% + 2, hl1%); CHRS$(1r%) ; 
i] 
'Center and print top title, if present 
IF wlieles <> 2S THEn 
LOCATE w.row, (W.col + col2% - LEN(wTitlesS) + 1) \ 2, O 
COLOR w.fgdTitle, w.bgdTitle 
PRINT wTitles; 
END IF 
t 
'Center and print prompt, if present 
fF wPrompts <> "" THEN 
LOCATE row2%, (w.col + col2% - LEN(wPrompts) + 1) We2Z7e 
COLOR w.fgdPrompt, w.bgdPrompt 
PRINT wPrompts; 
END IF 
t 
'Now make the new page visible and active 
SCREEN , , newpage%, newpage% 
| 
'Take next action based on action code 
SELECT CASE w.action 
CASE 1 
'Get a key code number and return it 
DO 


ILS 3S) 


w.returnCode = KeyCode%(Character$) 

LOOP UNTIL w.returnCode 

1 

CASE ELSE 
w.returnCode = 0 

END SELECT 

1 

'Reset the cursor position 


LOCATE cursorRow%, cursorcol% 
t 


END SUB 

PKKKKEKKEKKKKKKEKKKKKEKEKKKKEKKEKEKKEKKEKEKKKKEKEKEEKSE 
Vex Name: WindowsPop * 
"kx Type: Subprogram xx 
'k* Module: WINDOWS. BAS x 
kx Language: Quickbasic 4.50 * x 


PHKKKEKEKEKKKKEKEKEKEKKKKEKEKEKEEKEKKKKKKKKEKKKEKEKKEKKEEE 
i 


'Removes last displayed window 
1 


'Example of use:WindowsPop 


‘Parameters: (none) 

'Variables: mode% Current video mode 

columns% Current number of display columns 
pages Current display page 


edule level 
' DECLARATIONS: DECLARE SUB WindowsPop 
i] 
SUB WindowsPop STATIC 
VideoState mode%, columns%, page% 
IF page% THEN 
SCREEN 0, , pages - 1, pages - l 
END IF 
END SUB 


ow, 


B. SAMPLE DATA FILE "STARFILE. DAT!" 


"Sirius! 3"23:01.40", "O07O0: 0G" ooo. 
"Betelquese™, "03:05:00", "00.00.00. U2000" 
“Pollux Il 233421 2334. Od On 
"Dubhe" "250334259" , "56.34.59 200°C! 
"Castor" 14) POO US S37 PO.) a OCs 


f 

eee | 6tt ee Re Oe te 
id / I 

ete) | 6k ee Oh ee et 
/ Uy f 

rt) | 6ftet oat os ones 
I f f 

ont 6688 otk sets 
Uy / f 

ret | 6tt te ke lee ae 
f , f 

rte) | 60het | 6th ee Ot eS 
/ a id 

rete | 6th e Oke Ot ae 
f f f 

teet | offer 6k ee oft tt 
f / f 

eet 6th et 68k tt OSES 
f f I 

rrr | 6th te Othe One 
I f f 

mee 6fe te ose tk Ont at 
, f f 

tert 6th te oft te ote ae 
I a I 
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